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ABSTRACT 

The  stress  level  in  the  milled  skin  at  the  root  of 
a  tapered,  multicell  swept  wing  is  predicted  by  means  of 
a  matrix-force  method.   The  solution  is  achieved  by- 
minimizing  the  internal  strain  energy,  and  the  results 
compared  with  experimental  tests.   A  single  loading, 
consisting  of  a  nose-up  couple,  is  applied  to  each  tip  rib. 
The  loading  is  transferred  from  the  tip  rib  to  an  idealized 
structure  by  means  of  simple  torsion  theory.  The  idealized 
structure  represents  the  inboard  one-half  of  each  semi-span. 

Results  indicate  that  an  accurate  solution  of  the 
stress  distribution  in  the  actual  wing  can  be  achieved 
provided  that  the  root  boundary  conditions  are  preserved. 
An  extension  of  the  analysis  is  suggested  in  order  to 
more  closely  define  the  maximum  accuracy  inherent  in  the 
particular  matrix-force  method  of  solution. 
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INTRODUCTION 

Highly  swept  wings  of  low  thickness  to  chord  ratios 
were  first  built  with  conventional  stringer  sheet  methods. 
This  introduced  the  problem  of  how  to  successfully  predict 
and  then  to  relieve  the  severe  concentration  of  stress  along 
the  rear  spar  and  aft  root  section  skin. 

The  stress  concentration  was  reduced  somewhat  by  shift- 
ing to  multi-cell  construction.  Multi-cell  construction  as 
used  herein  describes  wings  having  closely  spaced  spanwise 
spars  or  webs,  and  relatively  few  streamwise  ribs. 

Wing  analysis  b3r  customary  beam  bending  and  torsion 
methods  was  unable  to  accurately  account  for  the  root 
area  stress  levels  unless  modification  factors  were  used. 
This  was  due  primarily  to  the  neglect  of  bending-torsion 
inter-action.   However,  as  flight  speeds  increased  and 
aspect  ratios  decreased,  the  inability  of  beam  theory  to 
properly  define  chordwise  deformations  caused  even  greater 
concern  since  thermo-aeroelastic  problems  replaced  static 
stressing  in  order  of  importance. 

In  order  to  properly  account  for  bending-torsion 
interaction,  chordwise  curvatures,  shear  deflections,  the 
increasing  use  of  large  cut-outs,  and  very  heavy  milled 
skin,  several  alternative  methods  of  analysis  have  been 
developed.  Without  considering  the  older  relaxation 
techniques  and  their  iterative  solutions,  each  depends 


upon  the  use  of  a  digital  computer  for  solution.   Each 
attempts  to  achieve  a  mathematical  model  which  exhibits 
a  properly  deformed  shape  having  absolute  displacements 
comparable  to  the  actual  case. 

The  most  prevalent  of  these  methods  treat  the  struct- 
ure as  an  assemblage  of  elastic  components,  which  allows 
matrix  formulation  of  the  solution  in  terms  of  the  various 
energy  theorems.   Argyris  (Ref.  1)  has  shown  that  these 
energy  theorems  derive  from  the  two  fundamental  principles 
of  virtual  displacements  and  virtual  forces,  both  of  which 
originated  in  the  work  of  Maxwell,  Mohr,  and  Engesser. 
Since  both  principles  are  independent  of  elastic  laws  or 
the  structural  material,  their  application  to  non-linear 
problems  or  to  structures  where  initial  thermal  strains 
exist  is  quite  easy. 

The  "matrix-force"  method  leads  to  an  analysis  in 
terms  of  forces  as  the  unknowns,  while  the  "matrix- 
displacement"  method  leads  to  an  analysis  in  terms  of 
displacements.  Either  method  provides  for  exact  satis- 
faction of  equilibrium  and  compatibility  in  a  structure 
made  up  of  discrete  elements.   Simplification  of  structural 
modifications  during  the  preliminary  design  stage  is  an 
additional  advantage  which  accrues  from  the  use  of  discrete 
elements. 


Each  method  is  again  categorized  as  to  whether 
:f lumped-parameter7'  or  "finite  element"  idealization  of 
plane  panels  is  used.   "Lumped-parameter"  idealization 
as  used  herein  allows  distinct  properties  of  the  structure, 
such  as  direct  stress,  shear,  bending  and  torsion,  to  be 
separated  or  concentrated  at  discrete  locations.   The  web 
or  panel  may  be  attached  to  the  adjacent  flanges  at 
corners,  mid-points  of  the  panel,  or  continuously.   "Finite- 
element"  idealization  prescribes  plane  panels  of  either 
triangular  or  rectangular  shape,  to  be  attached  at  corners 
or  nodes  of  the  structure.   In  either  idealization,  panel 
warping  must  be  avoided,  even  by  slightly  revising  the 
geometry  of  the  structure,  in  order  to  achieve  a  rigorous 
solution. 

Although  a  complete  solution,  giving  both  forces  and 
displacements,  may  be  obtained,  from  either  the  matrix-force 
or  matrix-displacement  method,  it  has  been  shown  in  Ref.  2 
that  if  only  forces  are  of  interest,  the  matrix-force 
method  is  the  most  efficient.   Conversely,  if  only 
deflections  are  of  interest,  the  matrix-deflection  method 
is  the  most  efficient. 

The  degree  of  redundancy  in  the  force  method  is  less 
than  for  the  displacement  method.   However,  the  stiffnesses 
of  the  elements,  as  used  by  the  displacement  approach,  are 
easier  to  obtain  than  are  the  flexibilities  of  elements 


required  by  the  force  approach.   The  force  method  has 
the  additional  disadvantage  of  neglecting  Poissonfs 
effect  for  that  portion  of  the  axial-load  carrying  skin 
which  is  luraced  with  the  beam  flanges. 

Argyris  (Ref.  3)  has  shown  that  while  the  force 
method  is  nearly  always  more  suitable  for  fuselage  analysis, 
the  displacement  method  may  often  be  more  convenient  for 
the  analysis  of  complex  win,p;s,  since  in  the  displacement 
method  structural  "node  points"  specify  the  matrix  size. 
This  generally  results  in  smaller  matrices  than  would  be 
required  to  obtain  similar  accuracy  by  the  force  method. 

Current  literature  abounds  with  instances  where  the 
two  methods,  each  formulated  several  ways,  have  been 
applied  to  laboratory  specimens.   This  allows  some  study 
of  parameter  effects,  but  allows  no  true  indication  of  the 
method  applicability  to  an  actual  wing  which  may  possess 
extreme  and  sometimes  abrupt  variations  in  skin  thickness, 
discontinuous  spars,  massive  ribs,  skin  and  web  cut-outs, 
and  taper. 

It  is  the  object  of  this  report  to  compare  predicted 
stress  and  strain  levels  near  the  root  of  an  actual  wing 
(which  possesses  all  of  the  aforementioned  analytical 
obstacles)  with  the  stresses  and  strains  obtained  experi- 
mentally.  The  specimen  to  be  analyzed  was  obtained  from  a 
Mach  2.5  all-weather  fighter  project  which  was  canceled 


by  the  U.  S.  Navy  before  the  structural  test  program  was 
completed.   The  matrix-force  method  of  analysis  is  used, 
and  consists  of  a  "lumped-parameter"  structural  idealization 
as  proposed  by  Wehle  and  Lansing  (Ref.  4).  The  matrix 
formulation  of  the  solution  is  similar  to  the  scheme 
proposed  by  Lang  and  Bisplinghoff  (Ref.  5)  and  organized 
in  considerable  detail  in  Ref,  6.  This  formulation  allows 
efficient  use  of  the  digital  computer  and  at  the  same  time 
allows  intermediate  results  to  be  printed  out  for  study. 
This  characteristic  makes  changes  to  the  structural  ideal- 
ization relatively  simple  to  incorporate,  and  causes  the 
effect  of  such  a  change  to  be  readily  apparent.  Thus  from 
the  outset  it  was  recognized  that  a  follow-on  analysis, 
using  the  same  digital  computer  program,  would  allow 
investigation  of  the  boundary  conditions,  and  indicate 
the  merit  of  the  idealization  with  only  limited  additional 
effort. 

This  work  was  conducted  at  the  U.  S.  Naval  Postgraduate 
School,  Monterey,  California,  during  the  1962-1963  school 
year.   The  authors  are  deeply  indebted  to  Professor  C.  H. 
Kahr,  who  conceived  the  need  for  such  an  investigation. 
Professor  Kahr,  together  with  Associate  Professor  U.  Haupt, 
rendered  invaluable  assistance  and  encouragement  throughout 
the  project.   The  authors  express  their  gratitude  to 
Aeronautical  Laboratory  Supervisor  R.  E.  McConnelland 


assistants  R.  A.  Besel  and  R.  0.  Cunningham,  whose 
cooperative  efforts  were  so  essential  in  the  procurement, 
fabrication  and  installation  of  all  experimental  apparatus. 
Appreciation  is  also  expressed  to  Chance  Vought  Aircraft, 
Inc,  for  their  cooperation  in  furnishing  essential  data. 


II   THEORETICAL  ANALYSIS 


1.   General 


The  structure  chosen  for  the  analysis  was  the  wing 
center  section  of  an  F&U-3,  as  shown  in  Fig.  1  below.  This 
section  consists  of  seven  spars  and  one  post  beam  per  semi- 
span,  upper  and  lower  milled  skins,  two  massive  pivot  ribs 
closely  spaced  about  the  aircraft  center  line,  two  inter- 
mediate ribs,  and  two  heavy  wing  fold  ribs  located  at  the 
extremities  of  the  center  section. 


Pivot  Ribs 


Post  Beam 


Intermediate 
Rib 


Fold 

Line  I 


Fig.  1 
Sketch  of  FPU-3  Wing  Center  Section 

The  first  requirement  for  the  analysis  was  to  study 
the  pictures  and  drawings  of  the  dis-assembled  wing  in 
order  to  establish  a  basis  for  simplifying  assumptions. 
These  simplifications  would  be  necessary  to  enable 
representation  of  the  structure  by  a  mathematical  model 


consisting  only  of  axial-load  carrying  bars  and  plane, 
shear  carrying,  constant  thickness  panels. 

The  first  simplification  was  introduced  by  considering 
the  pivot  rib  to  be  the  root  location.   The  pivot  rib  was 
thus  assumed  to  be  infinitely  stiff  and  free  of  chordwise 
bending.  This  step  was  taken  to  simplify  the  mathematical 
model,  reducing  by  16  the  number  of  components  to  be 
analyzed.   This  xvould  undoubtedly  reduce  the  accuracy  of 
the  solution.   However,  the  effect  should  be  predictable, 
and  once  the  reduced  solution  is  achieved,  incorporation  of 
the  center  section  (pivot  rib  to  aircraft  center  line) 
would  require  comparatively  minor  additional  effort. 

A  post  beam,  shown  in  Fig.  1,  extended  from  the  pivot 
rib  ("root")  to  the  intermediate  rib  in  the  actual  wing. 
Since  no  shear  web  was  involved,  the  post  beam  was  removed 
from  the  idealized  structure  and  its  flange  area  distributed 
to  the  surrounding  structure.  The  posts  themselves  were 
ignored.  The  remaining  structure  was  then  represented  by 
axially  loaded  bars  and  flat  shear  panels  of  constant 
thickness,  as  shown  in  Fig.  2. 

It  was  then  necessary  to  determine  the  cross-sectional 
areas  of  the  bars  and  the  thickness  of  the  shear  panels. 
This  was  done  for  the  bars  at  the  cross-sections  located 
at  stations  yw  =  25.095  and  yw  =  £1.93,  Fig.  3.   A  bar  end 
area  included  the  existing  flange  area,  and  the  total  skin 
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area  perpend icular  to  the  bar,  taken  from  the  center  of 

the  panel  just  ahead  to  the  center  of  the  panel  just  behind, 

as  shown  in  Fig.  4.   Spar  web  areas  were  not  included. 

The  concentration  of  axial  load  carrying  skin  at  the  flange 

locations  has  the  effect  of  removing  lateral  contraction 

of  the  covers  from  much  of  the  analysis  (  jj.   =  0). 


Beams 


y  =  81.98 


y-  yw=  25,095 


Fig.  4 
Distribution  of  Cover  Skin  to  Beam  Flanges 

Difficulty  was  experienced  in  determining  the  "average" 
skin  thickness  of  each  of  the  14  cover  panels  due  to  the 
span-wise  variation  of  thickness.   Chordwise  variation  in 
thickness  was  easier  to  account  for  since  the  aspect  ratio 
of  the  individual  panels  was  on  the  order  of  four. 
Representation  of  the  multiply  tapered  skin  by  panels  of 
average  thickness  will  undoubtedly  penalize  the  accuracy 
of  the  analysis. 
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The  un-tapered  bars  of  the  intermediate  rib  were 
formed  by  adding  to  the  existing  chordwise  flanges.   The 
added  area  consisted  of  a  skin  area  equal  to  the  local 
thickness  times  a  plate  width  (equal  to  the  panel  width) 
distributed  evenly  on  either  side  of  the  rib,  as  shown 
in  Fig.  5. 


rib  augmented  by  this  skin 


Fig.  5 
Distribution  of  Cover  Skin  Area  to  Rib  Bar  Area 

The  dimensions  of  the  resulting  elements  of  the 
idealized  wing  are  listed  in  Table  I. 

2.   Loading  Method 

After  establishing  the  idealized  model,  it  was 
necessary  to  establish  a  means  of  load  application.   It 
must  transfer  a  load  identical  to  the  tip-applied  load  on 
the  actual  wing  to  the  idealized  model  at  the  intermediate 
rib.   As  shown  in  Fig.  6,  a  nose-up  (the  wing  is  mounted 
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Fig.  6 


Application  of  Couple  to  Wing  Tip 


13 


CD 

O 

+> 

cc 

to 

*a 

li 

CD 

111*                                  \    \ 

* 

*•     ,  > 

•                                       ** 

v 

•      o  »a 
to           o> 

0L         / 
CD                          J 

$ 

Cu         O  -H 

—                       / 

\  y 

buti 
deal 

V 

•HM 

•   #     \  \ 

b 

■p 

CO 
•H 

Q 

O    /T    \\             \ 

\      ~/\    / 

80                    1          \\                \ 

\               Jm     l\       ' 

1  J\  w 

\\5*®\ciie  i 

^M<»> 

V*.     /"\\          \\     o- 

QcvJ 

CQ 

•     *-s         \\          \\       A 

v4/ 

/®  y\ 

w   / 

\  \J\  M 

yA 

\  tt    7 

\  \\/  \  /®  y 

\\A  -  /  X 

^ 

of" 

• 

* 

ii       ii 

• 

1tf 


inverted)  couple  perpendicular  to  the  C.I.B.  was  applied 
to  each  wing  tip. 

Considering  only  the  instrumented  wing,  this  couple 
must  be  converted  to  12  cover  shears ,  seven  vertical  web 
shears,  and  14  axial  flange  loads  on  a  streamwise  cross- 
section  (intermediate  rib)  at  station  y  s  £l.9#,  as  shown 
in  Fig.  7.  This  was  accomplished  by  assuming  simple  torsion 
theory  would  be  satisfactory  inboard  from  the  tip  to 
approximately  the  intermediate  rib.  This  assumption  was 
based  upon  tests  reported  in  Ref.  7. 

The  shear  flows  induced  at  two  cross-sections 
perpendicular  to  the  C.I.B. ,  and  located  at  stations 
yw  s  74.3  and  yw  =  9#.7>  were  determined  from 


=    Y?  A  q  (1) 


The  calculations  and  results  are  shown  in  Appendix 

A.   A  linear  interpolation,  using  the  planform  dimensions 

from  Fig.  8,   gave  the  shear  flow  for  each  cell,  as  oriented 

perpendicular  to  the  C.I.B.  These  individual  cell  shear 

flows  were  then  resolved  in  stresses  along,  and  perpendicular 

to,  the  streamwise  intermediate  rib  at  station  yw  s  &1.9S. 

From  Mohr's  circle  of  stresses,  shown  in  Fig.  9,  there  is: 

0  s  (90  -  2A) 
fn-  fs  sin  0 

fp=  fs  COS  0 
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FIG.  8 

CENTER  SECTION  OF  F8U-3  PORT  WING 
(Main  Structural  Elements) 
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Fig.  9 
Mohrfs  Circle  of  Stresses  at  Intermediate  Rib 

Applying  these  stresses  to  the  panel  edge  along  the  stream- 
wise  rib,  and  taking  average  panel  thickness  t  and  edge 
distance  wj  the  forces  Fn  and  Fp  due  to  streamwise  shear 
stress  and  stress  normal  to  the  rib,  fs  and  fp,  are: 

Fn  =  f n  w  t  =  fw  sin  G  wt  =  qs  w  sin  G 
Fp  =  f  p  w  t  »  f  s  cos  G  wt  ■  qs  w  cos  e 

It  can  be  seen  that  the  streamwise  component  may  be 
carried  03/"  the  rib  flanges.   However,  the  force  component 
perpendicular  to  the  rib,  Fp  in  Fig.  10,  cannot  be  carried 
by  the  shear  panel  and  must  be  separated  into  a  streamwise 
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yw  =25.09^- 

Fig.  10 

Component  of  Shear  Force 
Perpendicular  to  Intermediate  Rib 


component,  Aqn»  which  is  additive  to  the  force  Fn,  and 
a  component Ap  taken  axially  by  the  adjacent  flanges,  as 
shovm  in  Fig.  11. 


Fig.  11 
Vectorial  Relation  Between  AQn  and  AP 


One-half  of  theAp  component  was  assigned  to  each  adjacent 
flange.  The  various  load  components  and  the  resulting 
applied  loads  on  the  idealized  model  are  given  in  Table  II. 
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3.   Element  Loads 

Loads  were  assigned  to  the  individual  structural 
components  as  shown  in  Figs.  12  and  13.  Only  the  redundants 
and  sufficient  loads  to  define  the  individual  member 
flexibilities  and  general  state  of  stress  were  included. 

The  statically  determinate  system  was  defined  by 
making  cuts  inboard  of  the  intermediate  rib  at  station 
yw  =  #1.9$.   By  "cutting"  the  five  interior  beam  webs 
and  all  flanges  except  the  lower  front  and  two  rear  flanges, 
a  minimum  determinate  and  stable  structure  was  obtained. 
A  cross-section  of  the  determinate  structure  is  shovm 
in  Fig.  14. 


+  y 


4H-+  x 


Lower 
Rear  Flange 


Fig.  1*f 
Cross-section  of  Determinate  System 
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Since  the  intermediate  "tip"  rib  at  station  yw  =  £l.9#  was 
considered  to  be  determinate,  there  were  16  redundant  and 
33  determinate  loads.   The  load  numbering  system  shown  in 
Fie;.  13  was  established  in  order  to  facilitate  organization 
of  the  analysis  method,  and  to  assist  in  conditioning  a 
matrix  for  later  inversion.   Applied  loads  are  therefore 
numbered  P,  through  P33,  and  the  determinate  reactions  are 
numbered  a^  through  aoo.   The  redundant  loads  were  then 
assigned  numbers  aoi  through  a, g. 

Shear  panel  edge  loads  were  proportioned  by  first 
assigning  a  unit  shear  force  to  the  outer  edges  of  all 
cover  panels  and  beam  webs.   Unit  shear  loads  were 
assigned  to  the  forward  edge  of  the  rib  webs.   By 
representing  a  typical  panel  as  shown  in  Fig.  15,  where  a^ 
is  the  unit  load,  there  is: 


ab 

=  ai 

h 

wc 

ac 

=  a. 

W-j 

wc 

ad 

=  ai 

h 

wc 
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Fig.  15 
Proportioning  of  Panel  Edge  Loads 

The  geometry  of  the  panels  and  the  ratios  necessary 
for  determining  the  loads  are  shown  in  Table  III. 

4.   Bending  Torsion  Interaction 

Interaction  between  bending  stresses  and  cover  shears 
was  accounted  for  by  the  method  reported  in  Ref.  4.  Swept 
panels  were  assumed  to  have  parallel  edges,  as  shown  in 
Fig.  16.   A  rectangular  panel  was  then  formed  with  a  constant 
shear  flow  equal  to  the  average  of  the  end  shear  flows,  as 
given  in  Fig.  13. 
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Fig.  16 
Allowance  For  Interaction  of  Bending  Stresses  and  Cover  Shears 

The  triangular  segment  remaining  was  placed  in 
equilibrium  by  reacting  one-half  its  spanwise  component 
at  each  spar  cap  by  the  load  P,  where 

P  ■  q  w  sinAp 

Interior  spar  caps  receive  a  contribution  from  panels  on 
either  side.  These  loads  are  shown  by  dashed  arrows  in 
Fig.  13. 
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5.   Equilibrium  Equations 

Since  the  redundant  loads  were  to  be  obtained  by 
minimizing  the  internal  strain  energy,  33  independent 
equilibrium  equations  were  required  to  establish  the  load 
distribution  in  the  determinate  structure.  Because  there 
were  more  than  33  members  in  the  structure,  some  choice 
was  available  in  writing  the  equilibrium  equations.  It 
was  decided  to  avoid  where  possible  the  writing  of  equations 
for  those  members  or  combinations  of  members  where  it  was 
likely  a  small  load  would  be  applied  or  resisted. 

The  first  six  equations  were  obtained  by  expressing 
the  equilibrium  of  the  entire  model  in  relation  to  the 
axis  system  chosen: 

1.)  X>x  =  0 

2.)  J>y  =  0 

3.)  Df  =  0 

4.  )  ZMxx  =  0 

5.)   EMyy  =  0 

6.  )  5>zz  =  0 

Using  the  right-hand  rule,  each  applied  load  and  its  "external" 
reaction  may  contribute  incremental  moments  according  to 

AM^  =  (Fz)±  Y±  -  (Fy^Zi 
AMyy  =  (Fz)i  X±  -  (?x)iZ± 
AMZZ  =   (Fy)i  Xx  -  (Fx)i  Yi 
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Cover  shears  at  either  end  were  assumed  to  act  at  the 
forward  or  top  edge  of  the  panel  or  beam  web. 

The  delta  moments  are  then  summed  and  set  eoual  to 
zero. 

The  next  eleven  equations  were  obtained  from  axial 
equilibrium  of  the  11  redundant  beam  flanges,  elements  16, 
IS,  20,  22,  24,  29,  31,  33,  35,  37,  and  39. 

Ten  more  eauations  were  obtained  from  the  rib  at 
station  yw  -   Sl.9#  by  considering  the  equilibrium  of  inter- 
beam  segments  between  the  R.B.  and  R.I.B.,  R.I.B.-A.I.B. , 
A.I.B.-C.I.B.,  C.I.B.-FWD.I.B.,  and  FV/D.I.B.-F.I.B.   For  a 
rib  segment,  looking  inboard  from  the  tip,  the  moments  about 
the  y-axis  (lower  forward  corner  of  the  segment)  and  forces 
in  the  x-direction  were  set  equal  to  zero,  giving  two 
equations  per  segment. 

Following  Ref.  5,  "the  vertical  shear  loads  in  the  beam 
webs,  applied  at  station  yw  =  Si. 9^  by  loads  Pj_  through  Py , 
were  assumed  to  act  through  seven  imaginary  posts  of  unit 
cross-section,  as  shown  in  Fig.  12.  The  final  six  equations 
were  obtained  by  setting  all  but  the  F.I.2.  post  in  equilibrium. 

This  gave  a  total  of  33  equilibrium  equations.  Each 
equation  was  arranged  internally  so  that  only  the  non- 
redundant  (a]_  through  aoo)  loads  appeared  on  the  left  side 
of  the  equation.  The  applied  loads  (P^  through  Poo)  and 
redundants  (a^^  -  £49 )  appeared  on  the  right  side. 
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6.   Matrix  Formulation 

Following  the  matrix  formulation  given  by  References 
12  and  13,  all  at  the  equilibrium  equations  can  be  expressed 
by 

W33x33H3x1     =     H&»     ^  H*1  (1) 

Solving  equation  (1)  for  the  non-redundant  loads  anr 
requires  the  inversion  of  the  matrix  of  coefficients,  [A  . 
Therefore  the  row  order  of  the  33  equilibrium  equations  was 
arranged  to  favorably  condition  the  A  matrix  by  either 
placing  the  large  terms  on  the  main  diagonal  or  symmetrically 
placing  groups  of  terms  about  the  main  diagonal.  This  in 
turn  fixed  the  arrangement  of  the  [b|  matrix.  The  |A|  and 
the  IbI  matrices  are  shown  in  Tables  IV  and  V. 

Since  there  are  16  redundants,  16  additional  equations 
must  be  obtained  by  minimizing  the  strain  energy  of  the 
system. 

From  Refs.  4  and  5,  the  internal  strain  energy  of  the 
structure  may  be  written  as 


U  =   1    [^  ...  aj   [p]  (a,  ...  an} 


(2) 


2  E 
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where  elements  of  the  symmetric  IfI  matrix  are  the  sums  of 
the  various  member  flexibilities.   These  flexibilities  were 
obtained  as  shown  in  Appendix  B.   The  IfI  matrix  was 
designated  [UCj  for  the  digital  computer  program,  and 
together  with  the  |Aj  and  [b|  matrices,  formed  the  .basic 
structural  input  to  the  digital  computer  program  shown  in 
Appendix  C. 

Referring  again  to  equations  (1)  and  (2)  it  is  necessary 
to  express  the  strain  energy  in  terms  of  the  external  forces 
and  redundant  loads. 

Equation  (1)  was  re-arranged  by  taking  the  inverse  of 
matrix  IaI  ,  and  post  multiplying  the  inverse  1  A""-*-!  by 
matrix  |B  ,  giving  the  equation 


<■  nrJ  33x1    L  J33x33  L  J  33*b*  ar>     33x^9  \  aj 


fP 

The  column  matrix  {anrj  and  the  [CJj  matrix  were 
increased  to  49x1  and  49x49  dimensions  respectively,  by 
simply  adding  the  trivial  equations 


a3*f  =  a3>+ 
a35  =  a35 


(4) 


aif9  =  a^9 
The  augmented  matrix  |CJJ  was  then  redesignated  the 
IgI  matrix. 
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The  transpose  of  matrix      IgI   x^as   then  taken,    allowing 
the   strain  energy  matrix    I H I    to  be  formed : 

[h]  =    M  [P]  [g]  (5) 

49x^9  *+9x>+9     JU9x^f9       >+9x^9 

By  substituting  equation   (5)   into  equation   (2),   the  internal 
strain  energy  may  now  be  written  as 

U  =  ^If,    ...P33  a3lt..a39Ji^    [h]{pv.P33  a^.  .a^}  (6) 

The  strain  energy  is  then  minimized  by  differentiating 
equation  (6)  and  setting  the  result  equal  to  zero,  according 
to 

ir  =  0  <7) 

Since  H  is  a  symmetrical  matrix,  it  may  be  shown  that 
the  operations  of  equation  (7),  as  performed  upon  the  expand- 
ed form  of  equation  (6),  may  be  expressed  simply  by 


Wi«  =  -[H22"1UieN16x33W 


(B) 


6x33  l  J33x1 


The  sub-matrices  [H22J  anc*  F21J  come  ?rom   the 
partitioning  of  |H|  according  to: 
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[H] 


^9x^9 


H 


11 


pxp 


H 


21 


rxp 


H 


12 


"pxr 


H 


22 


rxr 


p  =  applied  load=  33 
r  =  redundant=l6   (9) 


h9xk9 


By  defining  the  product  of   pop"  J  ^   lH21 J  as : 


[en]    =  -[h22"11    [h21] 

L    J  16x33  i22     J  16x16  L  ^1J1 


6x33 


(10) 


tie  answers  to  the  redundant  loads,  as  given  in  the  computer 
program  of  Appendix  C ,  are : 


K)16x1  =  [»L_,  {'}. 


16x33      L     33x1 


The  total  loading  system  is  then  determined  from 


(db) 


W       =  fsl         (p)  (9) 

l  J  k 9X1     L  J  ^ 9x33    33x1 
where  the  matrix  [si  is  constructed  from  [EN]  and  [g]  by 

32? 


[■]  ■  H. 


^3x33 


[EN] 


16x33 


(10) 


^■9x^9 

^9x33 
Finally  the  internal  strain  energy  is  written  in  terms 

of  the  applied  loads  and  the  flexibility  influence 

coefficients  [cl  as: 
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U  =  i  [P,  ....P33J[c]33x33[Pl  ....  P33J     (11) 


By  substituting  equation  (&b)  into  equation  (6),  the 
matrix  of  flexibility  influence  coefficients  is  given  by 

M  =   i([Hii]  -  MM)  (12) 


7.   Computer  Programming 

The  solution  of  equations  (1)  through  (12)  was  adapted 
by  a  FORTRAN  program  to  the  CDC  1604  computer.  The  variable 
names  and  abbreviations,  flow  charts,  programs  and  sub- 
routines are  given  in  Appendix  C. 

Considerable  effort  was  required  to  establish  the 
desired  program  within  the  basic  computer  storage  capacity. 
Any  future  additions  to  the  main  program  would  require  either 
use  of  machine  language  or  an  input-output  program  on 
peripheral  equipment. 

Print-outs  throughout  the  program  were  made  in  order 
to  establish  a  means  for  continuing  accuracy  checks. 
Although  a  complete  hand  solution  of  one  set  of  dal  a  is 
possible  in  many  computer  programs  where  merely  repetitious 
iterations  are  done,  here  only  equilibrium  of  elements  and 
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standard  machine  inversion  checks  could  be  made. 

The  first  machine  check  was  performed  by  writing  a 
short  auxiliary  program  named  CHECK.   This  confirmed  the 
accuracy  of  the  inversion  of  the  matrix  A  as  accomplished 
by  the  single-precision  Gauss  3  subroutine.  The  matrix 
[CHECKJ  was  formed  from 


[check]  =  [a]  [a-1] 

L  J  33x33       LJ  33x33 L      J33x33 


•  •  • 


1 


1 


(13) 


33x33 
and  is  shown  in  Appendix  C. 

The  second  check  was  made  by  comparing  the  member 
loads  of  the  statically  determinate  structure  with  values 
from  the  ICJJ  matrix  shown  in  Table  VI.  For  example 
consider  the  statically  determinate  structure  shown  in 
Fig.  17. 


Y»  =  2S.  095 


Ky=8/.98 


Fig.  17 
Loaded  Determinate  Structure  For  [CJJ  Check 
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From  statics   a  value  of  a22  maY  be  obtained  for  a  unit 
load   Py.      This  should   compare  with  the  value     c^p2  n     ^n 
the   expression 

|anrl  =    fcj]  I—  I  (14) 

I  nri33x1  L     J  33x*+9 1  aj^xl 

where  all  terms   in  the  column  load  matrix  are  zero  except 

Setting  the  y-component  of  a22  equal  to    (a22)v,   by 
statics : 

(a22>v     =  p?  x  (81.98  -  25.095)   =     7-90069^  P7 
*  '  7T20  ' 

a22         =     7*90069^     Py       =     8.86583  P? 
.8911^ 

From  equation   (14),   and  the  matrix    I  CJJ 

a22  =     P1    (8.8660330)   +   ...   +     Py   (8.8660330) 

which  checks  well  to  the  fourth  decimal. 

By  similar  reasoning,  the  three  non-redundant  flanges, 
elements  14,  26  and  27,  should  react  directly  the  loads 

P27>  P26»  p33# 

Applied  loads  Reaction  by  [CJJ  matrix 

P2y  ■  1.0  a21  =  1.000055$ 

P26  =  1.0  a22  =  1.0000268 

P33  =  1.0  a23  =  1.000023$ 
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Again  considering  a  unit  load  at  P«,  the  total  shear 
load  should  be  carried  by  the  rear  beam  web,  ag.  From  the 
IcjJ  matrix,  a  unit  load  at  Pn   gives 

ag  =  1.0103419 

Consequently,  the  cover  panel  shear  loads  should  be 
zero.   From  the  |CJJ  matrix, 

a^  =  .0000646  al5  =  .0000437 

a10  =  .0000591  a16  =  .0000471 

an  =  .0000550  a1?  =  .0000511 

a12  =  .0000510  alg  =  .0000550 

a13  =  .0000471  a19  =  .0000591 

a14  =  'OOOOI+38  ^20  =  •000°645 

The  next  check  was  to  determine  the  state  of 
equilibrium  of  each  flange  and  post.  For  example,  the 
front  post  of  the  intermediate  rib,  station  y  =  $1.9#, 
Fig.  12,  gives 

1.0000000  -  .5981273  -  .^018727  =  .0000000 


for  a  unit  load  at  P, .  For  the  seventh  post,  the  result 
from  a  unit  load  at  P7  is  also  zero.  A  maximum  deviation 
from  true  equilibrium  of  7.9  percent  existed  at  the  second 
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post  from  the  front.   It  should  be  noted  that  no 
equilibrium  equation  was  written  for  this  post. 

Referring  to  Fig,  13,  a  unit  redundant  load  was  applied 
in  turn  to  each  redundant  flange.  The  greatest  deviation 
from  equilibrium  was  in  element  39,  where 

Error  r  ,  =  -  O.OOOOOV}  lb 

[cj] 

Table  VII  contains  the  results  of  the  CJ  matrix 
equilibrium  checks,  which  completes  the  check  of  the  basic 
equations  and  computer  operations  on  them,  as  well  as  the 
validity  of  the  matrix  inversion  sub-routine. 

A  print-out  statement  was  used  to  determine  that  proper 
partitioning  of  the  H|  matrix  occured.  Since  computer 
storage  space  was  not  available,  hand  calculations  were 
performed  to  confirm  the  proper  inversion  of  H^?]  • 

An  interim  calculation  according  to  equation  (£b) 
gave  the  redundant  load  matrix  {ar}  .  One  of  these  values 
was  computed  by  hand  as  a  random  check  and  found  to  be 
correct.  All  16  values  were  then  checked  with  the  final 
16  loads  of  the  {AT}  matrix  obtained  from  equation  (9). 

This  completed  the  computer  program  checks.  The 
program  was  then  run  twice  and  answers  compared  in  order  to 
detect  any  computer  core  errors.  There  were  none. 

The  flexibility  influence  coefficient,  stress  influence 
coefficient,  and  load  matrices,  [cl,  [s]  ,  and  {AT}  ,  are 
listed  in  Tables  VIII,  IX,  and  X. 
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8,        Theoretical  Results 

Compared  with  ordinary  beam  theory  solutions,  the 
accurate  determination  of  stresses  in  an  idealized  structure 
by  a  matrix  force  method  requires  considerable  engineering 
judgment.   It  has  been  pointed  out  by  V/arren  (Ref.  9)  that 
errors  of  one-hundred  percent  may  occur  in  interpreting  the 
results  of  such  a  solution.  Therefore  the  original 
formulation  assumptions  were  used  in  determining  the  results 
of  this  analysis. 

The  investigation  was  narrowed  to  determining  the 
average  panel  stress  at  station  yw  =  34.5  from  the  loads  of 
the  [kr.j   matrix.  For  each  panel,  the  maximum  normal  stress, 
the  angle  of  the  principal  axis,  and  the  maximum  shearing 
stress  were  determined.  From  the  components  shown  in  Fig, 
IS,  and  noting  that  the  direction  of  shear  flow  reflects 
the  minus  sign  indicated  by  load  an  in  the  {AT}  matrix, 


a^  +  .5685  an 


a^  +  .5685  aft 
+  .5358  a10 


Fig.   1$ 
Theoretical  Root   Panel  Stresses 
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the  axial  load  components  of  the  skin  v;ere  determined. 
Using  the  solution  formulation  data  of  Table  I,  the  bar 
loads  (excepting  &iq)   were  re-apportioned  by  area  ratios  to 
the  panels  from  which  the  skin  came.  For  panel  3#,  shown 
in  Fig.  IS, 


*97*+  a^  +  1.03  ai^  +  1.137  a9   =  axial  load 
.6^-93(5^61)  +  .3388(8682)  +  1.137(-3785)  =  2l8*f 


or,  axial  panel  load  =  21#4  lb  compression.  This  compressive 
force,  when  divided  by  the  panel  and  cross-section  (taken 
perpendicular  to  the  panel  sweep  angle  ,-^p ) ,  gives 


A, 


218^  =  218*+ 

tw      .lWx  20.^9  x  .73865 


p  =    978  psi  (compression) 


The  shearing  stress  was  found  by  taking  the  average 
shearing  force,  given  for  either  end  of  the  panel  03'  the 
{AT}  matrix,  and  dividing  it  by  the  same  tw  used  in  deter- 
mining the  axial  stress.  The  shearing  stress  and  axial 
stress  thus  found  were  resolved  by  Mohr!s  circle  to  a 
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maximum  principal  stress t(j  max,  a  maximum  shearing  stress, 
Tnax>  an(*  t^ie  an^^e  °^  tne  principal  axis,  (h  p.  These 
stresses  and  angles,  as  listed  in  Table  XI,  refer  to  the 
mid-plane  of  the  skin,  since  average  values  are  used  for 
their  determination.  These  values  are  shown  graphically  in 
Fig.  25. 

The  vertical  deflections  of  the  intermediate  rib 
(station  y  =  81 .98)  end  points  were  calculated  from  the 
deflection  influence  coefficient  matrix  [Cj  by  applying 
the  theoretical  vertical  loads  P1  through  ?„   •  At  the 
forward  end  of  the  rib(load  point  for  P1  )  the  deflection 
was  +  .0601  in.  At  load  point  seven,  the  deflection  was 
-. 06*+ 7  in.  These  values  are  plotted  in  Fig,  26. 
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Ill  EXPERIMENT 
1.   Introduction 

In  order  to  determine  the  actual  stresses  and  compare 
them  with  those  predicted  by  theory,  an  experimental 
investigation  has  been  carried  out  on  the  full  scale 
F#U-3  wing  acquired  from  Chance  Vought  Aircraft, 
Incorporated.  The  experimental  set-up  was  also  designed 
to  provide  the  Aeronautical  Structures  Laboratory  of  the 
U.  S.  Naval  Postgraduate  School  with  a  permanent  specimen 
of  modern  wing  construction  mounted  such  that  a  variety  of 
test-to-theory  correlation  experiments  could  be  academically 
demonstrated. 

The  F#U-3  wing  was  structurally  complete  in  every 
detail  and  was  representative  of  current  doubly  tapered 
multi-cell  swept  wing  configurations.  It  was  divided  into 
three  major  sections,  the  center  section  which  extended 
to  the  wing  fold  rib  at  wing  station  y  =  152.25  and  an 
outer  panel  on  either  side.  For  this  experiment  the  two 
outer  panels  were  removed  along  with  the  leading  edges  and 
those  portions  of  the  trailing  edges  outside  of  the  main 
structural  box  that  were  not  integrally  connected  with  its 
single  piece  skin.   Hereafter  the  term  wing  will  refer  to 
the  main  structural  elements  of  the  center  section  which 
consisted  of  seven  spars,  hereinafter  called  beams,  and 
two  streamwise  ribs.  The  section  was  a  typical  thick 
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milled  skin  construction  and  the  skins  were  tapered  along 
the  span  as  well  as  in  the  chordwise  direction. 

Structural  behavior  was  investigated  with  the  wing 
subjected  to  two  pure  torque  loads  applied  perpendicular 
to  the  elastic  axis  at  the  wing  fold  rib.   The  elastic 
axis  is  defined  by  the  curve  connecting  the  shear  centers 
of  all  sections  of  the  wing.   Its  position  was  taken  from 
manufacturer^  data  as  shown  in  Fig.  £.  The  torque  loads 
were  applied  to  each  wing  fold  rib  by  hydraulic  cylinders 
actuated  from  a  common  pressure  manifold.   The  load 
magnitudes  were  measured  by  dynamometers  connected  in 
series  with  each  cylinder  and  the  torque  was  then  easily 
computed  knowing  the  fixed  lever  arm.   Nominal  torque 
values  of  294000  in.  lb  and  336000  in.  lb  were  used.  These 
were  well  within  the  elastic  range  but  still  of  sufficient 
magnitude  to  give  adequate  strain  levels  throughout  the 
structure  for  repeatably  accurate  measurements. 

Some  electrical  strain  gages  were  already  installed 
by  Chance  Vought  Aircraft,  Incorporated,  at  various 
locations,  but  for  the  purposes  of  this  and  future  invest- 
igations additional  SR-4  strain  gage  rosettes  of  the  AR-7-2 
and  A-7  type  were  installed  in  particular  areas.   Those  of 
immediate  interest  were  located  in  the  skins,  beam  webs  and 
beam  caps  around  a  section  perpendicular  to  the  center 
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intermediate  beam  at  wing  station  yw  =  9S.7.   From  the 
wing  fold  to  this  point  the  elastic  axis  remains  nearly- 
parallel  to  the  center  intermediate  beam  which  was  used 
as  the  reference  axis  of  the  wing.   Other  rosettes  were 
located  around  a  streamwise  section  near  the  root  at  wing 
station  yw  =  34.5.   Here  rosettes  were  placed  on  both  sides 
of  the  very  thick  skin  to  check  for  differential  bending 
effects.   It  was  expected  that  the  stresses  that  are 
unpredictable  by  simple  beam  theory  in  this  area  would  show 
closer  agreement  with  the  predictions  of  the  matrix  force 
method  of  analysis.   Strain  gage  data  was  programmed  into 
a  digital  computer  to  obtain  principal  stresses  and  directions 
at  all  rosette  locations. 

Deflections  were  measured  by  means  of  scales  hung  from 
the  structure  along  the  front  and  rear  beams  and  read  by  a 
transit.   Jig  support  deflections  were  also  measured  in  this 
manner. 
2.   Equipment 

The  variable  incidence  swept  wing  was  of  a  conventional 
multi-cell  construction  with  a  quarter-chord  sweep  back  of 
1+2   degrees.   The  main  structural  elements  are  shown  in 
Fig.  $,     Fig.  19  is  a  photograph  of  the  wing  with  the  upper 
skin  removed.   Attachment  was  made  to  the  fuselage  at 
approximately  wing  station  yw  =  25  by  two  pivot  lugs  located 
approximately  5-3  inches  aft  of  the  main  box  rear  beam, 
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two  bumpers  located  at  the  main  box  front  beam  and  two 
incidence  actuators  located  approximately  17.  £  inches 
forward  of  the  main  box  front  beam.   The  wing  fold  at  wing 
station  yw  *  152.25  was  accomplished  by  fittings  with 
multiple  upper  and  lower  lugs.  The  upper  and  lower  surfaces 
were  single  piece  thick  skins  tapered  along  the  span  and 
chord.  The  front  beam  had  numerous  supporting  lugs 
integrally  machined  into  the  beam  giving  it  a  very 
discontinuous  variation  of  thickness.  The  structure 
primarily  consisted  of  7079-T6  thick  skins  and  forgings 
with  7075-T6  sheet  metal  beams. 

The  wing  was  mounted  inverted  on  a  rigid  support  jig 
and  fastened  to  it  at  four  points,  the  two  fuselage  pivot 
lugs  and  the  two  points  on  the  main  box  front  beam  directly 
opposite  the  bumper  points.   Plywood  pads  were  used  under 
the  jig  to  distribute  the  loads  to  the  laboratory  floor. 
A  photograph  of  the  mounted  wing  is  shown  in  Fig.  20. 

The  torque  loads  were  applied  through  a  fitting 
connected  to  the  wing  fold  lugs  as  shown  in  Fig.  6.   The 
fitting  assembly  was  designed  to  provide  a  loading  plane 
perpendicular  to  the  plane  of  the  wing  and  the  elastic  axis. 
It  consisted  primarily  of  mild  steel  parts  designed  to 
develop  the  ultimate  strength  of  the  wing  fold  lugs.   A 
limiting  torque  of  420000  in.  lb  was  selected  which  was  well 
below  the  design  ultimate  at  this  station  but  still  high 
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enough  to  ensure  adequate  strain  levels  in  the  wing.   A  42 
inch  lever  arm  was  used  to  permit  attainment  of  the  420000 
in.  lb  of  torque  with  loads  of  only  10000  lb.   This  enabled 
the  use  of  available  10000  lb  dynamometers  which  were 
graduated  in  100  lb  increments.   Instruments  of  larger 
capacity  had  250  lb  increments. 

Only  the  port  wing  was  instrumented  but  the  torque  was 
applied  to  each  wing  fold  station  for  balance  and  to 
minimize  warping  of  the  structure  and  support  jig.  Each 
of  the  four  loads  were  applied  through  a  series  linkage  of 
attachment  fittings,  clevises,  hydraulic  cylinder  and 
dynamometer  as  shown  in  Fig.  21.   The  four  hydraulic 
cylinders  were  high  pressure  wing  fold  actuators  and  were 
simultaneously  subjected  to  pressure  from  a  common  manifold 
to  insure  load  equalization.  Except  for  the  dynamometers 
each  linkage  was  designed  for  30000  lb  and  was  pre-tested 
to  15000  lb  before  installation. 

Load  magnitude  was  determined  and  monitored  by  two 
separate  means.   On  the  instrumented  wing  two  10000  lb 
Dillon  Dynamometers  were  used  to  monitor  equality  of 
individual  loads.  They  were  graduated  in  100  lb  increments 
but  could  be  accurately  read  to  10  lb.   The  final  load 
magnitudes  were  determined  by  hydraulic  manifold  pressure 
related  to  pressure-load  calibration  curves  for  the 
cylinders  as  found  in  Appendix  D.   The  hydraulic  gage  was 
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mounted  at  the  data  taking  station  and  provided  an 
excellent  means  of  setting  the  load  and  monitoring  the 
consistancy  of  its  magnitude  throughout  the  run.  It  was 
a  1500  psi  gage  graduated  in  10  psi  increments  and  could 
be  accurately  read  to  5  psi. 

The  dynamometers  and  the  hydraulic  system  were 
calibrated  on  a  300,000  lb  Riehle  Tensile  Testing  Machine 
(Appendix  D).  The  Reihle  machine  had  exhibited  a  maximum 
calibration  error  of  only  0.35  percent  over  the  entire 
300,000  lb  range,  but  showed  no  error  in  the  range  of 
7000  lb  to  9000  lb,  which  was  of  concern  in  this  case. 
Therefore  the  calibration  curves  represent  variation  from 
absolute  values. 

The  up-load  linkages  were  anchored  to  a  three  column 
supporting  structure  shown  in  Fig.  21.   The  base  of  this 
structure  was  arranged  to  accommodate  attachment  of  the 
down-load  linkage  thereby  minimizing  strength  requirements 
of  the  base  elements.   Considering  future  uses  of  this 
structure  the  base  was  located  approximately  in  the  plane  of 
the  load  fitting  such  that  single  point  loads  could  be 
applied  anywhere  along  its  reach  and  be  readily  anchored 
with  a  minimum  of  effort. 

Hydraulic  pressure  was  supplied  by  an  electrically 
driven  Vickers  V-line  Piston  Type  Pump  shown  in  Fig.  22. 
It  was  rated  at  5000  psi,  1800  rpm  and  1.72  cu  in.  per 
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revolution.   By  means  of  a  pressure  compensator  in 
conjunction  with  a  volume  limit  hand  control  and  relief 
valve,  pressures  could  be  held  within  5  psi  of  any  desired 
setting  for  long  period  of  time  and  with  very  few  adjust- 
ments.  This  feature  enabled  accurate  repetition  of  the 
same  loads  which  was  essential  in  the  data  taking  process. 

All  strains  were  measured  with  SR-4  strain  gage 
rosettes  manufactured  by  Baldwin-Lima-Hamilton  Corporation. 
In  Appendix  E  is  a  list  of  the  4#5  gaSes>  their  gage  factors, 
resistances  and  coordinate  locations  in  the  wing.  The  types 
of  gages  installed  by  Chance  Vought  Aircraft  were  deduced 
from  resistance  tests  and  inspection.   They  were  all  of 
the  AX-5  or  AR-1  type.   The  gages  installed  for  this 
experiment  and  for  future  multi-purpose  investigations  were 
either  of  the  AR-7-2  or  A-7  type.  The  gages  of  particular 
interest  to  this  investigation  are  depicted  in  Fig.  23. 
Selected  root  skin  rosettes  are  backed-up  on  both  sides  of 
the  skin.   Although  additional  back-up  rosettes  were 
desired  throughout  the  inboard  panels,  they  could  not  be 
installed  because  of  time,  laboratory  priority  and  monetary 
considerations.  Note  that  all  of  the  gages  have  a  gage 
factor  very  close  to  2.0  and  resistances  of  120  ohms.   This 
enabled  the  use  of  a  single  temperature  compensating  gage 
and  factilitated  reading  all  gages  without  adjusting  any 
of  the  other  instruments. 
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Selectivity  in  reading  any  combination  of  gages  was 
made  possible  by  routing  all  gage  leads  to  a  junction 
panel  as  shown  in  Fig.  24.   Any  number  of  gages  could  then 
be  connected  into  either  switching  and  balancing  units  or 
automatic  scanning  devices  as  future  needs  dictate.   All 
electrical  connections  were  soldered  to  minimize  contact 
resistance  except  for  the  banana  plug  connections  between 
the  junction  panel  and  the  switching  and  balancing  unit. 

A  20  channel  Baldwin -Lima- Hamilton  Switching  and 
Balancing  Unit  was  used.   It  was  connected  to  an  external 
Wheat stone  bridge  circuit  powered  by  a  Hewlett  Packard 
Power  Supply  (Model  721A).   A  constant  6  volts  was  used 
throughout  the  tests.  The  bridge  output  was  amplified, 
then  fed  into  a  voltage-to-frequency  converter  and  displayed 
as  strain  in  units  of  micro-inches  per  inch  by  an  electronic 
counter.   Calibration  of  this  equipment  proceeded  each  run 
and  the  method  is  described  in  Appendix  F. 

The  electronic  counter  was  a  Model  521DR  manufactured 
by  Hewlett  Packard  Corporation.   The  voltage-to-frequency 
converter  was  a  Model  DY-2210  manufactured  by  Dymec 
Incorporated  and  the  amplifier  was  a  Kintel  Model  111BF. 
By  using  the  electronic  counter,  readings  could  be  taken 
faster  and  more  accurately  than  with  the  common  strain 
indicators.   Strain  indicators  can  normally  be  read  to 
about  5  micro-inches  per  inch,  whereas  the  electronic  counter 
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indicates  strains  to  one  micro-inch  per  inch. 

Wing  deflections  were  measured  by  means  of  scales 
hung  at  intervals  along  the  front  and  rear  beams  and  read 
by  a  transit  manufactured  by  Keuffel  and  Esser  Company. 
The  scales  were  graduated  in  .02  inch  increments  but  could 
be  read  to  .005  inch.  The  locations  of  the  scales  are 
shown  in  Fig.  £.   Note  that  a  single  scale  could  be  paired 
with  one  of  two  others  such  that  either  a  streamwise  chord 
or  a  section  perpendicular  to  the  center  intermediate  beam 
would  be  defined.  For  instance,  scale  numbers  8  and  7 
define  the  test  section  perpendicular  to  center  intermediate 
beam  at  wing  station  yw  ■  93.7  and  number  £  and  11  terminate 
the  streamwise  rib  at  wing  station  yw  =  £l.9#.  This  was 
believed  to  be  of  use  in  future  investigations. 
3.   Experimental  Procedures 

The  experimental  procedures  consisted  of  two  tests, 
in  each  case  a  different  magnitude  of  torque  being  applied 
at  the  wing  fold  rib.   In  the  first  test  a  nominal  294000 
in.  lb  was  applied  and  336000  in.  lb  in  the  second  test. 
These  values  were  initially  attained  by  dynamometer  indi- 
cations of  7000  lb  and  #000  lb  respectively,  which  will  be 
the  names  of  the  two  tests  hereafter.  The  purpose  of  the 
7000  lb  test  was  to  establish  reliability  of  the  entire 
4#5  strain  gages  and  to  check  for  linearity  of  strain 
readings  in  the  structure.  The  #000  lb  test  was  accomplished 
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to  correlate  the  theoretical  results.  Therefore  in  the 
£000  lb  test  attention  was  focused  on  approximately  150 
gae-es  of  particular  interest. 

The  procedure  in  conducting  the  tests  was  straight 
forward,  but  because  of  the  large  number  of  strain  measure- 
ments the  tests  took  considerable  length  of  time.  Only  20 
gages  could  be  read  during  any  one  run  because  only  one 
switching  and  balancing  unit  was  employed.  Additional  units 
were  tried  but  the  added  wiring  and  connections,  resistance 
peculiarities  inherent  in  each  unit  and  prolonged  time 
required  to  adjust  initial  zeroes,  produced  an  unacceptable 
drift  in  the  zero  rechecks.  On  the  other  hand  load 
repeatability  was  very  accurate  and  therefore  a  single  unit 
was  used  and  the  number  of  runs  increased.  The  sequence  of 
operation  for  each  run  was  as  follows : 

1.  The  electronic  counter  was  calibrated  before  each 
series  of  runs  and  rechecked  after  completion. 

2.  The  switching  and  balancing  unit  leads  were  plugged 
into  20  gage  terminals  at  the  junction  panel. 

3.  All  strain  gages  and  dynamometers  were  zeroed. 

4.  The  hydraulic  loads  were  applied  and  adjusted  by 
reference  to  previously  known  values  of  hydraulic  pressure 
and  strain  gage  readings  from  the  same  check  gage  used  in 
each  test  run. 

5.  The  dynamometers  were  read  and  compared  with  each 
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other  to  insure  that  equal  couple  loads  were  being  applied 
to  each  wing. 

6.  The  strains  were  read  and  recorded  twice. 

7.  The  loads  were  removed. 

£.   All  zero  values  of  strain  were  re-checked  and 
recorded. 

Extreme  care  was  employed  to  achieve  load  and  strain 
repeatability.  The  electronic  counter  was  calibrated 
before  each  run  and  checked  afterwards.  A  maximum  departure 
from  linearity  over  a  range  of  1000  micro-inches  per  inch 
was  only  5  micro-inches  per  inch.  The  average  departure  . 
was  closer  to  2  micro  inches  per  inch.  This  effect  on 
maximum  strain  levels  encountered  (approximately  200  micro- 
inches  per  inch)  was  therefore  quite  negligible. 

In  addition  particular  care  was  exercised  to  minimize 
the  strain  gage  zero  drift.  For  instance,  the  variable 
temperature  effects  of  sunlight  and  outside  electrical 
interference  with  the  sensitive  instrumentation  were 
virtually  eliminated  by  conducting  all  tests  at  night.   At 
least  one  hour  warm-up  time  was  given  all  instruments  to 
insure  steady  state  conditions.   All  initial  zeroes  were 
re-checked  three  times  before  loading.  All  permanent 
electrical  connections  were  soldered  to  minimize  contact 
resistance  and,  as  a  double  precaution,  a  test  was  repeated 
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if  any  electrical  leads  had  been  disturbed  or  if  the 
voltage  supply  varied  any  time  during  the  run.  By  reading 
only  20  gages  per  run  the  time  per  run  was  minimized, 
thereby  reducing  the  critical  time  the  system  would  be 
exposed  to  outside  effects.   Two  strain  readings  at  load 
were  taken  also  as  a  double  check.  With  such  care,  average 
zero  rechecks  within  1.5  micro-inch  per  inch  were 
continually  demonstrated. 

Two  strain  gages  were  selected  to  act  as  cross-checks 
on  the  repeatability  of  the  loads,  results  and  test 
conditions.   Gages  32  and  53  were  monitored  in  each  test 
run.   The  magnitudes  of  strain  on  each  run  were  then 
compared  to  the  original  values.   In  over  40  runs  the 
average  departure  from  mean  values  of  the  check  gages  was 
only  1.39  micro-inch  per  inch.  The  mean  value  for  gage 
53  in  the  #000  lb  test  was  17S.6  giving  an  error  of  only 
0.7#  percent.   Since  the  hydraulic  pressure  reading 
reouired  to  produce  this  strain  level  was  always  the  same 
the  load  repeatability  was  also  within  an  accuracy  of  one 
percent. 

Accuracy  of  measuring  the  absolute  load  magnitude  was 
then  dependent  only  on  the  accuracy  of  the  dynamometers  and 
hydraulic  pressure  gage.   During  the  weeks  of  testing, 
difficulty  was  encountered  in  re-setting  the  dynamometers 
to  the  same  zero.   Also  there  was  noticable  differences 
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in  readings  under  load  conditions  which  verified  suspicions 
of  unreliable  zero  settings.   On  the  other  hand  the 
hydraulic  pressure  gage  readings  checked  consistantly  with 
strain  gage  32  and  53  readings  and  therefore  the  pressure 
gage  was  felt  to  be  more  reliable  in  establishing  magnitudes 
of  applied  loads. 

Deflection  measurements  were  made  for  only  the  #000 
lb  test  condition.   Since  this  loading  condition  and 
consequently  the  deflections  were  repeated  over  20  times 
it  was  deemed  necessary  to  only  sample  the  deflections 
periodically  throughout  the  test  series.   Although  more 
deflection  measurements  were  taken  to  ensure  repeatability 
in  the  series  only  three  complete  sets  of  data  were 
recorded. 
4.   Experimental  Results 

From  the  recorded  strain  data  the  only  calculations 
required  were  those  involved  in  determining  the  principal 
stresses,  maximum  shear  stress  and  principal  axes  at  each 
rosette  location.   These  calculations  for  all  rosettes  in 
the  7000  lb  test  and  those  particular  rosettes  of  interest 
in  the  #000  lb  test  were  performed  on  a  Control  Data 
Corporation  1604  Digital  Computer.   The  results  along  with 
the  FORTRAN  program  and  the  measured  strains  are  found 
in  Appendix  G. 
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The  data  for  the  7000  lb  test  was  pertinent  to  this 
experiment  only  to  show  linearity  of  strain  readings.   A 
sampling  of  net  strain  readings  from  various  locations  for 
both  loading  conditions  is  shown  in  Table  XII.  Satisfactory 
linearity  is  demonstrated  by  noting  that  the  average 
difference  between  measured  strains  for  the  #000  lb  test 
and  those  extrapolated  from  the  7000  lb  test  is  2.3  micro- 
inches  per  inch  regardless  of  the  magnitude  of  the 
measurement.   The  discrepancy  is  obviously  quite  small  but 
it  would  be  still  further  diminished  if  the  accuracy  of 
the  applied  loads  were  considered.  Therefore  it  can  be 
concluded  that  linearity  of  strain  readings  was  indeed 
achieved. 

Another  result  of  the  7000  lb  test  was  the  discovery 
of  9  faulty  rosettes.   However  none  of  the  faulty  rosettes 
hampered  the  results  of  this  experiment.  The  faulty 
gages  are  indicated  in  Appendix  E. 

Calibration  of  the  hydraulic  system  after  the  completion 
of  the  tests  revealed  that  instead  of  #000  lb,  actually 
#400  lb  loads  were  being  delivered  to  the  linkages. 
Considering  a  gage  reading  accuracy  of  5  psi  and  the  slope 
of  the  calibration  curve  as  10. #4  lb  load  per  psi  hydraulic 
pressure,  the  accuracy  of  the  applied  loads  would  be  plus 
or  minus  50  lb  or  0.6  percent.  The  difference  in  desired 
and  achieved  magnitudes  is  not  of  real  importance  however 
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since  linearity  of  strains  enables  interpolation  to  the 
desired  values. 

Since  a  degree  of  uncertainty  always  exists  concerning 
the  extent  to  which  a  pure  torque  loading  is  realized,  an 
equilibrium  check  was  made  at  the  instrumented  section 
perpendicular  to  the  CXB,  at  yw  ■  96.7.   At  this  section 
all  rosettes  were  oriented  such  that  the  diagonal  gages 
were  at  a  45  degree  angle  to  the  section,  which  facilitated 
shear  flow  calculations.   Therefore  theoretically 
calculated  shear  flows  could  be  readily  checked  with 
experimental  results  since  the  theoretical  calculations  were 
all  ready  completed  in  the  determination  of  the  loading 
for  the  idealized  structure  as  shown  in  Appendix  A. 

The  experimental  shear  flows  were  calculated  by 
considering  only  the  strains  in  the  diagonal  gages.  This 
essentially  eliminated  the  effects  of  any  normal  stresses 
caused  by  unknown  bending  loads  and  the  results  would 
better  represent  the  shear  flow  induced  by  only  the 
effectively  applied  torque.   From  the  experimental  data 
for  the  #000  lb  load,  shear  flows  were  calculated  using 
the  relation, 


q=  ^-L£2  (15) 

1   +  JJ 


I  . 
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These  results  are  listed  in  Table  XIII  and  are  labled  on 
the  schematic  diagram  in  Fig.  25.   The  values  of  shear 
flow  in  the  top  and  bottom  skin  of  the  same  cell  compared 
closely.  The  average  value  of  each  cell  was  used  as  also 
shown  in  Table  XIII  and  Fig.  25  and  compared  to  the 
theoretical  values.  The  theoretical  shear  flows  for 
T  ■  336000  in.  lb  were  3. #5  percent  higher  than  the  experi- 
mental values.  The  total  torque  resisted  was  then  computed 
as  the  sum  of  the  (2Aq)  for  each  cell.  This  value  was 
323,7^5.  in.  lb  which  is  S.2  percent  lower  than  the 
352,600  in.  lb  possible  with  an  #400  lb  force  acting  on  a 
42  inch  lever  arm,  and  3.6  percent  lower  than  that  used  in 
the  analysis. 

Neglecting  the  accuracy  of  the  experimental  loads  and 
shear  flow  determination  this  means  that  720  lb  of  the 
applied  load  of  #400  lb  was  not  effective  in  producing  pure 
torque.   In  this  regard,  during  the  tests  there  had  been 
visible  evidence  of  streamwise  twisting  of  the  wing  fold 
rib  and  a  measured  spanwise  tilt  of  the  entire  wing.  There 
was  also  an  element  of  uncertainty  regarding  the  exact 
angle  of  the  elastic  axis  since  it  was  scaled  from  a 
reproduced  drawing.   Coupled  with  understandable  inaccuracies 
of  maintaining  a  constant  loading  plane,  perfectly  parallel 
lines  of  force  and  a  perfectly  rigid  loading  jig,  the  amount 
of  the  bending  load  component  is  within  reason.   Although 
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no  calculations  were  made  to  determine  the  magnitude  of 
the  bending  component  its  existance  can  be  qualitatively 
confirmed  by  observing  the  magnitudes  of  the  strain 
readings  for  the  perpendicular  gages  at  the  section  being 
analyzed,  namely  rosette  numbers  131,  46,  49,  52,  55,  5#, 
etc. 

Using  a  torque  value  of  323  ^7&5.  in.  lb  the  theoretical 
shear  flows  were  obtained  by  interpolating  between  values 
computed  for  the  294000  and  336000  in.  lb  condition. 
They  are  also  shown  in  Table  XIII  and  Fig.  25.   As  would 
be  expected  the  comparison  with  experimental  values  is 
very  good.   The  average  error  in  the  theoretical  shear 
flows  is  only  1.95  percent.   Considering  the  loss  of 
effective  applied  torque  the  torsional  equilibrium  check 
was  believed  to  be  quite  close. 

The  computer  program  described  in  Appendix  G  solves 
the  well  known  rectangular  rosette  equations  for  tf^x, 
0^±n,    ^max  anc*  typ*      Tlie  inPut  to  the  program  was  the 
value  of  strain  for  each  gage  in  the  rosette.   The  lowest 
numbered  gage  was  consistantly  called  6-^,  the  diagonal 
and  perpendicular  gages  were  the  next  higher  numbered  gages 
respectively.   The  principal  axis  was  then  computed  with 
respect  to  the  €^  axis.   To  compare  these  results  with 
theoretical  values  the  orientation  of  the  principal  axes 
were  adjusted  to  the  yw  axis  of  the  wing  by  simple  arithmetic. 
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The  root  area  rosettes  used  to  compare  with  results  of 
the  theoretical  root  analysis  were  numbers  13,  19,  2#,  367, 
279  and  161  on  the  inside  surfaces  and  back-up  rosettes 
numbers  4$1,  47$,  472,  466,  457,  and  4$4  on  the  exterior 
surfaces.   Differential  bending  was  readily  apparent  in  all 
of  the  root  panels  as  seen  in  Table  XIV  by  noting  the 
differences  of  strain  in  matching  gages  on  opposing  surfaces. 
The  effect  is  less  in  panel  number  21  which  lies  near  the 
elastic  axis.   The  average  of  the  two  gage  readings  was 
taken  as  representative  of  the  strain  at  mid-thickness 
of  the  skin.   The  values  of  C^x*  "^max  anc*  iP   were  then 
hand  calculated  using  the  equations  in  Appendix  G  and  are 
listed  in  Tables  XI  and  XIV  and  plotted  in  Fig.  26. 

The  deflection  measurements  were  of  concern  primarily  in 
establishing  test  conditions.   In  this  regard  the  average 
values  of  the  total  deflections  are  plotted  in  Fig.  26. 
It  is  clearly  shown  that  the  support  jig  rocks  forward  and 
the  entire  structure  tilts  right-wing-down.   It  was  also 
noted  during  the  test  that  the  amount  increased  slightly 
with  each  cycle  which  indicated  plastic  yielding  of  the 
plywood  pads  under  the  support  jig.   Obviously  any  tilt 
would  adversly  affect  the  lines  of  action  of  the  applied 
loads  since  the  load  linkages  were  positioned  to  accommodate 
only  one  loading  plane.   This  would  account  for  some  of  the 
discrepancy  in  achieving  pure  torque  loads.   It  would  be 
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wise  to  consider  adjusting  the  load  fittings  and  linkages 
such  that  a  more  symmetric  loading  could  be  realized  in 
future  tests.   Furthermore  it  would  be  advantageous  to 
replace  the  plywood  pads  with  a  more  suitable  material. 
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IV  COMPARISON  OF  RESULTS 
Extreme  care  was  taken  to  ascertain  strain  and  load 
(iiioUdes  and  insure  repeatability  of  loads  to  minimize 
any  experimental  error  in  the  determination  of  the  compar- 
itively  low  stress  levels  found  in  the  critical  root  area. 
All  of  these  errors  were  kept  well  below  one  percent  and 
logically  do  not  appreciably  influence  any  comparison  with 
theory. 

The  only  experimental  error  that  would  noticeably 
affect  the  final  results  would  be  the  understandable 
discrepancy  in  developing  pure  torque  loading  on  which 
the  theoretical  analysis  was  based.  The  equilibrium  check 
at  the  section  perpendicular  to  CJ£«  at  yw  "  9&.7  indicated 
that  the  torque  developed  was  3.6  percent  lower  than  that 
used  in  the  analysis  and  nine  percent  lower  than  actually 
applied  at  the  tip.   This  meant  that  about  $.6  percent  of 
the  applied  load  introduced  unwanted  bending  effects  that 
were  not  accounted  for  in  the  analysis.   The  cause  of  the 
torque  discrepancy  is  correctly  attributed  mainly  to  the 
method  of  loading  and  tilting  of  the  support  jig.   Other 
than  reducing  the  accuracy  of  the  experimental  results 
to  a  comparatively  minor  degree  the  effect  of  the  small 
bending  loads  would  not  measureably  influence  the 
comparison  with  theory  in  such  a  complex  structure.   These 

factors  become  even  less  important  in  establishing  the 
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general  validity  of  this  analytical  method  to  an  actual 
complicated  wing  when  one  is  reminded  that  the  literature 
mentions  100  percent  disagreements  in  comparing  results  of 
other  methods  applied  to  simplified  thin-skinned  laboratory 
specimen. 

Experimental  and  theoretical  results  for  the  root  area 
are  listed  in  Table  XI.   Three  observations  may  be  made 
immediately.   The  sense  of  maximum  normal  stress,  the 
orientation  of  the  principal  axis,  and  magnitude  of  stresses 
at  the  critical  rear  panel  area  agree  rather  well.  Maximum 
shear  stresses  differ  by  only  25  percent,  and  although  one 
of  the  maximum  normal  stresses  differs  by  40  percent  the 
other  is  only  seven  percent.   On  the  other  hand,  the 
theoretical  stresses  obtained  at  the  leading  edge  are  much 
too  high.   This  is  shown  graphically  in  Fig.  27.   Since  the 
equilibrium  checks  made  during  formulation  of  the  problem 
indicated  excellent  equilibria]  existed,  and  agreement  of 
theoretical  and  experimental  results  at  the  rear  of  the 
wing  is  good,  the  leading  edge  divergence  was  likely  to  be 
the  result  of  an  initial  assumption. 

Attention  was  immediately  directed  to  the  root  boundary 
conditions.   It  was  assumed  that  the  pivot  rib  (considered 
to  be  the  root  of  the  idealized  model)  was  infinitely  stiff 
and  completely  resisted  chordwise  deformation.  This  is 
still  felt  to  be  a  reasonably  accurate  assumption.   However, 
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re -examination  of  the  drawings  of  the  pivot  rib  indicates 
that  the  rib  width  varies  by  descreasing  from  rear  to  front 
substantially.  Moreover,  the  skin  thickness  varies  the 
same  way.   It  would  therefore  appear  that  significant  twist 
about  the  streamwise  axis  occurs  towards  the  front  of  the 
pivot  rib.   This  has  the  effect  of  removing  the  root 
restraint  imposed  upon  the  theoretical  analysis,  and  allows 
stress  relief  to  occur  towards  the  leading  edge  in  the  actual 
case. 

The  sizeable  differences  in  back-up  gage  readings 
would  further  stimulate  interest  into  investigating  the 
plate  bending  and  twisting  effects  of  individual  panels 
about  their  respective  axes.  Unfortunately,  insufficient 
back-up  gages  were  installed  to  establish  any  definite 
conclusions,  but  examination  of  the  stresses  at  the 
available  back-up  locations  and  the  stresses  measured  by 
neighboring  rosettes  on  only  one  side  of  the  skin  evokes 
considerable  concern  in  this  area.   Panels  25  and  30 
would  be  particularly  well  suited  for  such  an  investigation 
in  the  future. 

It  is  therefore  believed  that  the  analysis  method 
can  achieve  accurate  results  for  this  wing.   The  boundary 
condition  in  question  could  be  removed  entirely  from  the 
analysis  by  moving  the  theoretical  root  to  the  aircraft 
center  line.   This  would  introduce  only  16  new  unknowns 
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to  the  analysis,  all  within  the  additional  rib  required 
at  station  y  =  25^095»   The  addition  of  16  unknowns  would 
in  turn  require  conversion  of  the  computer  program  to 
machine  language  in  order  to  accommodate  the  increased 
storage  requirements.   This  could  be  accomplished  by  use 
of  the  Fortran  MAP  program  available  at  the  U.  S.  Naval 
Postgraduate  School  computer  Facility, 

There  were  no  load  attachment  points  on  the  specimen 
corresponding  to  the  analytical  load  attachment  points. 
This  could  have  enabled  direct  experimental  verification 
of  a  greater  portion  of  the  flexibility  influence  coefficient 
matrix  |CJ  (Table  VIII),  However  two  node  points  can  be 
compared  here.  The  measured  net  vertical  deflections  at 
scale  numbers  8  and  11  (Fig.  26)  are  compared  to  deflections 
caused  by  the  analytically  calculated  vertical  loads  P* 
through  P„,  At  scale  number  8  the  theoretical  value  differs 
from  the  measured  value  by  only  16,8  percent  and  at  scale 
number  11  the  difference  is  only  1*f.8  percent.  The  better 
agreement  between  theoretical  and  measured  deflections  was 
expected(Ref .  9)» 

Improved  accuracy  in  the  solution  could  undoubtedly 
be  obtained  by  reducing  the  size  of  the  structural  grid. 
This  would  at  the  least  double  the  number  of  unknowns. 
With  the  additional  element  flexibilities  required, 
computer  programming  should  be  extended  to  do  more  of  the 
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labor.   As  has  been  pointed  out  by  Rattinger  and  Gallagher 
(Ref.  2),  up  to  three  man-years  may  be  required  for 
completely  programming  a  solution  using  a  displacement 
method.   It  is  felt  that  the  force  method  of  this  report 
could  be  extended  to  -twice  the  current  size  with  about  one 
man-month  additional  effort,  if  project  familiarity  were 
equal  to  that  of  the  writers. 
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V  CONCLUSIONS  AND  RECOMMENDATIONS 

1.  Conclusions 

It  can  therefore  be  concluded  that  a  valid  comparison 
was  achieved  between  experimental  and  theoretical  stresses 
as  predicted  by  the  matrix  force  method  of  analysis,  using 
a  structural  idealization  proposed  by  Wehle  and  Lansing. 
Further,  this  analytical  method  gave  remarkably  reasonable 
results  in  the  critical  rear  root  area  of  the  highly 
complex  structure. 

The  failure  of  the  method  to  agree  with  actual  stress 
conditions  in  the  forward  root  area  is  attributed  primarily 
to  the  simplifying  assumption  that  the  pivot  rib  was 
infinitely  rigid. 

2.  Recommendations 

It  is  recommended  that  all  cover  skin  rosettes, 
especially  those  between  the  intermediate  and  pivot  rib,  be 
backed-up.   This  would  greatly  enhance  the  experimental 
potentiality  of  the  laboratory  by  providing  means  to 
investigate  the  plate  phenomena  of  bending  and  twisting  of 
individual  panels  about  their  respective  axes. 

In  addition  greater  academic  value  could  be  realized 
from  the  analytical  results  if  means  were  provided  for 
applying  single  point  loads  at  the  intersections  of  the 
intermediate  rib  and  the  beams.   This  would  permit  further 
interesting  academic  demonstrations  by  utilizing  the  matrix 
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of  flexibility  influence  coefficients  produced  in  this 
analysis. 

It  is  further  recommended,  to  improve  the  accuracy 
of  the  analytical  predictions  along  the  entire  root  area, 
that  the  root  boundary  condition  assumed  in  this  analysis 
be  eliminated  by  moving  the  theoretical  root  to  the  wing 
center  line  rib.   This  would  not  require  substantial 
increase  in  analytical  complexity  and  would  greatly 
enhance  the  agreement  between  theory  and  actual  measurement. 

It  is  finally  recommended  that  the  experimental 
method  of  loading  be  improved  to  enable  better  development 
of  pure  torque.   In  this  regard  an  exterior  rib,  clamped 
around  the  section  perpendicular  to  the  center  intermediate 
rib  at  wing  station  yw  =  143 ,  would  eliminate  the  stream- 
wise  twist  of  the  wing  fold  rib.  Regardless  of  the  method 
of  loading,  the  plywood  pads  under  the  wing  support  jig 
should  be  replaced  with  a  more  suitable  material. 
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TABLE  II 
COMPONENTS  FOR  APPLIED  LOADS 


Cell 

<ls 
Shear  Flow 

Station 

81.98 

qssin9* 

Fp  = 
qsw  cos  9 

Ap  = 
Fp/cosA 

-(Ap  sLnAp) 

1 

265. 74 

99.84 

3581.83' 

4781.46 

3223.13 

2 

302.03 

113.76    . 

3790.74 

4869.70 

3056.91 

3 

398.78 

149.82 

4815.86 

5920.33 

3443.56 

4 

457. 81 ' 

171.99 

5335.69 

6330.00 

3405.79 

5 

472.81 

177.43 

5310.11 

6128.50 

3059.59 

6 

442.47 

166.23 

4840.75 

5476.53 

2561.21 

1 

*e»  (90- 2 A) 

(See  Appendix  C  for  resultant jpj  matrix) 
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TABLE  III 

• 

* 

THE  EDGE  SHEAR  FORCES  ON  THE  PANELS 

lii           >J 

Jl 

\  ' 

f 

ac     wc 

\ 

• 

1                     :*-■■ 

■■•d 

Piece 

wi 

wc 

ib 

Id 

ab/ai  = 

ac/ai  = 

ad/ai  = 

No. 

in. 

in. 

in. 

in. 

3-b/wc 

wi/wc 

ld/we 

1 

5.64 

6.82 

14.340883 

14.328741 

2.102769 

.826979 

2.100988 

2 

6.82 

7.44 

13.507743 

13.504185 

1.815557 

.916666 

1.815079 

3  - 

7.44 

7.66 

13.030764 

13.030300 

1.701144 

.971279 

1.701084 

4 

7.66 

7.54 

12.575610 

12.575467 

1.667853 

I.OI5915 

1.667834 

5 

7.54 

7.06 

12.137011 

12.129637 

1.719123 

1.067989 

1.718079 

6 

7.06 

6.30 

11.804608 

11.798490 

1.873747 

1.120635 

1.872776 

7 

5.64 

7.24 

79.3021 

79.3021 

10.953328 

.779006 

10.95333 

8 

6.82 

8.94 

74.9206 

74.9206 

8.380380 

.762864 

8.38038 

9 

7.44 

9.66 

71.3846 

71.3846 

7.389710 

.770186 

7.38971 

10 

7.66 

9.70 

68.5890 

68.5890 

7.071030 

.789691 

7.07103 

11     - 

7.54 

9.20 

66.4567 

66.4567 

7.223554 

.819565 

7.22355 

12 

7.06 

8.32 

64.9038 

64.9038 

7.800938 

.848558 

7.80094 

13 

6.30 

7.20 

63.8339 

63.8339 

8.865819 

.875000 

8.86582 

15 

14.328741 

20.870584 

79.3021 

74.9206 

3.799707 

.686557 

3.58977 

17 

13.504185 

19.099132 

74.9206 

71.3846- 

3.922723 

.707058 

3.73758 

19 

13.030300 

17.787945 

71.3846 

68.5890 

4.013088 

.732535 

3.85593 

21 

12.575^67 

16.511403 

68.5890 

66.4567 

4.154038 

.761623 

4.02490 

23 

12.129637 

15.239329 

66.4567 

64.9038 

4.360868 

.795943 

4.25897 

25 

11.798490 

14.141857 

64.9038 

63.8339 

4.589482 

.834295 

4.51383 

28 

11.804608 

14.174828 

64.9038 

63.8339 

4.578807 

.832787 

4.50333 

30 

12.132011 

15.245680 

66.4567 

64.9038 

4.359051 

.795767 

4.25719 

32 

12.575610 

16.513296 

68.5890 

66.4567 

4.153562 

.761544 

4.02444 

3^ 

13.030764 

17.787956 

71.3846    • 

68.5890 

4.013086 

.732561 

3.85592 

36 

13.507743 

19.102224 

74.9206 ; 

71.3846 

3.922088 

.707129 

3.73698 

38 

14.340883 

20.887886 

79.3021 

74.9206 

3.796559 

.686564 

3.58680 
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TABLE    IV 
Matrix     [A] 
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•OOOOCO 
•OOOOCO 

•OOOOCO 

.688888 

.000000 
.000000 
.000000 

.000003 

7.659698 
.000000 
.000000 
.000000 

.ooocco 

-7.459920 

7.539910 

.000000 

.000000 

1.000  190 

•OOOOOO 

-7.538520 

1.000320 

.000000 
-1.000000 

•OOOOOO 

•OOOOOO 

1 .000320 

.000000 

-1 .090990 

•OCCCOO 
•OCOCOO 
•OOOOOO 
•OOOOOO 
1.000]90 

006000 

OOOOOO 
OOOOOO 
OOOOOO 

999999 

.800888 

.000000 
.000000 
.000000 

•982222 

OOCOOO 
OOOOOO 
OOOOOO 
OOOOOO 

000909 

•OOOOOO 
•OOOOOO 
•OOOOOO 

•OOOOCO 
•OOOOCO 
.999120 

.000000 

.999770 

-  1.000000 

.999970 

-1.000000 

.000000 

-1.000000 
.000000 
•OOOOOO 

.000000 
.000000 

.000000 

.008800 

.000000 
.000000 

-1 .888008 

.000000  -1 
.000000 

999978 

OOOOOO 
OOOOOO 

.888888 

.999770 
-1.000000 

880000 

OOOOOO 
9991420 
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TABLE 

Mai  r  ix 

031 

. 

Inp 

ut  Form 

1  1 

1.0 

134 

0.0 

220 

1.0 

244 

-1.0      '       T 

3    3 

-1.0 

335 

1.0 

4    4 

-1.0 

436 

l.c            1 

5    5 

-1.0 

537 

1.0 

6    6 

-1.0 

6  38 

l.c            1 

7    7 

-1.0 

8    1 

55.2466 

8    2 

69.5754 

8    3 

83.C796       1 

8    U 

96. 1099 

8    5 

108.6853 

8    6 

120.8150 

8    7 

132.6135       1 

8    8 

7.4179 

8    9 

5.9191 

810 

4.2884 

811 

3.0187       1 

812 

1.8124 

813 

0.5448 

814 

0.5448 

815 

1.8124 

816 

3.0187 

817 

4.2884 

818 

5.9191 

319 

7.4172       1 

820 

2.5220 

821 

3.2013 

822 

3.5384 

823 

3.5707       1 

824 

3.3099 

825 

2.8758 

82£ 
830 

2.3642 

3-5707 

-15.9213 

827 

2.5220  1 
3.3099 

828 
832 

i-Ml 

829 

833 

3.5384 

831. 

2.3642 

834 

835 

-30.7843       1 

836 

-45.6107 

837 

-60.8687 

838 

-75.9532 

839 

-3.2013       1 

8U0 

-3.5384 

841 

-3.5707 

842 

-3.3099 

843 

-2.8758       1 

844 

-2.5220 

845 

-3.2013 

846 

-3.5384 

847 

-3.5707 

848 

-3.3099 

849 

-2.8758 

921 

1.0 

934 

8.3804       1 

9U5 

-1.0 

1022 

1.0 

1035 

7.3897 

1046 

-1.0              1 

1123 

1.0 

1136 

7.0710 

1147 

-1.0 

1224 

1.0              1 

1237 

7.2236 

1248 

-1.0 

1325 

1.0 

1338 

7.8009 

1349 

-1.0 

14    7 

0.0 

14    8 

7.0563 

1532 

-1.0               1 

1538 

-7.8009 

1543 

1.0 

1631 

-1.0 

1637 

-7.2236       1 

1642 

1.0 

1730 

-1.0 

1736 

-7.0710 

1741 

1.0              1 

1829 

-1.0 

1835 

-7.3897 

1840 

1.0 

1928 

-l.C 

193U 

-8.3804 

1939 

1.0 

20    8 

-0.99948 

20    9 

-0.99980     1 

2010 

-0.99999 

20  11 

-0.99996 

2012 

-0.99974 

2013 

-0.99915     1 

2014 

0.99915 

2015 

0.99974 

2C16 

0.99996 

2017 

0.99999    1 

2018 

0.99980 

2019 

0.9994  8 

2020 

-0.69666 

2021 

-0.65009     1 

2022 

-0.60391 

2023 

-0.55916 

2024 

-0.51788 

2025 

-0.48235    1 

2026 

-0.45373 

2027 

0.69666 

2028 

0.65009 

2029 

0.60391     1 

2030 

0.55916 

2031 

0.51788 

2C32 

0.48235 

2033 

0.45373    1 

2039 

-0.65009 

2040 

-0.60391 

2041 

-0.55916 

2042 

-0.51788     1 

2043 

-0.48235 

204U 

0.6966  6 

2045 

0.65009 

2046 

0.60391     1 

2047 

0.55916 

2048 

0.51788 

2049 

0.48235 

2120 

-0.71732     1 

1 

2121 

-0.75927 

2122 

-0.79688 

2123 

-0.82936 

2124 

-0.85597    1 

2125 

-0.87645 

2126 

-0.89114 

2127 

0.71732 

2128 

0.75927    1 

2129 

0.79688 

2130 

0.82936 

2131 

0.85597 

2132 

0.87645    1 

2133 

0.891  14 

2139 

-0.75927 

2140 

-0.79688 

2141 

-0.82936    1 

2142 

-0.85597 

2143 

-0.87645 

2144 

0.71732 

2145 

0.75927    1 

2146 

0.79688 

2147 

0.8293  6 

2148 

0.85597 

2149 

0.87645     1 

22    1 

56.885 

22    2 

56.835 

22    3 

56.885 

22    4 

56.885         1 

22    5 

56.885 

22    6 

56.885 

22    7 

56.885 

22    8 

1.83113    1 

22    9 

i. 12519 

2210 

0.27134 

2211 

-0.48011 

2212 

-1.30551     1 

2213 
2217 

-2.34025 

2214 
2218 

-2.34025 
1.12519 

2215 
2219 

-1.30551 
1.83113 

2216 
2220 

-0.48011  1 
2.59681     1 

0.27134 

2221 

3.39403 

2222 

3.84895 

2223 

4.02233 

2224 

3.93750     1 

2225 

»   3.64622 

2226 

3.2081 3 

2227 

2.59681 

2228 

3.39403    1« 

2229 

3.84895 

2220 

4.02233 

2231 

3.93750 

22  32 

3.64622     1 

2233 

3.20813 

2239 

-3.39403 

2240 

-3.84895 

2241 

-4.02233    1 

2242 

-3.93750 

2243 

-3.64622 

2244 

-2.59681 

2245 

-3.39403     1 

2246 

-3.84895 

2247 

-4.02233 

2248 

-3.93750 

2249 

-3.64622     1 

23    8 

56.8554 

23    9 

56.8736 

2310 

56.8844 

2311 

56.8827^     1 

2312 

56.8702 

2313 

56.8367 

2314 

-56.8367 

2315 

-56.8702       1  . 

2316 

-56.8827. 

2317 

-56.8844 

2318 

-56.8736 

2319 

-56.8554       1 

2320 

0.0 

2321 

-15.8462 

2322 

-31.8510 

2323 

-47.9019       1 

2324 

-63.5718 

2325 

-78.4493 

2326 

-92.3667 

2327 

O.C               1 

2328 

15.8462 

2329 

31.8510 

2330 

47.9019 

23  31 

63.5718 

2332 

78.4498 

2333 

92.3667 

2339 

-15.8462 

234C 

-31.8510       1 

2341 

-47.9019 

2342 

-63.5718 

2343 

-78.4498 

2344 

0.0              1 

2345 

15.8462 

23ti6 

31.8510 

2347 

47.9019 

2348 

63.5718       1 

2349 

78.4498 

2412 

6.8182 

2511 

7.4397 

2610 

7.6599      1 

27    9 

7.5385 

28    8 

1.0 

29    1 

1.0 

29    2 

1.0              1 

29    3 

1.0 

29    4 

1.0 

29    5 

1.0 

29    6 

l.C              1 

29    7 

1.0 

29    8 

0.03219 

29    9 

0.01978 

2910 

0.C0477    1 

2911 

-0.00844 
-0.02295 

29  12 

-C. 02295 

2913 

-0.04114 

2914 

-0.04114     1 

2915 

2916 

-0.00844 

2917 

0.00477 

2918 

0.01978    1 

2919 

0.03219 

2920 

0.01009 

2921 

0.01415 

2922 

0.01555     1 

2923 

0.01487 

2924 

0.01249 

2925 

0.00971 

2926 

0.00705     1 

2927 

0.01009 

2928 

0.0141 5 

2929 

0.01555 

2930 

0.01487     1 

2931 

0.01249 

2932 

0.0097  1 

2933 

0.00705 

2934 

-0.76286     1 

2935 

-0.77019 

2936 

-0.78969 

2937 

-0.81957 

29  38 

-0.84856     1 

2939 

-0.01415 

2940 

-0.0155  5 

2941 

-0.01487 

2942 

-0.C1249     1 

2943 

-0.00971 

2944 

-0.01009 

2945 

-0.01415 

2946 

-0.01555    1 

■' 

2947 

r0. 01487 

2948 

-0.01249 

2949 

-0.00971 

30    9 

1.0              1 

3018 

-1.0 

3110 

1.0 

3117 

-1.0 

3211 

1.0 

' 

3216 

-1.0 

3312 

1.0 

3315 

-1.0 

2 
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Table  VI 
Matrix    [Cj]33xlf9 

Seven  Columns  Per  Page 


1 

2 

3 

If 

5 

6 

7 

•598127 

-.327439 

-.257290 

-.189602 

-.124277 

-.061267 

.000022 

.401873 

.327439 

.257290 

.189602 

.124277 

.061267 

-.000022 

.665961 

.809153 

-.  146801 

-.104302 

-.063274 

-.023723 

.014746 

.596015 

.693787 

.785932 

-.154726 

-.097443 

-.042211 

.011512 

.566441 

.640150 

.709618 

.776648 

-.142992 

-.065500 

.009876 

.562138 

.621649 

.677734 

.731852 

.784081 

-.101880 

.009044 

.592316 

.649593 

.703573 

.755659 

.805927 

.854414 

.009229 

.663773 

.727960 

.788452 

.846822 

.903154 

.957490 

1  .010342 

-1.159430 

-.944685' 

-.742298 

-.547013 

-.358546 

-.176759 

.000065 

-1.060313 

-.863926 

-.678841 

-.500250 

-.327895 

-.161648 

.000059 

-.987361 

-.804486 

-.632135 

-.465832 

-.305335 

-.150527 

.000055 

-.916608 

-.746837 

-.586837 

-.432451 

-.283455 

-.139740 

.000051 

-.846247 

-.689508 

-.541790 

-.399255 

-.261696 

-.129013 

.000047 

-.786806 

-.641077 

-.503734 

-.371211 

-.243315 

-.119951 

.000044 

-.785232 

-.639794 

-.502726 

-.370469 

-.242828 

-.119711 

.000044 

-.846168 

-.689444 

-.541740 

-.399218 

-.261672 

-.129001 

.000047 

-.916801 

-.746994 

-.586960 

-.432542 

-.283515 

-.139769 

.000051 

-.987680 

-.804746 

-.632339 

-.465983 

-.305434 

-.150575 

.000055 

-1.060486 

-.864067 

-.678951 

-.500332 

-.327948 

-.161675 

.000059 

-1.158845 

-.944209 

-.741924 

-.546738 

-.358366 

-.176670 

.000065 

-.000000 

-.000000 

-.000000 

-.000000 

-.000000 

-.000000 

-.000000 

8.866033 

8.866033 

8.866033 

8.866033 

8.866033 

8.866033 

8.866033 

8.866033 

8.866033 

8.866033 

8.866033 

8.866033 

8.866033 

8.866033 

-.646568 

-2.2579  34 

-1.794824 

-1.347952 

-.916713 

-.500668 

-.095979 

-.497418 

-1  .65  1491 

-2.739149 

-2.036788 

-1.358953 

-.705097 

-.069092 

-.470755 

-1.275365 

-2.033670 

-2.765364 

-1.829741 

-.927268 

-.04  9437 

-.442605 

-.954483 

-1.436902 

-1.902393 

-2.351632 

-1.176225 

-.032913 

-.322946 

-.575915 

-.814326 

-1.044370 

-1.266384 

-1.480528 

-.017330 

.324519 

.577197 

.815333 

1.045112 

1.266870 

1.480768 

.017330 

.444256 

.955828 

1.437959 

1.903172 

2.352143 

1.176477 

.032913 

.472214 

1.276553 

2.034603 

2.766052 

1.830192 

.927491 

.049437 

.498556 

1.652419 

2.739873 

2.037325 

1.359306 

.705271 

.069092 

.647534 

2.258722 

1.795443 

1.348408 

.917012 

.500815 

.095979 

t-7-4 


[«]■■ 


8 

Q 

0 

10  " 

11 

12 

13 

1h 

•036899 

.037669 

.039542 

.o^ous 

.045770 

.051703 

.010636 

-.036899 

-.037669 

-.039542 

-.042045 

^-.045770 

-.051703 

-.010636 

.034379 

.034678 

.035642 

.037019 

.039139 

.024387 

-.022867 

.035882 

.036412 

.037889 

.039932 

.028518 

.015871 

-.026460 

.040567 

.041390 

.043551 

.033290. 

.022985 

.011318 

-.028936 

.049881 

.051159 

r~~  .039397 

.029728 

.019886 

.008558 

-.031406 

.062674 

.052079 

.039989 

.030006 

.019789 

.007952 

-.034162 

.070235 

.058362 

.044813 

.033626 

.022176 

.008912 

-.038283 

.106457 

.108677 

.114080 

.121303 

.132049 

.149167 

.030687 

.097356 

.099387 

.104328 

.110933 

.120760 

.136415 

.028063 

.090658 

.092549 

.097150 

.103301 

.112452 

.127030 

.026132 

.084161 

.085917 

.090188 

.095898 

. 104393 

.117927 

.024260 

.077701 

.079321 

.083265 

.088537 

.096380 

.108874 

.022398 

.072243 

.07  3750 

'   .077416 

.082318 

.089610 

.101227 

.020824 

-.064340 

-.062880 

-.059247 

-.054351 

-.047043 

-.035369 

.157176 

-.069333 

-.067760 

-.063845 

-.058569 

-.050694 

-.038114 

.169374 

-.075121 

-.073416 

-.069174 

-.063458 

-.054926 

-.041295 

.183512 

-.080928 

-.079092 

-.074522 

-.068364 

-.059172 

-.044488 

.197700 

-.08689*4 

-.084922 

-.080015 

-.073403 

-.06  35  34 

-.047767 

.212273 

-.094953 

-.092798 

-.087U36 

-.080211 

-.069427 

-.052198 

.231961 

-.764737 

-.764982 

.  -.765127 

-.76510U 

-.764936 

-.764486 

.764486 

.287108 

.177081 

.044001 

-.0731  19 

-.201765 

-.363040 

-.366458 

.902660 

.792830 

.  .659867 

.542729 

.413947 

.252310 

-.981807 

.183427 

.189783 

.203545 

.221340 

.247334 

-.711696 

.131170 

.343094 

.352019 

.372465 

.399356 

-.560988 

-.530593 

.117642 

.494713 

.506429 

.533983 

-.429501 

-.399895 

-.376442 

.098736 

.646520 

.661364 

-.303205 

-.278067 

-.257154 

-.239627 

.074731 

.810056 

-.171278 

-.152296 

-.1384  96 

-. 126648 

-.  116090 

.043173 

.053084 

.034691 

.015677 

.001871 

-.009961 

-.020463 

.093135 

.069U50 

.054833 

.019530 

-.005610 

-.026476 

-.043860 

.208500 

.061992 

.050451 

.022996 

-.0135U5 

-.04  3068 

-.066280 

.343761 

.042114 

.033310 

.012932 

-.013978 

-.053739 

-.083788 

.496563 

.017647 

.011353 

-.002375 

-.020182 

-,046239 

-.087331 

.667707 

7? 


[cj]  -  3 


15 

16 

17 

18 

19 

20 

21 

.004679 

.000945 

-.001560  ' 

-.003425 

-.004185 

.0546C4 

.01*3469 

-.004679 

-.000945 

.0C1560 

.003425 

*  .004185 

-.054604 

-.043469 

-.008142 

-.010273 

-.011651 

-.012606 

.054153 

.060804 

.057410 

-.01383B 

-.002433 

-.0C4U78 

-.005946 

.036489 

.054417 

.053040 

-.017294 

-.006997 

.003262 

.001109 

.027107 

.051717 

.051458 

-.020102 

-.010269 

-.0CC601 

.011169 

.021526 

.051323 

.051813 

-.022350 

-.012143 

-.002161 

.009937 

.020543 

.054078 

.054891 

-.025047 

-.013608 

-.002422 

.011136 

.023021 

.060603 

.061513 

.013498 

.002727 

-.004500 

-.009881 

-.012073 

-.105861 

.125411 

.012344 

.002493 

-.004115 

-.0,09036 

-.011041 

-.096811 

-. 140276 

.011495 

.002322 

-.003832 

-.008414 

-.010282 

-.090150 

-.130625 

.010671 

.002155 

-.003558 

-.007811 

-.009545 

-.083690 

-.121264 

.009852 

.001990 

-.003284 

-.007212 

-.008812 

-.077266 

-. 11 1956 

.009160 

.001850 

-.003054 

-.006705 

-.008193 

-.071839 

-. 104092 

.145616 

.138351 

.133461 

.129791 

.128262 

-.071677 

-.065837 

.156916 

.149087 

.143818 

.139863 

.138216 

-.0772U0 

-.070946 

.170014 

.161532 

.155823 

.151538 

.149753 

-.083687 

-.076868 

.183159 

.174021 

.167870 

.163254 

.161331 

-.090157 

-.082811 

.196660 

.186848 

.180244 

.175287 

.173223 

-.096803 

-.088916 

.214900 

.204178 

.196961 

.191545 

.189289 

-.105781 

-.097162 

.764936 

.765104 

.765127 

.764982 

.764737 

-1 .000056 

-.845400 

'  -.205185 

-.076540 

.040580 

.173661 

.283688 

.809454 

.956900 

-.820898 

-.692387 

-.575286 

-.442088 

-.331864 

.809454 

.705339 

.090120 

.064085 

.046284 

.032545 

-.233361 

-.059017 

.136139 

.087630 

.047890 

.020988 

.000602 

-.•145950 

-.045401 

. 100015 

.075564 

.046062 

.009533 

-.017927 

-.094124 

-.042970 

.059596 

.057389 

.036545 

.011416 

-.023888 

-.058456 

-.040403 

.024154 

.032708 

.020896 

.007100 

-.011911 

-.030289 

-.029481 

.001239 

.103704 

.115561 

.129370 

.148364 

.166700 

.029643 

.037004 

.226033 

.246955 

.272103 

-.692449 

.341841 

.040591 

.055083 

.367216 

.396829 

-.566650 

-.539078 

.536820 

.043161 

.064041 

.526955 

-.433390 

-.4C6492 

-.386028 

.760399 

.045587 

.071475 

-.291337 

-.265347 

-.247545 

-.233764 

1.032539 

.059212 

.086841 
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[cj]  -  k 


22 

23 

2k 

25. 

26 

27 

28 

.033989 

.024259 

.016002 

.007931 

.0CCC01 

.000001 

-.002875 

-.033989 

-.024259 

-.016002 

-.007931 

-.000001 

-.000001 

.002875 

"  .048118 

.038192 

.025181 

.012564 

.000833 

-.000002 

.006221 

.046621 

.036420 

.023823 

.011803 

.000650 

-.000001 

.007239 

.046558 

.038353 

.024920 

.012277 

.0CC558 

-.000000 

.007936 

.047802 

.040717 

.028631 

.014035 

.000511 

.000000 

.008626 

.051005 

.043961 

.031721 

.016898 

.0CC521 

.000000 

.009387 

.057158 

.049264 

.035548 

.018936 

.000584 

.000000 

.010519 

.098059 

.069989 

.046168 

.022880 

.0CC004 

.000004 

-.008294 

.089677 

.064006 

.042221 

.020924 

.000004 

.000003 

-.007585 

-.165678 

.059602 

.039316 

.019485 

.0CCC03 

.000003 

-.007063 

-.153806 

-.185426 

.036499 

.018089 

.000003 

.000003 

-.006557 

-.141999 

-.171192 

-.195711 

.016700 

.0CCC03 

.000003 

-.006054 

-.132025 

-.159168 

-.181964 

-.202870 

.000003 

.000002 

-.005628 

-.055302 

-.043870 

-.029019 

-.014187 

.0CCC03 

-.000016 

-.043634 

-.059593 

-.047274 

-.031271 

-.015288 

.000003 

-.000017 

-.047020 

-.064568 

-.051220 

-.033881 

-.016564 

.000004 

-.000018 

-.050945 

-.069560 

-.055180 

-.036501 

-.017845 

.000004 

-.000020 

-.054883 

-.074687 

-.059248 

-.0  39192 

-.019160 

.CCCC04 

-.000021 

-.058929 

-.081614 

-.064743 

-.042827 

-.020937 

.000004 

-.000023 

.214175 

-.682561 

-.51 1949 

-.338279 

-.166714 

-.0CCC04 

1.000056 

.845400 

. 1.048542 

1.093407 

1.094728 

1  .060437 

1.000027 

.000018 

.101080 

.703703 

.574792 

.406460 

.211091 

.0CCC07 

.000018 

.272641 

.081876 

.044594 

.023150 

.008466 

-.005421 

-.000010 

-.036391 

.258866 

.177913 

.111075 

.052196 

-.0C3902 

-.000010 

-.032660 

.172438 

.299429 

.190891 

.091747 

-.002792 

-.000008 

-.027402 

.096280 

.177961 

.268948 

.130310 

-.0C1859 

-.000006 

-.020722 

.036442 

.076770 

.122486 

.171152 

-.000979 

-.000003 

-.01 1957 

.040256 

.038491 

.030409 

.0174.71 

.000979 

-.000015 

-.026037 

.062795 

.061173 

.048327 

.026295 

.001859 

-.000032 

-.058223 

.075881 

.073921 

.0560  11 

.030209 

.0C2793 

-.000050 

-1095934 

.085651 

.080715 

.060050 

.032451 

.003903 

-.000071 

-.138535 

.098382 

.090858 

.066614 

.036122 

.005422 

-.000096 

-.186278 
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[cj]  -  5 


29 

30 

31 

32 

33 

3^ 

35 

-.003653 

-.003225 

-.002154 

-.001012 

.000001 

.667642 

.481468 

.003653 

.003226 

.002154 

.001012 

-.000001 

-.667642 

-.481468 

.007048 

.006581 

.004304 

.002289 

.000833 

-.275926 

.554455 

.009888 

.008097 

.0051 16 

.002595 

.000650 

-.188647 

-.368351 

.01 1660 

.01141* 

.007127 

.003516 

.000558 

-.142422 

-.279407 

.013178 

.013967 

.010959 

.005332 

.000511 

-.115168 

-.227115 

.014524 

.015766 

.013093 

.007721 

.000521 

-.110933 

-.219338 

.016277 

.017663 

.014673 

.008653 

.000534 

-.124316 

-.245798 

-.010538 

-.009303 

-.006215 

-.002921 

.000004 

1  .926192 

1  .389068 

-.009637 

-.008512 

-.005684 

-.002671 

.000004 

-.375191 

U270320 

-.008974 

-.007926 

-.005293 

-.002487 

.000004 

-.349377 

-.658480 

-.008331 

-.007353 

-.004914 

-.002309 

.000003 

-.324341 

-.611294 

-.007692 

-.006793 

-.004536 

-.002132 

.000003 

-.299444 

-.564370 

-.007151 

-.006316 

-.004218 

-.001982 

.000003 

-.278411 

-.524728 

-.083555 

-.121239 

-.156745 

-.190202 

.000003 

-.277603 

-.523376 

-.090040 

-.130643 

-.168909 

.029836 

.000003 

-.299146 

-.563992 

-.097555 

-.14155+ 

.065446 

.032326 

.000003 

-.324117 

-.611070 

-.105098 

.106843 

.070505 

.034826 

.000003 

-.349175 

-.658313 

.154708 

. 114719 

.075702 

.037393 

.000004 

-.374914 

1.270295 

.169057 

.12535? 

.082724 

.040861 

.000004 

1 .924836 

1.388115 

.682561 

.511949 

.338279 

.166714 

.000004 

.000000 

.000000 

.151244 

.160424 

.132661 

.076154 

.000004 

.000000 

.000000 

.496083 

.679040 

.820928 

.925501 

1.000024 

.000000 

.000000 

-.057173 

-.0614  33 

-.047151 

-.026497 

-.005421 

3.093957 

1.977363 

-.053981 

-.057958 

-.044993 

-.024997 

-.003903 

2.216014 

4.253180 

-.046122 

-.052093 

-.041572 

-.023065 

-.002792 

1  .545211 

2.967102 

-.035005 

-.040413 

-.034598 

-.019508 

-.001859 

.983284 

1.889707 

-.020057 

-.023213 

-.020309 

-.012482 

-.000979 

.486114 

.935407 

-.056322 

-.091663 

-.132169 

-.175678 

.000979 

-.485306 

-.934056 

-.123693 

-.198273 

r.  282192 

-.136654 

.001859 

-.982178 

-1.887978 

-.201806 

-.32080* 

-.204866 

-.098466 

.002792 

-1.543881 

-2.965149 

-.290150 

-.200233 

-.125686 

-.059240 

.003902 

-2.214480 

-4.251059 

-.122717 

-.073609 

-.042194 

-.017702 

.005420 

-3.092144 

-1.975265J 
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[cj]    -  6 


36 

37 

38 

39 

ko 

»f1 

k2 

.338325 

.224309 

.1152»40 

.002875 

.003653 

.003226 

.002154 

-.338325 

-.224309 

-.1152»40 

-.002875 

'-.003653 

-.003226 

-.002154 

.420912 

.276236 

.139517 

-.006221 

-.007048 

-.006581 

-.0043C4 

.469520 

.306469 

.154475 

-.007239 

-.009888 

-.008097 

-.005116 

-.424733 

.374470 

.188657 

-.007936 

-.011660 

-.01 1414 

-.007127 

-.345168 

-.498116 

.252915 

-.008626 

-.013178 

-.013967 

-.010959 

-.333312 

-.482222 

-.662407 

-.009387 

-.014524 

-.015766 

-.013093 

-.37  3523 

-.540398 

-.7U2320 

-.010519 

-.016277 

-.017668 

-.014673 

.976090 

.647145 

.332476 

.008294 

.010538 

.009308 

.006215 

.892646 

.591822 

.304054 

•  .007585 

.0096  37 

.008512 

.005684 

.831230 

.551104 

.283134 

.007063 

.008974 

.007926 

.005293 

-.930730 

.511612 

.2628U5 

.006557 

.008331 

.007358 

.004914 

-.859285 

-1.184810 

.242668 

.006054 

.007692 

.006793 

.004536 

-.798928 

-1.101588 

-1.478073 

.00  5623 

.007151 

.006316 

.004218 

-.797021 

-1.099057 

-1.474699 

.043634 

.083555 

.121239 

.156745 

-.858872 

-1.184347 

.242500 

.047020 

.090040 

.130648 

.T68909 

-.930564 

.511519 

.2627M2 

.050945 

.097555 

.141554 

-.065446 

.831290 

.551066 

.283055 

.054883 

.105098 

-.106843 

-.070505 

.892568 

.591687 

.303920 

.058929 

-.154708 

-.  114719 

-.075702 

.975353 

.646566 

.332108 

-.214175 

-.169057 

-.125359 

-.082724 

.000000 

.000000 

.ooooco 

-.84  54.00 

-.682561 

-.51  1949 

-.338279 

.000000 

.000000 

.  .ooooco 

-. 101080 

-. 151244 

-.160424 

-.  132661 

.000000 

.000000 

.ooooco 

-.272641 

-.496083 

-.679040 

-.820928 

1.228755 

.743506 

.360409 

.036391 

.057173 

.061430 

.047151 

2.884878 

1.836419 

.917554 

.032660 

.053981 

.057958 

.044993 

4.514167 

2.908448 

1.463261 

.027402 

.046122 

.052098 

.041572 

2.874906 

4.044536 

2.040715 

.020722 

.035005 

.040418 

.034598 

1.423008 

2.004500 

2.718373 

.011957 

.020057 

.023218 

.020309 

-1.421101 

-2.001970 

-2.7150C0 

.026037 

.056322 

.091668 

.132169 

-2.872586 

-4.041543 

-2.037511 

.058223 

.123693 

.198278 

.282192 

-4.51 1682 

-2.905548 

-1.460161 

.095934 

.201806 

.320804 

.204866 

-2.882332 

-1.833556 

-.914532 

. 138535 

.290150 

.200233 

. 125686 

-1.226285 

-.740776 

-.357518 

.186278 

.122717 

.073609 

.042194 
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[w]  -  7 


^3 


>4f 


.001012 

.001012 

.002289 

.002595 

.003516 

.005332 

.007721 

.008653 

.002921 

.002671 

.002487 

.002309 

.002132 

.001982 

.190202 

-.029836 

-.032326 

-.034826 

-.037393 

-.040861 

-.166714 

-.076154 

-.925501 

.026497 

.024997 

.023065 

.019508 

.012482 

.175678 

.136654 

.098466 

.059240 

.017702 


-.054604 

.054604 

-.060804 

-.054417 

-.051717 

-.051323 

-.054078 

-.060603 

.105861 

.096811 

.090150 

.083690 

.077266 

.071839 

.071677 

.077240 

.083687 

.090157 

.096803 

.105781 

1.000056 

-.809454 

-.809454 

.059017 

.045401 

.042970 

.0*40403 

.029481 

-.029643 

-.040591 

-.043161 

-.045587 

-.059212 


^5   1 

.043469 

.043469 

.057410 

.053040 

.051458 

.051813 

.054891 

.061513 

.125411 

.140276 

.130625 

.121264 

.1 11956 

.104092 

.065837 

.070946 

.076868 

.082811 

.088916 

.097162 

.845400 

-.956900 

-.785339 

-.136139 

-.100015 

-.059596 

-.024154 

-.001239 

-.037004 

-.055083 

-.064041 

-.071475 

-.086841 


»*6 


h? 


k8 


^9 


-.033989 

.033989 

-.048118 

-.046621 

-.046558 

-.047802 

-.051005 

-.057158 

-.098059 

-.089677 

.165678 

.153806 

.14  1999 

.132025 

.055302 

.059593 

.064568 

.069560 

.074687 

.081614 

.682561 

-1.048542 

-.703703 

-.081876 

-.258866 

-.172438 

-.096280 

-.036442 

-.040256 

-.062795 

-.075881 

-.085651 

-.098382 


-.024259 

• .024259 

-.038192 

-.036420 

-.038353 

-.040717 

-.043961 

-.049264 

-.069989 

-.0640C6 

-.059602 

.185426 

.171192 

.159168 

.043870 

.047274 

.051220 

.055180 

.059248 

.064743 

.511949 

-1 .093407 

-.574792 

-.044594 

-.177913 

-.299429 

-.177961 

-.076770 

-.038491 

-.061173 

-.073921 

-.080715 

-.090858 


-.016002 

.016002 

-.025181 

-.023823 

-.024920 

-.028631 

-.031721 

-.035548 

-.046168 

-.042221 

-.039316 

-.036499 

.195711 

.181964 

.029019 

.031271 

.033881 

.036501 

.039192 

.042827 

.338279 

-1.094728 

-.406460 

-.023150 

-.  11 1075 

-.190891 

-.268948 

-.122486 

-.030409 

-.048327 

-.056011 

-.060050 

-.066614 


-.007931 

.007931 

-.012564 

-.011803 

-.012277 

-.014035 

-.016898 

-.018936 

-.022880 

-.020924 

-.019485 

-.018089 

-.016700 

.202870 

.014187 

.015288 

.016564 

.017845 

.019160 

.020937 

.166714 

-1.060437 

-.211091 

-.008466 

-.052196 

-.091747 

-.130310 

-.  171 152 

-.017471 

-.026295 

-.030209 

-.032451 

-.036122 
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TABLE  VII 
CHECK  OP     CJ     EQUILIBRIUM 


- 

< 

Element 

Equation 

Error      > 

16 

3.58977(a2Q)   -  3.922723(alq)    + 

a39 

-.00000006 

18 

3-73758(a19)   -  4.013088(al8)    + 

a40 

+.00000018 

20 

3.85593(al8)   -  4.154038(a17)   + 

a4l 

+.00000028 

22 

4.02490(a1?)   -  4.360868(al6)    + 

a42 

+.00000033 

24 

4.25897(al6)   -  4.589482(a15)    + 

a43 

-.00000077 

29 

4.578887(al4)  -  4.25719(a13) 

a49 

+.00000053 

31 

4.359051  (a13)  -  4.02444(a12) 

a48 

• +.00000082 

33 

4.153562(a12)  -  3.85592(an) 

a47 

+.00000047 

35 

4.013086(an)  -  3.73698(a10) 

a46 

-.00000068 

37 

3.922088(a10)  -  3.58680(39) 

a/45 

--.00000026 

39 

3.796559(a9)   +  10.95333(a2) 

a44 

-.000001+30 

2nd  post 

.826979C-. 3274394)  +  1.0  -  .8091530 

-  0 

-.0799385 

1 

3rd  post 

.916666(-. 1468006)  +  1.0..  -  .7859316 

-  0 

+.0795013 

4th  post 

.971279C-. 1547262)  +1   -  .7766476 

-  0 

+.0730701 

5th  post 

1.015915(-.142922)  +1   -  .7840807 

-  0 

+.0707227 

6th  post 

1.067989C-. 1018803)  +1   -  .8544135 

-  0 

+.0367795 

7th  post 

1.120635(-. 0092,286)  +  1   -1,0103419 

-  0 

+  0.00000 
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TABLE    VIII 
ANALYTICAL     DEFLECTION     INFLUENCE      COEFFICIENTS 


TABLE   IX 
STRESS  COEFFICIENT  MATRIX    fSJ 


TABLE     X 
Matrix  {at| 


1  -92.32476270                   i      26  -1100.864  89153 

2  1592.32476283  27  -939.50998485 

3  -60.24970900  28  -859.546678C4 

4  -222.63401679  29  -3524.82157683 

5  -253.07605061  30  -3361.354731*30 

6  -94.44129715  31  -3227.876  56528 

7  442.2377934C  32  -2977.44333845 

8  -2286.41063917  33  -2714.08887357 

9  -3785.32614148  34  767.79220127 

10  -4119.10557437  35  9C4.957  34172 

11  -4764.02371061  36  592.9339C626 

12  -4985.30822517  37  -17.28751687 

13  -4969.72981596  38  -720.49754617 

14  -4492. 37862730  39  8592.00505352 

15  -4354.58987951  40  8725.85882282 

16  -4887.04344642  41  7870.6858331 C 

17  -5C13. 05519950  42  4969.65661299 

18  -4765.986935C2  43  1010.45895720 

19  -4219.25686240  uu  5460.77078927 

20  -3867.37382425  45  8681.69078517 

21  5143.74844420  U6  8356.9425C32 1 

22  200.58214441  U7  7980.72271180 

23  -379.38246968  J48  4504.14537561 

24  -1720.23525162  U9  769.32658458 

25  -1354.25763476 
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TABLE  XI 

RESULTS  AT  STATION  yw  =  3*f.  5" 

Comparison  of  Maximum  Normal  Stress,  Maximum  Shear  Stress 
And  Principal  Axis  Deviation  From  Y-Axis.   Left  Wing, 
Forward  Of  Y-Axis  Is  Positive.   Theoretical  Loading 
Corresponds  To  Experimental  Loading  Of 
336000  in.' lb  Torque 


8 

<D 


Panel  No. 

Maximum  Normal  Stress 

Maximum  Shear  Stress 

Principal 

Experi- 
mental 

Theo- 
retical 

Experi- 
mental 

Theo- 
retical 

Experi- 
mental 

Theo- 
retical 

15 

+89^ 

+2234 

+759 

1757 

+  3°36'    -5°15' 

17 

+2025 

1328 

-9°l8» 

19 

+1420 

525 

1025 

5°28> 

-iV  1 

21 

+780 

+1009 

690 

844 

+  7°8' 

+6°49' 

23 

-827 

723 

+10°29 ' 

25 

-576 

-994 

562 

706 

+10°56f 

+4°55' 

Rear  Beam 


CO 

u 


28 

+868 

+813 

800 

593 

+10°37' 

+6°9' 

30 

1 

+714 

596 

+0°11' 

32 

-835 

670 

+2°24T 

34 

-525* 

-1068 

525 

773 

+  5°l8f 

+I°l6f 

36 

-1575 

1075 

+0°08T 

38 

-800 

-1974 

535 

1493 

-15°19 ' 

+7°4l' 



*  Max.  compressive  stress  in  panel  3^ 
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,  •  TABLE  XII 

COMPARISON  OP  LINEARITY  OF  STRAIN  READINGS 
AT  VARIOUS  LOCATIONS  FOR  THE  7000  LB  AND  8000  LB  TESTS 


General  Location 
in  Wing 

Gage 
Number 

Strain  ( 

/J  in/in) 

Measured 
7000  lb 

8000 

1 
lb  Load 

Measured 

Extra- 

Differ- 

• 

Test 

8000  lb 

polated 

ence 

test 

from  7000 

Between 

lb  test 

Measured 
And  Extrap 

Upper  Skin 

Outboard  Near 

131 

-  52 

-  52 

-59 

7 

Leading  Edge 

132 

+123 

+147 

+141 

6 

133 

+  31 

+  44 

+  35 

9 

Outboard  Near 

**3 

-21 

-  22 

-23 

1 

Trailing  Edge 

44 

+  46. 

+  53 

+  53 

0 

45 

+  2 

+  4 

+  2 

2 

Inboard  Near 

16 

-  11 

-  12 

-  13 

,  1 

Leading  Edge 

17 

+  72 

+  83 

+  82 

1 

18 

+  28 

+  30 

+  32 

2 

Inboard  Near 

89 

+  8 

+  8 

+  9 

1 

Trailing  Edge 

90 

+  50 

+  53   x 

+  57 

4 

91 

-  11 

-17 

+  13 

4 

Lower  Skin 

Outboard  Near 

169 

+  55 

+  66 

+  63 

3 

Leading  Edge 

170 

+141 : 

+162 

+161 

1 

171 

-  36 

-  26.  • 

-  41 

-  41 

0 

Outboard  Near 

394 

-  30 

-  30 

0 

Training  Edge 

395 

+150 

+170 

+171 

1 

396 

+  40 

+  49 

+  46 

3 

Inboard  Near 

161 

+  12 

+  17 

+  14 

3 

Leading  Edge 

162 

+  10 

+  15 

+  11 

4 

Inboard  Near 

367 

+  13 

+  14 

+  15 

1 

Trailing  Edge 

368 

-42 

-  50 

+  48 

2 

- 

369 

+  21 

+  24 

+  24 

0 

FB  Web 

152 

-  7 

-  10 

-  8 

2 

153 

+157 

+176 

+179 

3 

154 

+  25 

+  25 

+  29 

4 

CIB  Web 

273 

+  2 

+  3 

+  2 

1 

274 

-  25  ' 

-  28 

-  29 

1 

275 

-  2 

-  1 

-  2 

1 

RB  Web 

409 

-  20 

-  22 

-  23 

l' 

410 

+134 

+155 

+153 

2 

111 

+  7 

+  5 

+  8 

3 

Averag 

b   Difference 

;  Between.  Mea 

sured  and  Ex 

trapolated 

=  2.3 

ftp 


TABLE  XIII 
EXPERIMENTAL  AND  THEORETICAL  SHEAR  FLOWS  AND  TOTAL 
TORQUE  REACTED  AT  THE  SECTION  PERPENDICULAR  TO 
CIB  AT  yw  =  98.7 


Cell 
Mo. 

Rosette 

Mo. 

i  Twice 
Cell 

Shear  Flow. a.( lb/in.) 

Torque  =  2Aa     (in.   lb) 

Fxp^rlme 

pfcai 

Thporet 

cal 

Exp. 

Theoretical 

Area 
(2A) 

Initial 
8000  lb 

Average 
8000  lb 

323,785 
in.   lb 

336,000 
in.   lb 

8000  lb 
test 

323,785. 
in.  lb 

336,000 
in.  lb 

in2 

test 

test 

load 

load 

load 

load 

Web 

152 

222* 

211 

206.2 

213.95 

1 

a 

131 
169 

86.2 

211 
211 

211 

206.2 

213.95 

18188.2 

17774.4 

18442.1 

Web 

181 
196 

224 
175 

200 

232.5 

241.35 

b 

46 
216 

105.8 

415 
415 

415 

438.7 

455.30 

43907.0 

46239.0 

48170.3 

Web 

225 
237 

48.8 
53.8 

51.3 

49.9 

51.78 

c 

49 
252 

127.6 

495 
481 

488 

488.6 

507.08 

62268.8 

62345.4 

64703.9 

Web 

273 

261 

10.9 
Bad  Gage 

10.9 

4.7 

4.83 

d 

52 

291 

140.0 

512- 
505 

509 

493.3 

511.91 

71260.0 

69062.0 

71667.4 

Web 

300 
312 

5.8 
7.8 

7 

17.8 

18.52 

i 

e 

55 
333 

143.4 

543* 
478 

478 

475.5 

493.39 

68545.2 

68186.7 

70752.4  i 

1 

; 
1 

Web 

342 
355 

35.5 

41.3 

38.4 

40.7 

'42.20 

• 

f 

58 
394 

138.0 

414 
450 

432 

434.8 

451.19 

59616.0 

60002.4 

62263.9  I 
i 

• 
• 

Web 

409 

405* 

432 

434.8 

451.19 

Total 

Reactei 

3  Torque  = 

l2Aq 

= 

323785.2 

323609. 9    [ 

336000.0 

*  Considered  Stray  Values 
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TABIE  XIV 
EXPERIMENTAL  RESULTS  AT  yw  =  34.5 


Pane] 
No 

Gage 

No. 

Strain 

headings 

Average 
Strains 

^max 

'  max 

$P** 

Interior 

Exterior 

Interior 

Exterior 

juin/in. 

/j-in/in. 

/-in/in 

psi 

psi 

15 

161* 
162 

486 
485 
484 

+12, 
+14.5 

-43.5 
+103. 

+52.5 

-15.8 
+103. 
+33.5 

+894 

759 

+37°34' 

21 

279 
280 
281 

459 
458 

457 

+14. 

+92. 

+2. 

+22. 

+95. 
-14. 

+18. 
+93.5 
-6. 

+780 

690 

+41006' 

25 

367 
368 
369 

468 
467 
466 

+14. 

-49. 
+24. 

-46. 

-95. 

+7. 

-16 
-72 
+15.5 

-576 

572 

-44°54' 

28 

28 
29 

30 

474 

473 
472 

-22. 

+136. 

+25. 

+28, 
+78. 
-13. 

+3 

+107 
+6 

+868 

800 

-440351 

3^ 

19 
20 
21 

480 

479 
478 

-6. 
+72. 
+26. 

-16, 

+64. 

+4. 

+11 

+68 

+5 

+555 

525 

-39°26' 

38 

13 
14 

15 

483 
482 
481 

-78. 

+59.5 

+35. 

-66, 
+59. 
+39. 

-72 
+59 
+37 

-800 

535 

+18°39' 

*   Two  gage  rosette. 

**  Angle  measured  from  sweep  angle  of  CIS.  Plus  angles  are  measured  in  the 
direction  of  the  diagonal  gage  of  the  interior  rosette. 
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Fig.  19 


F8U-3  Wing  With  Upper  Skin  Removed 


S5 


Fig,  20 


F8U-3  Wing  Mounted  Inverted 
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Fig.  21 
Loading  Frame  and  Linkage 
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Fig.  22 


Vickers  V-Line  Piston  Type  Pump  Assembly 
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Fig.  21* 
Instrumentation 

a.  Hydraulic  Pressure  Gage 

b.  Strain  Gage  Lead  Junction  Panel 

c.  2C  Channel  Switching  and  Balancing  Unit 

d.  Wheats tone  Bridge  Circuit 

e.  Power  Supply 

f .  Electronic  Counter  Assembly 
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'  Fig.  25 

EXPERIMENTAL  AND  THEORETICAL  ShEAR  FLOWS 
CUT  SECTION  PERPENDICULAR  TO  C.I.B.  AT  yw  =98.7 


EXPERIMENTAL  RESULTS* 
8000  lb  Test  (T  =  323785  in.  lb) 
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THEORETICAL   RESULTS 
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FRONT  AND  REAR  BEAM  DEFLECTIONS 

Left  Wing,  Mounted  Inverted 
8000  lb  Test  (T  =  323785  in.  lb) 
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APPENDIX  A 

THEORETICAL  DETERMINATION  OF  SHEAR  FLOWS 

AT  SECTIONS  PERPENDICULAR  TO  CIB  AT 

y  -  98o7  AND  y  =  7if03 
w  w 

The  theoretical  shear  flows  at  sections  perpendicular 
to  the  CIB  at  yw  =  9S07   and  yw=  7*+°3  were  calculated  for 
two  reasons •   First  the  shear  flows  at  both  sections  would 
he  used  to  establish  the  theoretical  loads  at  streamwise 
rib,  y  =  8t,980   Secondly,  the  shear  flows  at  the  section 
perpendicular  to  the  CIB  at  y  =  98*7   would  be  determined 
experimentally  and  comparison  with  theoretical  values 
would  facilitate  a  cheek  of  the  degree  to  which  pure 
torque  loading  was  achieved..   Both  cross-sections  are 
shown  schematically  in  Figo  A1 „   The  thicknesses  were 
average  mid-panel  values  and  the  web  heights  were  taken 
be Ween  the  mid-panel  of  the  upper  and  lower  skins 0 

A  box  structure  with  several  cells  wil.l  have  one  less 
redundant  than  the  number  of  cells c  In  this  case  there  is 
five  redundant  webs0  It  is  desired  to  write  six  equations 
in  the  six  unknown  shear  flows «  This  was  done  by  equating 
the  angle  of  twist  of  one  cell  with  the  remaining  five9 
which  gives  five  equations,  and  then  writing  an  equilibrium 
of  torsional  moments  equation0 

The  first  five  equations  were  obtained  by  equating 
the  angle  of  twist  per  unit  length,  Q  ,  of  one  cell  with 


9^ 


the  other  five  using  the  well  known  expression  for  a 
box  beam s 

.   _  Y'-q  As  L 
L_2   A  t  G 


The  equilibrium  of  torsional  moments  may  be  written 


as 


T  =   ^2Anqn 


(A2 


In  using  equation  (A2)  the  summation  is  carried  out 
around  the  entire  perimeter  of  each  cell0   The  unit  length 
L  and  the  constant  factor  2G  drop  out  leaving,,  for  cells 
a  and  b„ 


"q  As   _   \   q  As 


a 


Aa   t 


Ai 


The  value  of  q  for  any  exterior  web  of  cell  a  is  qat 
for  the  interior  web  is  (qa^  -  %t J  °   Using  the 
abbreviations , 

<5».  =  I  ^ 


(A3 


and 


aa 


6 
6 


bb 


ab 


a   t 
V  As 

As 

t 


a=~b 
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equation  A3  is  rewritten „ 

Vt  ^aa  -  ^bt  6ab  =  *bt  (5bb  .  Oat  6ab     (Ai+) 
Aa         Aa         Ab         Ab 

The  terms  0aa  andO^  represent  summation  around  the  entire 
perimeter  of  their  respective  cells  and  Q  ■.    the  value  for 
the  interior  web. 

For  each  cell  the  values  of  0  were  easily  calculated 
and  are  shown  in  Table  Al.   The  enclosed  areas  were  taken 
as  the  average  web  height  times  the  distance  between  webSo 

Equating  the  angle  of  twist  per  unit  length  of  cell 
a  to  the  remaining  cells  gives  the  following  equations? 

-L  [q*it  <5aa  "  qbt  Oabi 
Afl  / 

^bb  =  q-at    ^ab  °  ^ct     -;bc 


_  1 


^b 


qb 


Ac 


qct    (5ec   ~  <lbc    '-be  ~  qdt    ^cd 
qdt    ^dd   "  qct    *cd  "  <*et    ^de 


-   1 


(A5) 


Ad 


r     |  qet    5ee   -   qdt    C5d@   -   qft    5ef] 
I  qf t    <5ff  -   qet    rSef  J 


AQ 

A 


Substituting  values  of  6  from  Table  Alf,  performing  the 
indicated  arithmetic  and  rearranging^,  one  gets  five  equations 
in  terms  of  the  six  shear  flows  with  constant  coefficients „ 


c 


t". 


These  are  shown  in  Table  A2  for  both  sections u   The  sixth 
equation  is  also  shown  in  Table  A2,  in  which  case  the 
coefficients  are  simply  two  times  the  ceil  areas „ 

The  solution  to  these  six  simultaneous  equations  was 
obtained  on  the  CDC  1604  Digital  Computer  using  a  FORTRAN 
program  with  a  CO-OP  identification  of  F2  UTEX.  LINE  3N 
which  is  included  herein  as  pages  102  through  110.   It 
utilizes  Gauss ?s  method  of  elimination  with  row  pivoting 
and  back  substitution  and  is  designed  to  give  solutions 
for  one  or  more  column  vectors  forming  the  right  side  of 
the  set  of  equations,,   This  means  it  would  give  solutions 
for  one  or  more  values  of  applied  torque  in  this  case. 
Three  torque  values  were  used,  as  shown  in  Table  A3o   The 
values  of  294000  in,  lb  and  336000  in,  lb  were  used  in  the 
analysis o   The  value  of  177&25  ln*    Id  was  used  to  check 
the  computer  results 0   The  check  values  of  q  were  hand 
calculated  using  a  iterative  procedure  suggested  by  Bruhn 
(Ref0  11)  and  within  slide  rule  accuracy  agree  well  with 
the  computer  solutions . 


97 


Fig.  At 

SCHEMATIC  DRAWING  OF  SECTIONS 

PERPENDICULAR  TO  THE  C.I.B. 
AT  yw  «  9Q.7  AND  7*+.3 
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TABLE  -AI 
SECTION  PROPERTIES 


Section 

yw  =  98.7 

Section  yw  =  74. 

3 

Cell 

1  Area 

6nn= 

^W/l 

Area 

5nrT 

Umn" 

No. 

in2 

j  =  i   J 

in^ 

a 

43.1 

147.804 

27.940 

58.38 

333.606 

200.625 

b 

.52.9 

139.371 

67.600 

70.46 

455.417 

178.250 

c 

63.8 

250.743 

82.26 

369.016 

- 

133.60 

120.000 

d 

70.0 

324.737 

137.40 

84.39 

304.293 

119.206 

e 

71.70 

327.360 

132.40 

.  83.ll 

256.572 

78.556 

f 

69.0 

211.850 

75.11 

162.757 
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TABLE  AIII 
SHEAR  PLOW  RESULTS  FOR  VARIOUS  TORQUES 
AT  SECTIONS  PERPENDICULAR  TO  CIB  AT 
yw  =  98.7  ANDyw  =  7^.3 


Cell 

Shear  Flows   (lb/in.) 

Hand 
Calculated 

Computer  Results 

T=177,825 
in.  lb 

T=177,825 
in.  lb 

T=294,000 
in.  lb 

T=336,000 
in.  lb 

Yw=  7^.3 

a 

134.28 

135.289 

223.675 

255.628 

b 

166.09 

167.416 

276.790 

316.332 

c 

196.21 

197.292 

326.185 

372.782 

d 

222.05 

222.^50 

367.779 

420.319 

e 

229.65 

229.231 

378.990 

433.132 

f 

203.48 

201.905 

333.811 

381.500 

yw  =  98.7 

■ 

a 

' 

187.202 

213.945 

b 

398.384 

455.296 

c 

4^3.699 

507.084 

d 

447.921 

511.910 

e 

431.719 

493.392 

f 

• 

394.789 

451.188 

I  A.     IDENTIFICATION 

TITLE:    Solution  of  Simultaneous  Linear  Algebraic  Equations 

CO-OP  ID:    F2  UTEX  LINEQN 

CATEGORY:    Simultaneous  Linear  Equations 

PROGRAMMER:    C.   B.    Bailey 

DATE:    August  14,   1961 

B.  PURPOSE 
Solve  one  or  more  sets  of  linear  algebraic  equations  using  Gaussian 

elimination  with  row  pivoting  and  back  substitution. 

C.  USAGE 

1.  Calling  sequence: 
The  program  is  called  by  the  program  execute  card,    i.e.   LINEQN.  , 

in  the  normal  sequence  of  Fortran  control  cards. 

2.  Arguments: 
The  following  parameters  and  data  are  read  in  on  cards.     (See  9a). 

a.  N  -  the  order  of  the  matrix  in  the  equations  Ax  =  b,   that  is,  the 
number  of  linear  equations. 

b.  M  -  the  number  of  vectors  b  for  which  solution     vectors    x    are  to 
be  obtained,   that  is,   the  number  of  sets  of  linear  equations. 

c.  EP  -  Matrix  condition  parameter.     (See  Mathematical  Method.  ) 

d.  A  -  the  elements  of  the  matrix  of  coefficients  of  the  equations. 

e.  B  -  the  elements  of  the  column  vectors    b, ,   b_,    .  .  .  ,   b     , 

i        c  m 

3.  Space  required:     Undetermined. 

4.  Temporary  storage  required:    Space  is  reserved  for  solving  50  equa- 
tions in  50  unknowns  for  60  vectors    b.     This  amounts  to  8500  locations 

5.  Alarms  or  print-outs:    If  the  equations  are  inconsistent  or  dependent, 
MATRDC  SINGULAR  is  printed. 

6.  Error  returns:    None. 

7.  Error  stops:    None. 

8.  Input  and  output  tape  mounting:    Not  applicable. 

9.  Input  and  output  formats:    For  a  more  complete  description  of  the 
I,    E,    and  F  formats,    see  15. 
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Input (data  cards) 
(1)    First  card 

Columns  1-5      contain  the  value  of    N    inX5  format.     Thus    N    is 
is  punched  a  right- justified  fixed-point  integer. 
Leading  zeros  need  not  be  punched. 
Columns  6-10    contain  the  value  of    M    in  15  format. 
Columns  11-20  contain  the  value  of  EP  in  E10.4  format. 
To  enter  the  value  of    10         one  can  punch 
.1E-07  right-justified  in  the  field. 
(2)    The  coefficients  of  the  matrix  are  read  a  row  at  a  time.     Each 
row  is  begun  on  a  new  card.     Five  coefficients  per  card  are  in 
16  column  fields  (1-16,   17-32,    33-48,   49-$4,  -55-80)  using  F16.8 
format. 

The  matrix  of  coefficients  has  the  following  form: 

all  a12  •"  alN 

a21  a22  '  '  *  a2N 


aNl  aN2  "•  aNN* 

Thus,   for    N  =  6,  the  input  cards  would  contain  the  following 
information. 


• 

Columns 

Card  No. 

1-16 

17-32 

33-48 

49-54 

55-80 

2 

all 

a12 

>       a13 

a14 

a15 

3 

al6 

4 

a21 

a22 

a23 

a24 

a25 

5 

a26 

V 

6 

a31 

a32 

a33 

a34 

•  a35 

7 

a36 

12  a6l  a62  a63  a64  a65 


13  a66 
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(3)    The  elements  of  the  column  vectors    b.,   b~,    .  .  .  ,  b        are  also  read 

1       c.  m 

row  wise  in    F  16.8  format.     For    M  =  6,  the  vectors    b. ,   b.,,    ...,   b, 
are  assumed  to  be  arranged  in  the  following  manner: 

bll  b21  •••  b61 

b12  b22  "•  b62 


bl6  b26         •••  b66 


where    b. .    denotes  the    j-th    component  of  the  column  vector    b., 
Thus,  for    N  =  6  and    M  =  6,   the  data  cards  would  contain  the 
following  information: 

Columns 


Card  No. 

1-16 

17-32 

33-48 

49-54 

55-80 

14 

bll 

b21 

b31 

b41 

b51 

15 

b6l 

16 

b12 

b22 

b32 

b42 

b52 

17 

b,. 

24  b16  b26  b36  b46  b56 

25  b66 

(4)    Several  sets  of  data  may  be  processed  at  one  run.     To  terminate 

the  run,   place  one  blank  card  after  the  last  set  of  data- 
Output 

The  components  of  the  x  vectors  are  printed  row  wise  six  per 
line  in  E20.ll  format.  For  the  above  example  of  N  =  6  and  M  =  6, 
the  output  values  would  be  arranged  in  the  following  manner: 
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Line  No. 

Values 

1 

xll 

x21 

x31          x41 

x51 

X61 

2 

^2 

x22 

x32         x42 

x52 

x62 

b  X16  x26  x36         x46  z56  x66 

where     x..    denotes  the  j-th  component  of  the  column  vector    x, . 

10.  Selective  jump  and  stop  settings:    Not  applicable. 

11.  Timing:    Undetermined. 

12.  Accuracy:    Not  applicable. 

13.  Cautions  to  user:    None. 

14.  Equipment  configuration:    Not  applicable. 

15.  References: 

(a)  Fortran  System  for  the  Control  Data  1604  Computer,   Control  Data 
Corporation,   Computer  Division  Publication  08 7A,   Minneapolis, 
Minnesota  (1961). 

(b)  Fortran  Automatic  Coding  System  for  the  IBM  Data  Processing  System, 
International  Business  Machines  Corporation  (1958). 

(c)  Fortran  II  for  the  IBM  704  Data  Processing  System,  International 
Business  Machinea  Corporation  (1958). 

(d)  Kunz,   K.   S.  ,   Numerical  Analysis,   McGraw-Hill  Book  Co.  ,  Inc., 
New  York,   1957. 

(e)  Fadeeva,   V.   N.  ,    Computational  Methods  of  Linear  Algebra,  trans- 
lated by  Curtis  D.   Benster,   Dover  Publications,   Inc.,   New  York, 
1959. 

D.     MATHEMATICAL  METHOD 

Gauss's  method  of  elimination  with  row  pivoting  and  back  substitution  is 
used. 

Suppose  we  wish  to  solve 

Ax  =  b  and  Ax  =  c 

for  the  same  matrix  A.     This  routine  is  designed  to  give  the  solutions  for 
one  or  more  column  vectors    b,   c,   etc.    simultaneously. 
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Consider  a    3  by  3    matrix    A    augmented  by  two  column  vectors    b 
and    c. 


all 

a12 

a13 

bl 

cl 

a21 

a22 

a23 

b2 

C2 

a31 

a32 

a33 

b3 

C3 

The    A    matrix  is  triangularized  by  adding  multiples  of  one  row  to  another 
row.     For  column  1  find  the  row  containing  the  largest  element  in  column  1 
on  or  below  the  principal  diagonal.     Interchange  this  row  with  the  first  row. 

41 

Subtract  from  row    i  (i  =  2,    3)      times  row  1. 

all 


a 


•     - 


il 


=    a.,    -     - —  a,.  j  =  1, 


^    =    *ij    "     "^U   *lj 


This  eliminates  the  elements  in  column  1  below  the  principal  diagonal. 
This  process  of  pivoting  and  then  eliminating  elements  below  the  diagonal 
may  be  repeated  now  on  the  second  and  third  rows.     The  method  is  ex- 
tendible to  any  number  of  equations. 

Although  only  the    A    matrix  is  triangularized,  the  process  of  row 
interchange  and  addition  of  multiples  of  rows  is  applied  to  the    n    by    n+m 
matrix  consisting  of  the    n  by  n    coefficient  matrix  augmented  by    m 
column  vectors. 


where 
and 


Thus  we  have 

*11           ^                 ^3 

ell           e12 

0              d22                d23 

C21          e22 

o          o              d33 

e31          e32 

Back  substitution  is  used. 

e..  -  S 

xi  r      -       1J 

j  =  1,    ...,   M 
i  =  N,    ....   1 

l,J                    d.. 
ii 

(e 

S    a    0       for  i  =  N 

i 

N 

S    =     ^       dikXki 
k=i+l       liC    KJ 

for 

i  =  N-l,    .  .  .  , 

1. 


The  back  substitution  is  also  extendible  to  any  order  matrix. 
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V-1 


If  the  equations  are  either  inconsistent  or  dependent  (or  nearly  so), 
then  some  diagonal  element  of  the  triangularized  matrix  will  be  zero  or 
nearly  zero.  In  this  case  exponential  overflow  or  even  division  by  zero 
may  occur  in  the  back  substitution  phase.     The  alarm  print-out  (see  C.  5) 
is  set  off  when  a  diagonal  element  of  the  triangularized  matrix  is  less 
than  or  equal  to  argument  EP  (see  C.2.c).     Thus,   EP    is  that  value  which 
the  user  is  willing  to  call  zero.     It  should  be  larger  than  the  round-off 
error  which  can  occur. 
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PROGRAM  LINEQN 

DIMENSION  A(50,1  10),X(50,60) 

1  READ  101,N,MtEP 

101  F0RMAT(2I5,E10.4) 
IF(N)99,99,10 

10   PRINT  105 
105   FORMAT( 15H1PR0GRAM  LINEQN5X20HGAUSSI AN  ELIMINAT IDN// ) 
NPM  =  N+M 
NPO  =  N+l 
DO  2  K=l ,N 

2  READ  102, t A(K,J) ,J*1,N) 
DO  3  K=1,N 

3  READ  102, (A(K, J) ,J=NPO,NPM) 

102  F0RMAT(5F16.8) 

CALL  GAUSS2(N,M, EP,A,X,K1) 
GO  TO  (60,50) ,K1 

60  DO  61  K=l ,N 

61  PRINT  104, (X(K,J ),J=1,M) 
GO  TO  1 

104   F0RMAT(/(6E20.11 )) 
50    PRINT  103 

103  F0RMATM6H  MATRIX  SINGULAR) 
GO  TO  1 

99    STOP 
END 

SUBROUTINE    GAUSS2( N, M, EP, A, X, KER  ) 
DIMENSION  A(50,l 10),X<50,60) 
NPM=N+M 

10  DO  34  L=1,N 
KP  =  0 

Z  =  0.0 

DO  12  K=L,N 

IF(Z-ABSF(A(K.L) ) )11,12,12 

11  Z=ABSF(A(K,L) ) 
KP=K 

12  CONTINUE 
IF(L-KP) 13,20,23 

13  DO  14  J=L,NPM 
Z=A(L,J) 

A(L, J)=A(KP,J> 

14  A(KP,J)=Z 

20   IF(ABSF( A(L,L) )-EP )50, 50 ,30 

30  IF(L-N)31,40,40 

31  LP1=L+1 

DO  34  K=LP1,N 
IF(A(K,L) )32,34, 32 

32  RATIO=A(K,L)/A(L,L) 
DO  33  J=LP1,NPM 

33  A(K, J)=A(K, J)-RATIO»A(L, J) 

34  CONTINUE 

40  DO  43  1=1, N 
I I=N+1-I 

DO  43  J=l ,M 

JPN=J+N 

S  =  0.0 

IF( I I-N)41 ,43,43 

41  IIPlaII+1 

DO  42  K=IIP1,N 

42  S  =  S+A( II ,K)»X(K,  J) 

43  X(  II  ,J)  =  (A(II,JPN)-S)/A(II,II) 
KER=1 

RETURN 
50  KER  =  2 
END 
END 
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APPENDIX  B 

Di/.      &TION  OF  INDIVIDUAL  MEMBER,  _ 

LEXIBILITIES  REQUIRED  FOR  MATRIX  [f] 

I„   Dc ■  •   |  m      I   f  Ec 

The  symmetri<  matrix  |F|  of  equation  (2)  is  obtained 
bj   formulating  the  strain  energy  in  term-  of  internal 
ces.   Foi  example ,  the  an   tapered  bar  (A-j  -  A2)  of 
Figo  BJ  is  attached  to  a  web,  and   I  "    j     7  an  external 
axial  Jo ado  F„  and  an  internal  axial  load,  a-j0 


Fig.  Bl 
The  strain  energy  of  such  a  tar  may  be  expressed 
(Refe  B)    in  terms  of  the  end  loads  as 


r 
U   =   :  p(y)  dy  (Bla) 


Jr       2A.E 


For   a  load   equal   t< 


p(y)     =     u     +     (v  •=  u] 


y 


the  strain  energy  becomes 


U 


3AE 


p  p 

u  +  uv  +  v  )    (B' 


Upon  placing  the  bar  of  FIga  Bl  In  equilibrium  the 
end  loads  may  be  expressed  in  terms  of  the  internal  forces 
as 


u  —  "<^  a.,   t   a<5 
v  =   0  a*   +   a0 


(B2a) 

(B2b) 


By  substituting  equations  (B2)  into  (Bl)  the  internal 
strain  energy  can  be  written  as 

2 


2  EU 


3A 


2ai2  -6  ia^  +  3a3«^)  (B3) 


Now  referring-  again  to  equation  (2) 


2  EU  —  e^ • o  j  oS^ 


F 

31  F33. 


SI  i    o  o  o  OL  q 


1 1  2 


it  can  be  shown  by  taking  partial  derivative  of  the 
energy  that  the  F  matrix  terms  are  merely  coefficients 
of  squared  and  cross  terms  such  as 


2  EU  =  a,2  F„  +  F   V  +  F^a^  ♦  a^  (2a) 


Due  to  the  symmetry  required  of  |F   /Fn,  is  exactly 
one-half  the  coefficient  of  the  cross-term,  a-^a^o 

Three  factors 9  fxjs  ^l?3  an(*  ^22  are  now  introduced, 


11 


f- 


I 
3A 


6A„ 


f, 


22 


where  £■*-.    and  f?«,  will  refer  to  the  squared  terms,  and 


f1?  refers  to  the  cross-terms 


The  coefficients  of  a.,  and  a^  in  equations  \B2) 
are  now  identified  by 


(B5) 
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so  that  the  member  flexibilities  can  be  found  from 


F11  £  W  *  2  f12°^  *  f22J3 

f33  -   fl1(x   *   2f12qf^   +f22/^ 

F13  =  flAa/f22^4^f12CA°^^ 


^  (B6; 


Using  the  values  of   and   from  equations  (B2) 


P     s   If  I 

■11    ™ 

3A 


A3 


F„ 


I2 


=  F 


P 


33 


A 

1 
3A 


31 


II o   Application  to  Tapered  Bars 

For  linearly  tapered  bars  the  coefficients  f^i*  f^p: 
and  f.„  were  modified  by  the  functions  of  A3/A29  $ii? 
0  2,2?  and  0  22?  according  to  the  method  shown  in  Ref0  4? 
giving? 


f«      =    f      0 

11  1111 


ft    s   f    0 

12     12  12 


fU       B     f       0 

22     22  ^22 


B? 


11 


Equations  (B?)  are  then  used  directly  in  equations 


(B6) 


Consider  the  tapered  bar  shown  as  element  number  29 
in  Figo  13 5  and  re-drawn  below  in  Figo  B2Q   The  bar 


<nA- 


\05& 


A 


^^        q2 


q... 


*&**> 


f. 


,  -  3o8?  in< 


A 


A, 


A2  =  8o36  in? 


=  0V63 


Figo  B2 

geometry  is  obtained  from  Table  I0   Using  the  ratio 
Al/A2*  t^ie  P  functions  are  obtained  from  Fig0  4  of  Ref  0  4: 
giving 

f    »   ^11 

11 


12 


22 


3A-, 

VoV}87 

l012        = 

1 o 8001 

6A1 

f022        . 

3,0188 

3A- 


From  the  loading  system  upon  the  bar,  as  shown  in  Fig0  13 
the  alpha  (CX)  and  beta  {fi)   equations  are  respectively- 

ql  «  -  ho?055   a   +  >fo1502  a1lf  +  7o8009  a  g  -  a^ 
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M83  a       -  oV286  a1l( 


%9 


Substituting  into  equation    (B6),   ¥\>>   ^   becomes 


2 


F  -  kMS7i^o?055)  +  1  •  8001  (-V0  7055)  (  =  «M83) 

~?         +  C3o0l88){-o^83)2  =       102,6^85 


XII 0  Application  to  Symmetrical  Shear  Panel 

For  the  case  of  the  symmetrical  shear  panel ,  shown 
in  Figa  B3,  the  strain  energy  can  be  written  as  in  equation 
(B$): 


Fig.  B3 
Symmetrical  Shear  Panel 
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2EU 


,  '  ■■  •  J  I  S\     w         »■  (  1  H 


ri 


-/  t3 


2 


element,  number  13  s  Figo  3 


F    =   'tE 

8a8 


63o8339(13«>50) 

51  o^f 


|l/6(^|S2_)2  +  U32 


1o0! 


or 


E 


F 
898J 


E{  109o7173) 


IV„   Application  to  Skev/ed  Shear  Panels 

Skewed  panels  are  treated  according  to  the  method 
of  Garvey  (Ref0  10),  and  as  used  in  Ref„  6„ 


Fig,  Bh 


11? 


The  energy  of  the  swept  panel  shown  in  Fig,  B4  may  be 

wri  •  i  as 


2  EU 


I  i   )(  a.|  +  a2)  sin'  .  h  + 
"  It 


3TH 


Qi 


2 


^    (B9) 


where   gamma  is  wrj  as 


?  :  /a  .  2 

f  cot  o  cot  P  +■  cot 


Take  for  example  element  number  3#«   There  is 


T? 


D 


999 


giving 


F 


9,9 


io32(35e2288K0576I0(  1  *  .505C5  x  2.51a 


9 . 95-9 


(3 o 5868) 


9*9 


E 


C  78o5635) 


!  .  should  be  noted  that  for  both  the  bars  and  the 
panels,  trie  terms  of  the  [F  matrix  are  usually  made  up 
from  flexibilities  of  several  elements.,   In  addition,  all. 
main  diagonal  terns  (F  <)   are  posita   -   All  non-zero  terms 
of  the  flexibility  matrix  |Fj  are  given  in  the  computer 
pr     i  data  Input s  Appendix  C, 
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APPENDIX  C 

/OUTRAN  PROGRAM  FOR  CDC  1604 

DIGITAL  COMPUTER 

iES,  SYMBOLS,         I  .aTIONS 

A      Matrix  of  coefficients  of  determinate 

loads,,  basic  input  data„  size  33x33 
AI      Inverse  of  M    k  A.,  obtained  by  G3uss 

3  s  i  n  gl  e  -  pr  e  c is  i on  s  u  br-o  ut in  e 
AR.      Matrix  of  redundant  loads  f  output  s 

size  33x1 
AT      Matrix  of  all  loads s  output s   size  49x1 
B       Matrix  of  coefficients  of  applied  and 

redundant  loads,  basic  input,  size  49x33 
C       Matrix  of  deflection  influence  coefficients 

out put  5  size  33x33 
CJ      Matrix  product  of  (AI)(B),  size  33x49 
EN      Negative  of  matrix  product  (H22I)(H2l)s 

size  16x33 
ENP     Matrix  product  of  (H22I)(H21),  size  16x33 
G         r)  matrix  augmented  by  trivial  relations 9 

size  49x49 
GT      Transpose  of  [G)   matrix 
H      Matrix  triple  product  of  (GT)(UC)(G),  size 

49x49 
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Mali    i  first  33  columns  and  rows  of  (H), 
bained  d.    .  part      n   size  3  3x33 
Last  16  columns  (through  first 

3  3  rows)  of  (H)  matrix s  size  33x16 
H21    Matrix  of  last  16  rows  (through  first  33 

columns)  of  (H)  matrix,  size  16x33 
H22    Remainder  of  (H)  matrix,  size  16x1.6 
H22I   Inverse  of  (H22)  matrix 
H12N   Matrix  product  of  (H12)(EN),  size  33x33 

Gauss  3  subroutine  error  flag,  where 

assigned  value  of  2  indicates  singular  matrix 
LIST   Numbering  index 

P  Column  matrix  of  applied  loads,  size  33x1 
S  Matrix  of  stress  coefficientss  size  49x33 
UC     Matrix  of  flexibilities,  referred  to  in 

text  as  (F)  matrixs  basic  input  data. 

size  49x49 
UNITN  Matrix  formed  by  placing  unit  diagonal 

matrix  on  top  of  matrix  (EN),  giving 

matrix  saze  of  49x33 
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PROGRAM  ENERGY 


Start 


Clear  Arrays 

I   z 


Form  A,  B,  UC  Matrices 


Call  Gauss  3 


Return 
■ 1 


Yes 


-^Ja 


Print 
trixS  ingular 


) 


c 


Write  Out 
The  A I  Matrix 


[  Invert  A  Matrix 
r/f [Form  AI 


i 


Call  M&TMPY 


Form 
CJ  Matrix  =  A-1B 


I 


Form  G  Matrix  From  CJ  By 
Adding  Trivial  Relations 


( 


Write  Out 
The  G  Matri 


> 


Transpose  G 


I 


Form  H 


/   Write  Out^V 
V^The  H  Matrix/*" 


Form  H21 


i. 


Form  P  | ^vQ 

(     Write  Out   \ 
V^he  H22  Matrix  J 


Partition  H 

t — ; — ; — r~ 


Form  H22 


Form  H11 


Form  H12 


© 


(z> 


©- 


(^ 


Return 

zrz 


Call  Gauss  3 


Print 
atrlx  Singula 


Invert  H22 


Form  ENP 
ENP  =  H22I  x  H21 


Form  H22I 


EN  =  -ENP 


I 


\2 


Wt*  i  "fcf^ 
?he  H22I  Matrix 


D 


UNITN  =  I/EN 

5 


Form  S 
S  =  G  x  UNITN 


-fi 


Write  Out 
The  S  Matrix 


3 


3 


Form  AR 
AR  =  EN  x  P 


/     Write  Out  ^\ 
\The  AR  Matrix  J 


Form  AT 
AT  =  S  x  P 


3 


O 


f     Write  Out     ^v 
"VThe  AT  Matrix  y 


Form  C  ^ 

C'=  1/E(H11    +H12I  EN)   Hjg 


Write  Out 
he  C  Matrix 


; 


END 


PROGRAM  ENERGY 

COEFFICIENT  MATRICES  A,  E,  AND  UC  ARE  BASIC  STRUCTURAL  DMA 
OC I  MENS  ION  A(5C,50),AI(50,5C),B(5O,5O)tCJ(50,5O) ,UC(5C,50), 
lG(50,5C),H(50150)fGT(50t5C) ,H22 (50,50  )  ,H  1  2 ( 50 , 50 )  ,UCG ( 50  .  5C  ) . 
2H21(50,50),h22l(5C,50),EN<16,33),  AR(  16),  UNITN(  50,50  ).S(50,50), 
3AT(5C),H12N(50,50),C(50150) ,LIST(50),P( 35), ENP( 16,33) 
0ECUIVALENCE(A,CJ,UCG1H2  2,H2WS)t(B,G«Hl2,C),(UC,F), 
1(AI,GT,H22I.UNITN,H12N5 

CLEAR  ARRAYS  TG  RECEIVE  INPUT 

CO  10  1-1  ,5C  ' 

CO  1C  J=1,5C 

At  I,  J)  =  O.C 

E( I  ,J)  =  C.C 

LCI  I, J)  =0.0 
10  AI(  I, J)  =  0.0 

REAC  NCN-RECLNCANT  MATRIX  A 
31  OREAD  21,  WJ,(A(  I,J)),I1  ,J1,(A(I1,J1))»I2,J2,(A(I2,J2)), 
113,  J3, (A(  13, J3) ) .NEXT 

21  FORMAT  (M2I2,F12.6),I1  ) 
GC  TC  (31, 32), NEXT 
READ  LOAO  AND  REDUNDANT  MATRIX  B 

320REAC  22,I,J,(E(I,J)),I2,J2,(B(I2,J2)),I2,J2,(B(I2,J2)), 
113, J3, (BU3.J3)  ),NEXT 

22  FCRMAT(U(2I*.F12.6),I1) 
GO  TC  (32, 33), NEXT 

♦C      FORM  ENERGY  COEFFICIENT  MATRIX  UC 

330READ3U,I,J, (UC ( I . J ) ) , 1 1 , J 1 , (UC ( 1 1 , J  1 ) ) , 12, J2 ,  (UC ( 12  ,  J2  ) )  , 

113, J3, (UC(I3,J3)),NEXT 
3*4  FORMAT  (M2l2,F12.6),I1  ) 
GO  TO  (33.35),NEXT 

35  DO  36  I  =1,li9 
DO  36  J  =1,»i9 

36  UC( J, I )  =  UC(l, J) 

C      INVERT  THE  fi    MATRIX 

CALL  GAUSS3(33,1.E-10,A,AI,.KER) 
GO  TO(365,36«4)  ,KER 
36U  PRINT  151 

151  F0RMAT(2UH     MATRIX  A  IS  SINGULAR) 
365  PRINT  200 

37  PRINT  38 

380F0RMATI6UH  A  MATRIX  INVERSE  BY  GAUSS3t  9  COLUMNS,  IE  1-9,  10-18 
1PER  PAGE//)  ,. 

PRINT  100,  (  (Aid  ,J)  ,J=1  ,9),  1  =  1,33) 
PRINT  200 

PRINT  100, (  (AI(  I,J),J=1 0,18), 1  =  1,33) 
PRINT  200  •    .  ,,. 

PRINT  100, ( (AI( I, J) ,J=1 9,271,1=1,33) 
PRINT  200 

PRINT  300, ( (AI( I , J ) , J=28, 33 ) , 1*1 ,33 ) 
100  F0RMAT(/9F12.7) 
150  F0RMAT(/11F1C6) 
200  FCRMAT(lHl) 
C      POST  MULTIPLY  AI  WITH  B  MATRIX 

CALL  MATMPY(33,33,U9,AI ,B,CJ) 
C      FORM  G  MATRIX  BY  ADOITICN  OF  TRIVIAL  RELATIONS  TC  CJ 
CO  MC  1=34, M9 
UO  CJ(  I,I)=  l.C 
DO  Ml  1  =  1, I4S 
DO  Ul  J=1,MS 
U1  G(IiJ)  =  CJ(I,J) 
PRINT  200 

U3  F0RMAT(50H   G  MATRIX  .FCRMED  FROM  CJ  WITH  TRIVIAL  RELATIONS//) 
PRINT  100, ( (G(I»J)t J=1,9),I=1,U9) 
PRINT  200 

PRINT  100, ( (G(I,J),J=10,ie),I*1,U9) 
PRINT  200 

PRINT  100,  (  (G(I,J),J=19,27),I  =  1,I«9) 
PRINT  200 

PRINT  100, ( (G(ItJ)t J=28,36)tI*1tU9) 
PRINT  200 

PRINT100i((G(I,J)#J  =  37,ll5),I*1tU9   ) 
PRINT  200 


300 
U00 


i»U 


U5 


500 
600 


510 
57 
58 


581 


5U 

521 
522 


524 
525 

53 


5U1 
540 


56 

99 
62 

63 


PRINT 
FORMA 
FORMA 
TRANS 
CO  Ml4 
CO  UU 
GT(  1, 
FORM 
CALL 

PRINT 

PRINT 

FORMA 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

FORMA 

FORMA 

READ 

CO  51 

P(  I  ) 

REAC 

FORMA 

GO  TG 

PARTI 

CO  54 

CO  5M 

II  = 

JJ  = 

H22U 

PRINT 

PRINT 

FORMA 

PRINT 

PRINT 

INVER 

CALL 

GO  TO 

PRINT 

PRINT 

PRINT 

PRINT 

CC  53 

CO  53 

KK  = 

H21  (K 

FORM 

DO  5U 

00  5U 

SUM  = 

CO  5M 

SUM  = 

ENP(  I 

CO  56 

CO  56 

EN(  I, 

FORM 

CO  99 

CC  99 

UNITN 

DO  62 

UNITN 

CO  63 

CO  63 

KK  = 

UNITN 


400, ( 
T  (/  6 
T  (/  4 
POSE  G 

1  =  1, U 

J=1,4 
J)  =  G 
H  MATR 
MATMPY 
MATMPY 

200 

M5 
T(UUH 

500,  ( 

200 

500, ( 

200 

500, ( 

200 

500,  ( 

200 

500, ( 

200 

600, ( 
T  (/9F 
T  (/  4 
LOAD  M 
0  I  = 
=  0.0 
58,  I. 
T  (  12 

(57,5 
TION 
1  =  34 
J  =  34 
1-33 
J-33 
ItJJ) 

200 

522, ( 
T(/8F1 

200 

522. ( 
T   H22 
GAUSS 
(525,5 

150 

200 

522, ( 

522, ( 

I  =  34 
J=l, 
1-33 
K,J)  = 
REDUND 
0  1  =  1, 
0  K=l, 
0,0 
h   J=l, 

SUM  ♦ 

,K)   ; 
1=   , 
J=l, 

J)   =  - 

MATRIX 
I  =  1 
J  =  1 

(I, J) 

(1,1) 
1  =  34 
J- It 

1-33 

(I, J) 


(G(I,J),J=46,49),I*1,M9) 

F13.7) 

F13.7) 

9 

9 

(J, I) 

IX  BY  TRIPLE  PRODUCT  (GTMUCMG) 

<49,*»9,49,UC  ,G,UCG) 

(49f49t49,GT  ,UCG,H) 


MATRIX  H  FORMED  BY  TRIPLE 
(H(  I, J)  ,J  =  1,  9),  1  =  1,  »49) 

(H(1,J),J=10,18),I=1,49) 

(H(I,J),J=19,27),I=1,49) 

(H( I,J),J=28,36), 1*1,49) 

(H(I,J),J=37,45),I*1,49) 

(H(ITJ),J=46,49),Is1,49> 

13.  U) 

F13.4) 

ATRIX   P 

1,35 

(P(I)),NEXT 
•  F12.4,  I  1  ) 
81), NEXT 
F  MATRIX 
,49 
»M9 


PRODUCT,  GTXUCXG//) 


=  H(I,J) 

(H22(I,J),J  =1,8), I  =1,16) 
2.4) 

(H22(I,J),J  =9,  16), I  =1,  16) 

MATRIX 
2    (16,1.E-10,  H22,  H22I,KER) 
24),KER 

(H22I(I,J),J=1,8).I=1,16) 
(H22II  I,J),J-9,16)f'I*1t16) 

H(I,J) 
ANT  LOAD  COEFFICIENT  MATRIX,  EN. 
16 
33 

16 

H22K  I,J)*H21U,KI 
SUM 
16 

lie  '•  ENPU.J) 

UNITN  BY  PLACING  EN   BF.LOW   UNIT  MATRIX 
,33 
,33 
=  0.0 
1,33 
=  1.0 

»  EN(KK,J") 


124, 


800 


81 
80 

59 


600 
1 

700 


641 
640 


9000 
1 

1000 
52 


95 
90 


70 


730 
1 


FORK 
CALL 
PRINT 
PRINT 
FORMA 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
RECUN 
CO  80 
SUM  = 
DO  81 
SUN  = 
AR(  I 
CO  59 
LIST( 
PRINT 
PRINT 
FORMA 
16H 
PRINT 
FORMA 
TOTAL 
DO  64 
SUM  = 
CO  64 
SUM  = 
ATU  ) 
PRINT 
PRINT 
FORMA 

16H 
PRINT 
FORMA 
CO  52 
DO  52 
H12  ( I 
E  =  1 
INFLU 
DO  90 
CO  90 
SUM  = 
CO  95 
SUM  = 
H12N( 
SCALE 
CC  70 
DO  7C 
C(  I, J 
PRINT 
PRINT 
FORMA 
35H 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
ENC 


II 


S  MATRIX  . 
MATMPY  (M9 

200 

800 
T  CU3H 

150, ( ( 

200 
150*  t  (SUf 

200 
150,((S(I, 
DANT  LCAD 

I  =  1,16 

0.0 

J=  1,23 

SUM  ♦  EN( 

)  =  SLM 

1=1.49 

I)  =  1«I 

200 

60 
T  (50H  REC 

AR(  1-16)  = 

700, (LIST 
T  ( /  13X,  12 

UNKNOWN  I 

0  1=1, U9 
0.0 

1  J=1,33 
SUM  ♦  S(  I 
=  SUM 
200. 

900 
T  (  45K   M 
ATU) 

1000,  (LI 
T  (/13X.I2 
1=1,23 
J=34,M9 
,J-33)  =  H 
0300COC.C 
ENCE  CCEFF 
I  =  1,33 
K  =1,33 
0.0 

J  =  1,16 
SUM  +  H12 
I,K)  =  SUM 
=  1C0C0. 
1=1,23 
J=1,33 
)  =  KMIi 
200 
73 
T  (45  H  MA 
C   =((  H11 
150, ((C(I 
200 

150, ( (C(I 
200 
150, ((C(I 


S  =  (G  )(UMTN) 
,49,33,  G,  UNITN,  S  ) 


ESS  COEFFICIENT  MATRIX 
,J),J»1,  11),  I  -1,49) 

J),J=12,22),I  =1,49) 

J),J=23,33),I*1,49) 


(G)UMTN)  //) 


MATRIX,  AR 


I,J)»P(J) 


UNCANT  LCAD  MATRIX. 

A(   )) 
(J),AR(J  ),J=1,16) 
.F2C.8) 
NTERNAL  LCADS.   AT  = 


,J)«P( J) 


ATRIX    OF    ALL    INTERNAL    LOACS. 

A(        )     ) 

ST(I).AT(I    ),    1*1,49) 

,F20.8) 


AR    =    -      F22I 


(S)(P) 


K21 


AT    =    (S)(P)       // 


C    *    (    H11    +    (H12)     (EN))    /    E 


(I,  J) 

ICIENT   MATRIX. 

(I,J)«EN(J,K) 


J)       ♦      H12N(  I,J))/E)«SCALE 

TRIX    OF    CEFLECTION    INFLUENCE    COEFFICIENTS 

+  (H12)(N))/    E1MSCALE)) 
,J),J=    1 ,11), 1=1, 33) 

,J),J=12, 22), 1*1,331 

,J), J=23,33),I=1,33) 


12? 


2 
10 


11 
12 
13 


1U 


15 
20 
30 
31 


32 

33 

35 
36 
3U 
U0 


Ul 

U2 
1*3 


50 


20 
10 


SUB 

CIM 

DO 

CO 

X(  I 

CO 

X(K 

DO 

KP  = 

Z  =  C 

DO 

IF( 

Z  =  A 

KP  = 

CON 

IF{ 

DO 

Z  =  A 

A(L 

A(K 

CO 

Z  =  X 

X(L 

X(K 

IF( 

IF( 

LP1 

CC 

IF( 

RAT 

DO 

A(K 

CO 

X(K 

CCN 

CON 

DO 

11  = 

DO 

S  =  0 

IF( 

IIP 

DO 

S=S 

X(I 

KER 

RET 

KER 

END 

SUB 

DIM 

CC 

CO 

SUM 

CO 

SUM 

C(I 

END 

END 


■  X  »KER ) 
N  A(50,50),X($0,5C) 


ROUTINE    GAUSS3(NtEPiA 

TO 

1  J=l 

J)=0 

K=l 


J 

,K)=1 
3M  L= 
C 

.C 

12  K= 
Z-ABS 
BSF(  A 
K 

TINUE 
L-KP) 
1U  J  = 
(L.J) 
,J)=A 
P,J)  = 
15  J= 
(L,J) 
tJ)  =  X 
P,  J)  = 
ABSF( 
L-N)3 
=  L+1 
36  K  = 
A(K,L 
IO  =  A( 
33  J  = 
,J)=A 
35  J  = 
f  J)=X 
TINUE 
TINUE 
M3  1  = 
N+l-I 
«43  J= 
.0 

II-N) 
1  =  11* 
U2  K  = 
+  A(II 
I,J)  = 

URN 
=  2 


iN 
,N 
.0 
iN 
.0 
lfN 


L,N 

F(AIK.L'iniltl2i  12 
(K,LM 


13, 

L,N 

(KP 

Z 

1,N 

(KP 
Z 

A(L 
1*3 


20,20 
»J) 

,J) 

,L))-EP)50, 50,30 
M,2M 

N 


LP1 

)  )3i 

K.L 

LP1 

(K, 

1,N 

(K, J)-RATIO»X(L,J) 


5i> 
,36,32 
)/A(L,L) 
tN 
J)-RATIO«A(L,J) 


i;n 

1VN 

Ul,ii3,M3 

I  IP  1  ,N 

,K)«X(K.J) 

(X( II,J)-S)/A(I I, II) 


ROUTINE  MATMPY  (L,  M,  N,  A,  B,  C) 
ENSION  A(50,50),  B(5C,50),  C(50,5C) 
10  1=1, L 

in  k=i .m 


E.. 

10  1=1, L 

10  K=1,N 

=  CO 
20   J=1,M 

'  SUM  ♦  A(ItJ) 
,K)  =  SUM 


•  B(  J,K) 


26 


[F]  Matrix 
Input  Form 


1  1 

U0.5060 

2  2 

4712.54 

2  9 

568.242 

220 

656.335    1 

221 

-334.075 

244 

-289.5395 

3  3 

36.7518 

4  4 

35.4129 

5  5 

36.0195 

6  6 

37.5382 

7  7 

41.9716 

8  8 

1554.5632 

814 

-395.2744 

815 

-438.  193 

822 

-109.0759 

823 

-121.1411   1 

9  9 

415.0763 

910 

-123.270 

934 

-335.6280 

944 

-70.3038   1 

9U5 

63.0954 

1010 

271.1617 

1011 

-90.9668 

1034 

236.8452   1 

1035 

-207.059 

1045 

-38.5243 

1046 

42.6268 

11  11 

208.6649   1 

11  12 

-79.0163 

1135 

150.650 

1136 

-162. 1008 

1146 

-27.5895   1 

UU7 

34.1926 

1212 

194.3713 

1213 

-80.2663 

1236 

123.1598   1 

1237 

-156.6048 

1247 

-23.5277 

1248 

31.7581 

1313 

199.0123   1 

131U 

-85.8209 

1337 

124.342 

1338 

-169.2272 

1348 

-23. 198 

13U9 

•  31.517 

1414 

327.3802 

1422 

64.8044 

1438 

137.6857   1 

1449 

-23.8269 

1515 

371.9172 

1516 

-102.9527 

1523 

72.1931   1 

1538 

165.0686 

1543 

-28.6228 

1616 

239.0326 

1617 

-97. 174 

1637 

150.4141 

1638 

-202.5909 

1642 

-28.2054 

1643 

37.8843   1 

1717 

235.8671 

1718 

-96.6455 

1736 

150.466 

1737 

-189.7857   1 

17U1 

-28.9633 

1742 

38.8136 

1818 

257.8147 

1819 

-113.6383   1 

1835 

188.001 

1836 

-198.8093 

1840 

-34.6218 

1841 

42.5174   1 

1919 

330.3452 

1920 

-146.3721 

1934 

280.9922 

1935 

-259.2942   1 

1939 

-45.7466 

1940 

53.9308 

2020 

485.5535 

2021 

-81.2121   1 

203U 

-398.6962 

2039 

75.1214 

2121 

60.9998 

2222 

21.9326 

2323 

24.5350 

2424 

2.4861 

2425 

0.5188 

2525 

1.9818 

2526 

0.4721 

2626 

1.7829 

2627 

0.4193 

2727 

1.6776   1 

2728 

, 0.4194 

2828 

1.6758 

2929 

1.9267 

2930 

0.4860   1 

3030 

2.0048 

3031 

0.5164 

3131 

2.062 

3132 

0.5146   1 

3232 

2.0661 

32  33 

0.5184 

3333 

2.7071 

3434 

1997.6131 

3439 

-133.6859 

3445 

-112.4532 

3535 

1 194.6992 

3540 

-83.7719   1 

3546 

-66.4386 

3636 

822.3501 

3641 

-62.6632 

3647 

-50.6107   1 

3737 

808.2336 

3742 

-57.4131  T 

3748 

-47.0957 

38  38 

804.7401   1 

3843 

-58.4006 

3849 

-48.6683 

3939 

31.4094 

4040 

21.5756 

4141 

16.5636 

42  42 

14.3990 

4343 

13.3832 

4444 

52.8680   1 

4545 

26.3247 

4646 

16.8506 

4747 

13.0860 

4848 

11.6485  1 

4949 

11.0577. 

t 

2 

127 


MATRIX  [p] 
Input  Form 


1 

1500.00   1 

2 

252.80   1 

3 

714.00   1 

4 

528.00   1 

5 

109.00   1 

6 

-210.50   1 

7 

-2782.00   1 

8 

-4523.56 

9 

-5212.17   1 

10 

-5568.74   1 

11 

-5395.86   1 

12 

-4593.58   1 

13 

-4654.94   1 

14 

-4654.94   1 

15 

-4593.58   1 

16 

-5395.86   1 

17 

-5568.74   1 

18 

-5212.17   1 

19 

-4523.56   1 

20 

2390.73   1 

21 

4825.58   1 

22 

5395.02   1 

23 

6125.17   1 

24 

6229.25   1 

25 

5802.52   1 

26 

2738.27   1 

27 

2390.73   1 

28 

4825.58 

29 

5395.02   1 

30 

6125.17   1 

31 

6229.25   1 

32 

5802.52   1 

33 

2738.27   2 

128 


SUBPROGRAM 
CHECK 


Start 


n 


Clear  Arrays 


Form  Product 
CHECK  *   BA  x  AI 


,        I        s 

(   Write  CHECK  J 
"^ 

END 


■(write  A  J 


yesQ 


Print 
is  Singular, 


} 


!29 


PROGRAM   ECHECK 
C  INVERSE   CHECK  FOR  (A) 
C  IS  PRODUCT  OF  (A) (AI) 

DIMENSION  A(50,50) 
C  CLEAR  ARRAYS  FOR  INPUT 

DO  TO  I  =1,50 

DO  10  J=1,5C 

A( I  ,J)  a  O.C 

AI (I,  J)  =  0.0 

CHECK(  I, J)  =  0.0 
10  BA( I, J)  =  0.0 
C  READ  NON-REDUNDANT 
310READ  21«IV JV(A( 


MATRIX.   VARIABLE  CHECK 
AND  SHOULD  GIVE  UNIT  MATRIX 
,AI (50,50 ), CHECK ( 50,50 ),BA( 50  ,50) 


RIX  A,  THEN  PRINT  IT. 

)),  II  ,J1,(A( I1,J1 )),I2,J2,(A(I2,J2)), 
113, J3, (A( I3,J3) )fNEXT 

6), II  ) 


MATI 
ItJ 


21  FORMAT  (4(2I2,F12. 
GC    TC    (31, 32), NEXT 

32  PRINT  33 

33  F0RMAT(42H  (A)  M 
PRINT  150,((A(I,J) 
PRINT  200 

PRINT  15C,( (A(I,J) 
PRINT  200 

PRINT  150,((A(I,J) 
PRINT  200 
DO  40  1=1,33 
DO  40  J=l,32 
40  BA(  I, J)  =  A(I 


tJ) 

AND 
(33,1 


INVERT  MATRIX  ( /> ) 
CALL  GAUSS  3 
GO  TO  (42,41),KER 

41  PRINT  201 

42  CALL  MATMPY(33,33, 
PRINT  150, (  (CHECM 
PRINT  200 

PRINT  150, ((CHECM 
PRINT  200 

PRINT  150, ((CHECM 
150  FORMAT( 11F1C.6) 

200  FORMAT! IHl) 

201  FCRMAT(23H    MATRI 
•  END 

SUBROUTINE  KATMPY 
DIMENSION  A(50,50) 


20 
10 


DO  10  1=1, L 

DO  10  K=1,N 

SUM  =  0.0 

DO  20   J=1,M 

SUM  =  SUM  ♦  A(ItJ)  »  B( J,K) 

C(I,K)  =  SUM 

END 


ATRIX,  ELEVEN  COLUMNS  PER  PAGE////) 
,J=  1  ,111,1  =  1,33) 

,J=12,22),I=1,33) 

,J=23,33),I=1,33) 


FORM  PRODUCT  (A)(AI) 
.E-10,A,AI,KER) 


33,BA,AI,CHECK) 

I, J), J=  1,11), 1  =  1, 33) 

I, J), J=12,22),I=1,33) 

I, J), J*23v33)tl»1«33) 

X  A  IS  SINGULAR) 

(L,  M,  N,  At  B.  C) 
,  B(50,50),  C(50,50) 
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APPENDIX  D 

CALIBRATION  CURVES  FOR  DILLON  DYNAMOMETERS 
AND  HYDRAULIC  PRESSURE  SYSTEM 


•a  ♦«■ 

o 

n  *2 


I 

'  -2-1 

o  -4- 


— I — i.i.    i — i — i — i — 

130000  lb  Riehle  Tensile  Test  Machine 


-O 


O 0-~ — () 


2<?0O  4000  $000  fiboo 

Recorded  Load  —  lb 


fOOOO 


J         4         5  6         7         6^9 

Riehle  Reading  ~  lb  X  10""-> 


■ 


APPENDIX  E 


Note: 


1. 
2. 


Coordinate  Location  of  Strain  Gages 
F8U-3  Wing  Center  Section 

■  0  ■  Intersection  of  the  center  section  droop  leading 

edge  and  the  center  line  of  the  airplane. 

■  0  -  The  center  line  of  the  airplane. 


Upper  Skin  Inside  Gages 


Gage  Number 

Gage  Type 

X 

w 

Y 
w 

1-2 

AX- 5 

63.15 

5.82 

3-4 

it 

83.00 

5.94 

5-6 

•• 

101.95 

5.91 

7-8 

M 

119.25 

5.88 

9-10 

ii 

138.38 

5.94  • 

*L1-12 

•i 

153.78 

5.94 

13-14-15 

•  AR  7-2 

70.66 

33.03 

16-17-18 

•i 

94.83 

33.27 

19-20-21 

ii 

112.88 

33.06 

22-23-24 

•i 

129.34 

33.09 

♦25-26-27 

AR  lor  2 

153.54 

33.15 

28-29-30 

ii 

161.02 

33.33 

31-32-33 

AR  7-2 

127.60 

75.04 

34-35-36 

•i 

142.40 

75.13 

37-38-39 

•i 

155.64 

75.07 

40-41-42 

ii 

168.56 

75.10 

43-44-45 

" 

180.00 

75.25 

46-47-48 

.»» 

153.48 

106.54 

49-50-51 

ti 

161.20 

101.26 

52-53-54 

it 

169.22 

95.79 

55-56-57 

•i 

176.49 

88.92 

58-59-60 

ii 

185.31 

84.92 

61 

A-7 

144.86 

66.10 

62 

ii 

146.09 

65.22 

63 

•i 

150.48 

62.22 

64 

ii 

151.74 

61.44 

65 

ii 

146.67 

68.68 

66-67-68 

AR  7->2 

147.90 

67.87 

69-70-71 

ii 

149.34 

66.94 

72-73-74 

ii 

150.78 

65.92 

75-76-77 

it 

152.22 

64.92 

133 


Upper  Skin  Inside  Gages  (continued) 


Gage  Number 

Gage  Type 

• 

X 

w 

Y 
w 

78 

A- 7 

153.48 

64.05 

79 

•i 

151.20 

72.70 

80 

M 

152.61 

71.74 

81 

II 

154.05 

70.72 

82 

II 

155.52 

69.76 

83-84-85 

AR  7-2 

137.72 

28.53 

86-87-88 

H 

140.00 

32.25 

89-90-91 

II 

142.55 

30.48 

92-93-94 

II 

145.25 

28.74 

95-96-97 

II 

147.99 

45.40 

♦98-99-100 

II 

'   150.30 

43.84 

101-102-103 

II 

152.73 

42.25 

104-105-106 

II 

155.95 

58.68 

107-108-109 

II 

158.11 

57.21 

110-111-112 

II 

160.24 

55.74 

113-114-115 

'  II 

163.84 

71.68 

116-117-118 

•  1 

165.77 

70.51 

119-120-121 

II 

167.72 

69.19 

122-123-124 

II 

169.23 

.   78.83 

125-126-127 

II 

170.60 

78.83 

128-129-130 

II 

173.00 

78.83 

131-132-133 

II 

147.30 

110.81 

Notes 

1.  a)  Gages  1-2,  3-4,  5-6,  7-8,  9-10*  11-12  mounted  by  Chance 
Vought  and  oriented  with  one 
other  leg  perpendicular  to  Y 


Vought  and  oriented  with  one  leg  parallel  to  Y  and  the 


w 
b)  No  gage  factor  or  other  information  available. 

2.  a)  All  AR  7-2  gages  mounted  at  USNPS  and  oriented  with  one 

leg  perpendicular  to  C. I. B. ,  45°  gage  pointing  outboard 
and  toward  trailing  edge,  other  leg  parallel  to  C.  I.  B. 

b)  All  gages  from  lot  B-31,  120.5  +  .5  ohms,  1.97  +  2%, 
b  factor  -  -200. 

3.  a)  All  A-7  gages  mounted  at  USNPS  and  oriented  parallel  to  C.I.  B. 
b)  All  gages  from  lot  B-31,  120.0  +  ,3  ohms,  1.99  4  27.. 

4.  a)  Gages  25-26-27,  28-29-30  mounted  by  Chance- Vought  and 

oriented  approximately  the  same  as  USNPS  mounted  AR  7-2  gages. 

b)  No  gage  factor  or  other  information  available. 
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Bottom  Skin  and  Beam  Assembly 


Gage  Number 


Gage  Type 


w 


w 


134- 
137- 
140- 
143- 
146- 
149- 
152- 
155- 
158- 
161- 
163- 
166- 
169- 
172- 
175- 
178- 
181- 
184- 
187- 
190- 
193- 
196- 
199- 
202- 
204- 
207- 
210- 

*213- 
216- 
219- 
222- 
225- 
228- 

*231- 
234- 
237- 
240- 
243- 
246- 
249- 


135-136 

138-139 

141-142 

144-145 

147-148 

150-151 

153-154 

156-157 

159-160 

162 

164-165 

167-168 

170-171 

173-174 

176-177 

179-180 

182-183 

185-186 

188-189 

191-192 

194-195 

197-198 

200-201 

203 

205-206 

208-209 

211-212 

214-215 

217-218 

220-221 

223-224 

226-227 

229-230 

232-233 

235-236 

238-239 

241-242 

244-245 

247-248 

250-251 


AR  1 


AR  7-2 


AR  1 

AR  1 

AX  5 

AR  1 

AR  1 

AR  7- 

2 

AR  1 

AR  1 

AR  7- 

2 

AR  7- 

2 

AR  1 

AR  7-2 
AR  1 


AX 
AR 


5 
7-2 


AR  1 
AR  1 
AR  7-2 


AR 

7« 

2 

AR 

1 

AR 

1 

AR 

7- 

2 

AR 

1 

AR 

1 

AR 

1 

AR 

7- 

2 

AR 

1 

AR 

7- 

2 

AR 

1 

AR 

1 

72.82 
78.08 
92.18 

104.02 

114.54 

131.50 

142.46 

148.86 

155.01 
70.60 

120.00 
back-up 

146.64 
91.38 

118.29 

134.35 

149.67 

164.74 
back-up 
back-up 
back-up 
back-up 
back-up 
95.13 

127.63 

133.09 
back-up 

141.50 

152.91 

110.60 

133.24 

156.76 

174.23 
back-up 
back-up 
back-up 
back-up 

113.15 

146.31 
back-up 


gage 


gage 
gage 
gage 
gage 
gage 


gage 


gage 
gage 
gage 
gage 


gage 


41.26 
46.75 
61.29 
73.54 
84.50 
102.19 
113.48 
120.06 
126.39 
33.33 
82.28 
to  163-165 
110.90 
40.18 
71.86 
90.78 
. 108.86 
126.64 
to  172-174 
to  175-177 
to  178-180 
to  181-183 
to  184-186 
33.60 
75.43 
82.22 
to  207-209 
92.64 
106.70 
42.61 
72.67 
104.11 
127.33 
to  219-221 
to  222-224 
to  225-227 
to  228-230 
33.18 
82.07 
to  246-248 
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Bottom  Skin  and  Beam  Assembly  (continued) 


Gage  Number 

Gage  Type 

Xw 

Y 
w 

252-253-254 

AR 

7-2 

160.63 

101.50 

255-256-257 

AR 

126-07 

41.11 

258-259-260 

AR 

146.67 

72.07 

♦261-262-263 

AR 

7-2 

164.53 

98.86 

264-265-266 

AR 

174.23 

127.57 

267-268-269 

•  AR 

back-up 

gage 

to 

255-257 

270-271-272 

AR 

back-up 

gage 

to 

258-260 

273-274-275 

AR 

7-2 

back-up 

gage 

to 

261-263 

276-277-278 

AR 

back-up 

gage 

to 

264-266 

279-280-281 

'ar 

7-2 

129.19 

33.78, 

282-283-284 

AR 

7-2 

154.98 

74.65 

285-286-287 

AR 

159.04 

82.16 

288-289-290 

AR 

back-up 

gage 

to 

285-287 

291-292-293 

AR  ' 

7-2 

168.71  . 

96.19 

294-295-296 

AR 

140.63 

39.40 

297-298-299 

AR 

160.60 

72.55 

300-301-302 

AR 

7-2 

172.61 

93.42 

303-304-305 

AR 

192.55 

126.85 

306-307-308 

AR 

back-up 

gage 

to 

294-296 

309-310-311 

AR 

back-up 

gage 

to 

297-299 

312-313-314 

'ar 

7-2 

back-up 

gage 

to 

300-302  . 

315-316-317 

AR 

back-up 

gage 

to 

303-305 

*318-319-320 

.AR 

7-2 

142.55 

34.41 

321-322-323 

AR 

7-2 

147.42 

34.62 

324-325-326 

AR 

7-2 

167.90 

74.74 

327-328-329 

AR 

170.00 

82.16 

♦330-331-332 

AR 

back-up 

gage 

to 

327-329 

333-334-335 

AR 

7-2 

177.00 

90.60 

336-337-338 

AR 

L 

155.64 

41.08 

339-340-341 

AR 

173.60 

73.93 

342-343-344 

AR 

7-2 

180.99 

87.84 

345-346-347 

AR 

201.95 

126.43 

348 

omitted  by  number j 

Lng  error 

349-350-351 

AR 

back-up 

gage 

to 

336-338 

352-353-354 

AR 

back-up 

gage 

to 

339-341 

355-356-357 

AR 

7-2 

,   back-up 

gage 

to 

342-344 

358-359-360 

AR 

back-up 

gage 

to 

345-347 

361-362-363 

it 

153.81 

33.27 

364-365-366 

ii 

157.59 

33.15 

367-368-369 

ti 

159.78 

33.45 
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Bottom  Skin  and  Beam  Assembly   (continued)    . 


Gage  Number 


Gage  Type 


w 


w 


370- 
373- 
376- 
379- 
382- 
385- 
388- 
391- 
394- 
397- 
400- 
403- 
406- 
409- 
412- 
415- 
418- 
421- 

.  424- 
427- 
430- 
433- 
436- 
439- 
442- 

♦445- 
448- 
451- 


371-372 
374-375 
377-378 
380-381 
383-384 
386-387 
389-390 
392-393 
395-396 
398-399 
401-402 
404-405 
407-408 
410-411 
413-414 
416-417 
419-420 
422-423 
425-426 
428-429 
431-432 
434-435 
437-438 
440-441 
443-444 
446-447 
449-450 
452-453 


AR 

1 

AR 

1 

AR 

7-2 

•i 

ti 

•i 

AR 

1 

•  AR 

1 

AR 

7-2 

AR 

1 

•  it 

•• 

it 

AR 

7-2 

AR 

1 

AR  7-2 
ii 


162.58 
160.90 
165.88 
171.08 
176.64 
179.49 
184.92 

back-up 
185.10 
166.57 
170.21 
179.67 
184.65 
190.93 
208.77 
214.08 
141.17 
156.06 

back-up 

back-up 
147.39 

back-up 

97.51 

back-up 

61.05 

back-up 

back-up 

back-up 


g«ge 


gage 
gage 

gage 

gage 

gage 
gage 
gage 


33.06 

30.00 

40.81 

50.93 

69.04 

74.20 

82.10 

to   388-390 

85.13 

36.18 

43.45 

62.64 

72.70 

84.71 

120.12 

130.93 

25.34 

25.09 

to  418-420 

to  421-423 

8.56 
to  430-432 

10.12 
to  436-438 

25.68 
to  442-444 
to   152-154 
to  409-411 


The   following   listed  gages  are  mounted  externally  on 
the  upper  and   lower  skins,  and  are  back-up  gages  as 
indicated,  or  mounted  over  an  internal  beam  member: 


454-455-456 
457-458-459 
460-461-462 
463-464-465 
466-467-468 
469-470-471 
472-473-474 
475-476-477 
478-479-480 
481-482-483 
484-485-486 


AR  7-2 


121.45 

, 

33.61 

back-up 

gage 

to 

279-281 

136.59 

33.61 

151.86 

33.61 

back-up 

gage 

to 

367-369 

165.64 

33.61 

back-up 

gage 

to 

28-30 

back-up 

gage 

to 

89-91 

back-up 

gage 

to 

19-21 

back-up 

gage 

to 

13-15 

back-up 

gage 

to 

161-162 

37 


Bottom  Skin  and  Beam  Assembly  (continued) 

Notes  • 


I.   BEAM  GAGES 


1.   a)   All  AR  1   gages  mounted  on  vertical  webs  of  beams, 

outboard  of  pivot  rib,  were  mounted  by  Chance-Vought 
and  are  oriented  with  center  leg  of  strain  gage 
pointing  inboard  and  on  the  center  line  of  the  web. 
Other  legs  45°  each  side.   (Note:   Gage  178-179-180 
is  an  AR  7-2  mounted  at  USNPS ,  but  oriented  same  as 
above  gages.) 

b)  All  AR  1   gages  mounted  on  vertical  web  of  pivot  rib 
were  mounted  by  Chance-Vought  and  are  oriented  with 
center  leg  of  strain  gage  pointing  aft  and  on  the  center 
line  of  the  web.   Other  legs  45°  each  Side. 

c)  All  AR  1   gages  mounted  on  vertical  webs  of  center 
section  beams  were  mounted  by  Chance-Vought  with  center 
leg  of  strain  gage  pointing  inboard  and  on  the  center 
line  of  the  web.   Other  legs  45  each  side. 

d)  All  AR  7-2   gages  mounted  on  vertical  webs  of  beams, 
outboard  of  pivot  rib,  were  mounted  at  the  USNPS  and  are 
oriented  with  one  leg  perpendicular  to  CIB,  and  other 
legs  pointing  outboard  and  up  or  down  depending  on 
whether  it  is  a  "walk-around"  gage  or  a  "back-up"  gage. 

e)  All  AR  1  gages  mounted  on  vertical  web  Of  intermediate 
Rib  were  mounted  by  Chance-Vought  and  are  oriented 
with  the  center  leg  of  the  strain  gage  pointing  forward 
and  on  the  center  line  of  the  web.   Other  legs  45° 
each  side.' 


II.   LOWER  SKIN 

1.      a)     AX  5     gages  were  mounted  by  Chance-Vought  and  are 

oriented  with  one   leg  parallel   to  the  center   line  of 
the   skin  between   the  beams  and  other  leg  rotated 
90°  away. 
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Bottom  Skin  and  Beam  Assembly  (continued) 


b)  All  AR  1   gages  (361  thru  375)  were  mounted  by  Chance- 
Vought  and  are  oriented  with  one  leg  (the  first  number 
on  each  rosette)  parallel  to  the  center  line  of  the  skin 
between  the  beams;  the  second  leg  (second  number  of 
sequence)  rotated  45°  counter-clockwise;  and  the  third 
leg  (third  number  of  sequence)  rotated  another  45° 
counter-clockwise ,  making  it  perpendicular  to  first 
leg.   The  entire  rosette  points  inboard. 

c)  All  AR  7-2  gages  were  mounted  at  USNPS  and  oriented  with 
one  leg  perpendicular  to  CIB  and  the  entire  rosette 
pointing  outboard  and  forward. 

2.  No  gage  factor  or  other  information  available  on  gages  mounted 
by  Chance-Vought  on  either  the  beams  or  the  lower  skin. 

3.  All  AR  7-2  gages  mounted  at  USNPS  are  from  lot  B-31, 
120.5  +  0.5  ohms,  1.97  +  27.  gage  factor,  b  factor  *  -200. 

it 
The  following  gage  elements  are  inoperable  or  give  questionable 

results: 

11 

25-26-27 

98-99-100 

213 

233 

261 

318-319-320 

332 

445-446-447 
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DIX  F 

STRAIN  :  .  ■ 

The  method  of  taking  strain  -  loyed  common 

'.'heat-stone  bridge  circuitry  shown  in  Fij   F]    The  uniq 
part  of  the  method  \.   1  %e  use  of  a  digital  counter  to 
indicate  bi Ldge  unbalance  and  give  a  visual  display  of  the 
strain  directly  in  units  of  -        ties  per  i.nch„ 

For  the  special  case  considered  3   four 

legs  of  the  bridge  have  the  same  known  resistances  Rs  and 
for  a  small  change  of  resistance  in  one  leg  Ars  it  can  be 
shown  that  the  standard  V/heatstone  bridge  equation  for  out- 
put voltage  V0  reduces  to 

V  R0  AR 


Vo  = 


k   R  (R  +  Rn) 


(F1) 


where  RQ   is  a  constant  resistance  across  terminals  AC  and  V 
is  the  constant  voltage  source „   Further s  the  expression 
for  the  gage  factor  (GF)  that  relates  strain  and  change  of 
resistance  is 


(GP)6 


_  Ar 


(F2) 


R 


By  substituting  equation  (F2)  into  (Fl)  it  is  immediately 
apparent  that  the  output  voltage  is  a  linear  function  of 
strain, 

V0  =  |"  v  (QF)  Ro 
h   (R  +  Ro) 
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Utilizing  this  fact  it  was  then  easy  to  calibrate  the 
electronic  counter  in  units  of  strain0 

Considering  again  the  basic  bridge  circuit,  a  shunt 
resistance  Rc  across  one  of  the  legs  effectively  unbalances 
the  bridge.   Its  effect  on  the  basic  balanced  brid- 
equations  is  developed  below,, 

(GP)e  =  ~ 

R 

Ar  =  R  (GF)-  (pif) 


R 

-  Ar  = 

1 

R    Br 

1         1 

R         Rr- 

R  +   Rc 

A 

R  =  R  - 

R  Rc 

_        R2 

(P5) 


R  +  RC     R  +  Rc 
Substituting  equation  (F5)  into  (F4 )  and  solving  for  Rc 
gives, 


R   =  R    _J_  -  1 

Le(GF) 


(F6) 


From  equation  (F6)  the  shunt  resistance  necessary  to 
produce  an  equivalent  strain  of  0.001  in*  per  in0  when  R 
is  120  ohms  and  (GF)  equals  2.0  is  then, 


I  If  2 


Rc  =  120 


.■ 


'59780  ohms 


_(0o001)(2c0) 

If  on  the  other  hand  a  more  standard  60000  ohm  resistor 
were  used  for  R^,  the  equivalent  strain  would  be  993 
micro  inches  per  inch„   This  feature  was  used  to  calibrate 
the  electronic  counter,, 

A  Wheatstone  bridge  circuit  was  permanently  constructed 
into  which  precision  120  ohm  (*  ,25r5)  resistors  v/ere 
readily  connected  to  form  the  four  legs.  A  precision 
60000  ohm  {t   .057$)  shunt  resistor  was  connected  with  a 
switch  as  shown  in  Fig,  Fl„   In  addition  a  precision  50000 
ohm  potentiometer  was  provided  through  connection  CE  as  a 
variable  resistance  for  any  delicate  balancin 

With  the  switching  and  balancing  unit  and  the  switch  F 
disconnected,  the  amplifier  connected  across  terminals  AC 
and  the  variable  resistance  connected  across  CE,  we  have 
the  special  bridge  circuit  described  previously„   The 
variable  resistance  was  then  adjusted  to  exactly  balance 
the  bridge  and  consequently  give  a  zero  indication  on  the 
electronic  counter <,   Switch  F  uo.s   then  closed  and  the 
amplifier  gain  adjusted  to  give  a  993  reading  on  the  counter,, 
The  counter  had  been  previously  set  to  repeatably  indicate 
the  Integrative  result  of  a  one  second  sampling  of  V0  which 
greatly  enhanced  smooth  readings « 
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Repeating  this  procedure  and  noting  the  consistancy 
of  returning  to  zero  with  switch  7   open  and  to  99$  with 
the  switch  closed  insured  the  establishment  of  a  linear 
calibration  of  the  counter  to  read  in  units  of  micro-inches 
per  inch* 

After  calibrating  the  counter s  switches  E  and  F  were 
opened  and  resistances  R,  and  R~  were  replaced  by  the 
switching  and  balancing  unit  leads  for  the  active  and 
compensating  gages  respectively  as  shown  in  Fig0  Fl.   The 
active  and  compensating  gages  also  had  120  ohm  resistances c 
Each  of  the  20  gages  connected  through  the  switching  and 
balancing  unit  was  then  zeroed  on  the  switching  and  balancing 
unit  in  the  conventional  manner  prior  to  subjecting  them  to 
load.   At  any  time  during  the  test  run  switch  F  could  be 
closed  and  a  change  of  99$  observed  on  the  counter.   This 
p-ave  a  continual  check  on  the  calibration  of  the  counter „ 
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Fig,  F1 

SCHEMATIC  DIAGRAM    . 
STRAIN  GAGE  INSTRUMENTATION 


20  Channel 
Switch.  & 
Bal.  UnJt 
Baldwin  - 
Lima-Hamilton 


DC  Pwr  Supply- 
Model  721A  - 
Hewlett- 
Packard- 
O — Q 


O 


-*o 


Amplif j er 
Model  111BF 

Kintell 


I 


Volt.-to-f req. 

Converter 
Model  DY-2210 

Dymec,  Inc. 


Electronic  Counter 
Model  521 DL 
Hewlett  Packard 


R1  =  R2  =  R3  =  Rlf  =  CRCA  120  ohm  {.25%)    21  5-RL 

Rc  =  CRCA  60000  ohm  (.05%)    21 5-RL 

RT  =  50000  ohm  Helipot  TP  Precision  Potentiometer 
Resistance  Tolerance:  5% 
Linearity  Tolerance:  .5% 
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APPENDIX  G 


EXPERIMENTAL  PRINCIPAL  STRESSES  AND 
AXIS  ORIENTATION 

The  magnitudes  of  principal  stresses  and  the  orientation 
of  the  principal  axes  were  calculated  for  all  rosette 
locations  of  interest.   Experimental  strain  readings  are 
listed  in  Table  Gl.   Calculations  were  performed  on  the 
C.D.C.  1604  Digital  Computer  utilizing  a  FORTRAN  program 
named  ROSRED  shown  in  Tables  G2  and  G3o   The  computer  output 
is  listed  ih  Table  G40 

The  equations  solved  by  the  program  were  the  relations 
between  rectangular  rosette  readings  and  principal  stresses 
found  in  any  standard  text  on  the  subject 


(j. 


max 


<j. 


mm       — 


F,+G 


„d-/i)      (1+//) 

€i+€3    -  _J_ 
(V/i)       Of  (!) 


(F  -F  ) 


2tV<£,  +^3) 


<e,  -f3)2  +[2e2-(e/+eJ) 


r, 


max 


2(1+//) 


V  <e,-e5)2  +[2F2-(c„  +e3)] 


0     =  -§-  arc  tan 


262°  (e,  +e3^ 
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0p  is  the  angle  from  the  axis  of  €]_  to  the  maximum  normal 
stress  axis  and  6j_,  €3  and  ^3  are  strain  readings  in 
the  respectively  numbered  gages  shown  in  Fig.  Gl.   Rosettes 
used  in  this  experiment  were  numbered  with  three  con- 
secutive numbers  and  the  lowest  number  identified  the 
rosette.   Logically  then  (JlL  is  the  angle  from  the  axis  of 
the  lowest  numbered  gage  in  each  rosette  to  the  maximum 
normal  stress  axis,   A  positive  value  indicates  an  angle 
in  the  direction  of  €2* 

The  input  to  program  ROSHED  (Table  Gl)  was  designed 
to  accommodate  this  notation.   Only  the  rosette  identifying 
number  is  listed  which  is  also  the  gage  number  for  which 
the  first  column  of  strains  are  listed.   The  second  and 
third  columns  are  then  the  values  of  strain  on  the  next 
two  consecutive  numbered  gages  comprizing  the  rosette. 

The  results  listed  in  Table  G4  are  identified  by  the 
rosette  numbers  which  quickly  orients  the  axis  reference 
for0p. 

It  should  be  pointed  out  that  the  program  was  designed 
for  specific  values  of  E  and  jx  ,  namely  10.3  x  10°  psi 
and  0.32  respectively. 
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Fig.  G2 
FLOW  CHART  FOR  PROGRAM  ROSRED 


Return 


Start 


READ  Ros 

e 


osette  Number,   ^ 
1i>  €2j,  ^31.  J 


Ai  =€U  +63. 


Bi  =  eU  -  €3i 


Ri  =   /Bi2   +    (2  €2jL   -  Ai)2' 

♦ 

cr     -  i&£ 

"  max         2 

'  Ai  .   +     1        Ri 

L0.68        1,32      j 

•  r                           • 

cr-    =  10-3 

wmm        .p 

'  Ai_  -     1        Ri 

L0.68       1.32     J 

i ' 

%ax  =  l^% 

4>P 

=  £  arc tan 

[2  e2i  -  Ai] 

L          Bi          J 

no 


yes 


Q 


PRINT  Rosette  Number, 
^max'    ^min'    Tmax?  <Pp 


STOP 
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Table   G1 
EXPERIMENTAL  STRAIN  READINGS 
(Input  To  RO  SITED  Progranj) 


8000  It 

Test 

Rosette 

1  i 

C 

e, 

Rosette)  f 

€0 

^ 

Number 

' 

2 

3 

Number 

» 

*■  2 

13   -76. 

60. 

34. 

193 

131. 

8.  ' 

-130. 

16   -12. 

83. 

30. 

196 

18. 

-1  10. 

-20. 

19 

-  5. 

76. 

27. 

199 

140. 

32. 

-124. 

22 

-  1. 

68. 

6. 

204 

33. 

135. 

0. 

28 

-22. 

136. 

25. 

207 

-19. 

0. 

34. 

31 

-10. 

124. 

42. 

210 

46. 

-13. 

-53. 

34 

-12. 

124. 

9. 

216 

32. 

134. 

-26. 

37 

1. 

136. 

1. 

219 

-141. 

-11. 

130.  , 

40 

23. 

143. 

-17. 

222 

-179. 

-10. 

191. 

•  43 

-22. 

53. 

4. 

225 

5. 

68. 

-5. 

46 

-29. 

134. 

31. 

228 

-117. 

15. 

111. 

49 

-23. 

166. 

29. 

234 

166. 

46. 

-188. 

52 

16. 

179. 

18. 

237 

1. 

75. 

2. 

55 

16. 

200. 

-16. 

240 

90. 

2. 

-101. 

58 

42. 

158. 

-21. 

243 

24. 

94. 

-10. 

66 

5. 

98. 

3. 

246 

-17. 

1. 

26. 

69 

20. 

141. 

-18. 

249 

40. 

-10. 

-46. 

72 

-26. 

132. 

22. 

252 

23. 

160. 

-18. 

75 

-9. 

93. 

6. 

255 

-1. 

-104. 

101. 

83 

5. 

51. 

-  7. 

258 

-97. 

-19. 

-6. 

86 

7. 

60. 

-11. 

264 

-51. 

5. 

49. 

89 

8. 

53. 

-17. 

267 

110. 

20. 

-99. 

92 

7. 

46. 

-21. 

270 

9U. 

20. 

-74. 

95 

5. 

83. 

-4. 

273 

3. 

-28. 

-1. 

101 

14. 

83. 

—  6. 

276 

29. 

0. 

-61. 

104 

11. 

106. 

-6. 

279 

14. 

93. 

3. 

107 

10. 

105. 

-11. 

282 

11. 

162. 

14. 

110 

16. 

103. 

-13. 

285 

2. 

1. 

3. 

113 

9. 

132. 

-6. 

288 

0. 

-6. 

-5. 

116 

15. 

130. 

-17. 

291 

3. 

175. 

7. 

119 

18. 

130. 

-14. 

294 

-48. 

2. 

50. 

122 

18. 

148. 

-15. 

297 

-40. 

2. 

35. 

125 

12. 

146. 

-17. 

300 

1. 

-15. 

-6. 

128 

15. 

144. 

-23. 

.  303 

73. 

-8. 

-71. 

131 

-52. 

147. 

44. 

306 

-36. 

11. 

-45. 

134 

164. 

-14. 

-194. 

309 

34." 

6. 

-34. 

137 

200. 

-32. 

-192. 

4. 

-20. 

2. 

140 

213. 

0. 

-207.         318 

-7. 

73. 

13. 

,143 

215. 
135. 

-11. 
-2. 

-219.         333 

-23. 
-2. 

173. 
91. 

23. 
-12. 

rl46 

-134. 

<  342 

149 

163. 

-13. 

-162. 

355 

-9. 

106. 

31. 

152 

-10. 

176. 

25. 

367 

14. 

-U9. 

24. 

155 

189. 

-4. 

-182. 

394 

-30. 

170. 

49. 

158 

214. 

75. 

-176. 

1 

409 

-22. 

155. 

5. 

163 

12. 

2. 

-3. 

454 

-15. 

89. 

31. 

166 

15. 

-9. 

-20. 

457 

-14. 

95. 

22. 

169 

66. 

162. 

-U1.         466 

11. 

~95- 

-46. 

172 

-221. 

-47. 

153.       ,  «*69 

18. 

42. 

-21. 

175 

-241. 

-4. 

240. 

472    -11. 

83. 

35. 

178 

-134. 

3. 

135. 

475    -23. 

59. 

14. 

181 

-7. 

-141. 

-2. 

1 

478 

4. 

64. 

-16. 

184 

-164. 

-11. 

158. 

481 

39. 

59. 

-66. 

187 

159. 

33. 

-120. 

'• 

484 

53. 

103. 

-44. 

.190 

223. 

1. 

-132. 

(Table   G1    Continued) 


7000  lb 

Test 

Rosette 
Number 

1    e, 

^2 

£3 

JTosette 

Number 

C| 

^2 

13 

i 
-64. 

53. 

28. 

237 

+   2. 

+  69. 

+  2 . 

16 

-11. 

72. 

28. 

243 

19. 

183. 

-  1 0. 

19 

-7. 

64. 

22. 

246 

-  19. 

0. 

22 

22 

-2. 

58. 

2. 

249 

35. 

-  10. 

-  4  3. 

28 

-20. 

18. 

25. 

252 

18. 

138. 

-  18- 

31 

-8. 

106. 

35. 

255 

-   3. 

-100. 

1  \.J   • 

90. 

34 

-13. 

107. 

6. 

258 

-  85. 

-  15. 

70. 

37 

1  . 

1  19. 

1. 

264 

-  45. 

3. 

44; 

40 

18. 

125. 

-15. 

267 

98. 

'  18. 

-  87. 

43 

-21. 

46. 

2. 

270 

80. 

•18. 

-  65. 

46 

-29. 

1  17. 

26. 

273 

2. 

-  25. 

-   2. 

49 

-21. 

143. 

24. 

276 

25. 

-   2. 

-  50. 

52 

12. 

157. 

14. 

279 

12. 

86. 

5. 

55 

14. 

174. 

-15. 

282 

13. 

144. 

14. 

58 

35. 

138. 

-20. 

285 

1  . 

2. 

3. 

66 

-  2. 

+  84. 

-   4. 

288 

1. 

-   4. 

-   4. 

69 

+  15. 

+  126. 

.  -  17. 

291 

4. 

156. 

8. 

72 

-27. 

+  1  17. 

+  16. 

294 

-  40. 

3. 

43. 

75 

-12. 

+  84. 

+   4. 

297 

64. 

2. 

29. 

83 

+  6. 

+  48. 

-   5. 

300 

1. 

-  17. 

-   7. 

86 

+  7. 

+  56. 

-   4. 

303 

70. 

-   7. 

-  66. 

89 

+  8. 

+  50. 

-  11. 

306 

41. 

8. 

-  40. 

92 

+  8. 

+  45. 

-  15. 

309 

28. 

3. 

29. 

95 

+  5. 

+76. 

-2. 

312 

3. 

-  20. 

.  4. 

98 

+  12. 

+  77. 

-   3. 

315 

-  76. 

-  20. 

70. 

101 

+  14. 

+  77. 

■  -   4. 

321 

-   6. 

71. 

19. 

104 

+  8. 

+  95. 

-   4. 

324 

-   8. 

144.  ' 

32. 

107 

+  8. 

+  94. 

-   7. 

327 

10. 

3. 

-  10. 

no 

+  14. 

+  94. 

-  10.^ 

333 

20. 

151. 

22. 

113 

+  7. 

+  1  16. 

-   3. 

336 

46. 

4. 

-  43. 

116 

+   4. 

+  105. 

-  21. 

339 

68. 

10. 

-  66. 

119 

+   6. 

+  105. 

-19.  , 

342 

-   4. 

85. 

-   7. 

122 

+  22. 

+  120. 

-  21. 

345 

115. 

2. 

-104. 

125 

0. 

+  121  . 

-  23. 

349 

-  47. 

-.  1. 

44. 

128 

+   2. 

+  1  16. 

-  23. 

352 

-  74. 

1. 

75. 

131 

-  52. 

+  123. 

+  31. 

355 

-   9. 

96. 

27. 

134 

+  151. 

+  91  . 

-166. 

358 

110. 

-  18. 

-128. 

137 

+  181. 

-  28. 

-166. 

361 

-  13. 

-  60. 

15. 

140 

+  194. 

0. 

-181  . 

364 

-   6. 

-  44. 

23. 

143 

+  198. 

-   7. 

-187. 

367 

13. 

-  42. 

21. 

146 

+  128. 

0. 

-1  17. 

370 

-   2. 

-  58. 

18. 

149 

+  153. 

-11.. 

-141. 

373 

9. 

-  20. 

21. 

152 

-   7. 

+  157. 

+  25. 

376 

S-  41. 

50. 

7. 

155 

+  169. 

-   5. 

-161. 

379 

-  31. 

77. 

26. 

158 

+  191. 

69. 

-156. 

382 

i-  20. 

105. 

21. 

163 

+  10. 

0. 

-   6. 

385 

I-  21. 

124. 

36. 

166 

+  16. 

-   7. 

-  20. 

388 

-  13. 

1. 

1  1. 

169 

+  55. 

+  141. 

-  36. 

391 

43. 

-   4. 

-  46. 

172 

-187. 

-  41  . 

+  137. 

394 

-  26. 

150. 

40. 

175 

-216. 

-   3. 

+  213. 

397 

182. 

14. 

-190. 

178 

-121. 

+   7. 

+  121  . 

400 

165. 

0. 

-167. 

181 

-   8. 

-128. 

-   2. 

403 

}    216. 

30. 

-209. 

184 

-147. 

-  11. 

+  143. 

406 

224. 

0. 

-230. 

187 

+  143. 

+  24. 

-101. 

409 

-  20. 

134. 

7. 

190 

+  197. 

0. 

-203. 

412 

229. 

83. 

-189. 

193 

+  121. 

+   8. 

-1  18. 

415 

240. 

15. 

222. 

196 

+  16. 

-100. 

-  16. 

418 

-   9. 

-   3. 

6. 

199 

+  126. 

+  30. 

-112. 

421 

-  10. 

-   5. 

4. 

204 

+  32. 

+  121. 

-   3. 

424 

7. 

-   1. 

-   5. 

207 

-  15. 

+   1. 

+  33. 

427 

0. 

-   1. 

-   2. 

210 

+  43. 

-  11. 

-  46. 

430 

-  18. 

-  10. 

1  1. 

213 

+  26. 

+  109. 

-  14. 

433 

-  15. 

-   9. 

9. 

216 

+  30. 

+  121. 

-  19. 

436 

-   4. 

5. 

5. 

219 

-124. 

-   8. 

+  1  14. 

439 

-   4. 

10. 

15. 

222 

-157. 

-   8. 

+  165. 

442 

77. 

31. 

-  73. 

225 

+   3. 

+  69. 

-   4. 

448 

10. 

159. 

0. 

-108. 

+  13. 

+  98. 

j   451 

7. 

162. 

3. 

J 

234 

+  143. 

+  39. 

-165.. 

' 

1 

1 

Computer 
Coded 

Name 

EP1 

EP2 

EP3 

SIGMAX 

SIGMIN 

TAUMAX 

PHIP 

GAGE 

N 


TABLE  G2 

LIST  OF  SYMBOLS  USED  IN 
ROSRED  PROGR 

Definition 

€^»  strain  in  lowest,  numbered  gage  of  rosette 

C2>  strain  in  diagonal  gage 

Co,  strain  in  perpendicular  gage 
Cj^^  maximum  principal  stress 
Otnmi   minimum  principal  stress 
Tmaxi   maximum,  shearing  stress 
0  ,  angle  from  the  Ci  axis  to  C*max  axis « 
Rosette  number  (lov/est  numbered  gage  in  rosette) 
Number  of  rosettes  furnished  as  data  input 
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Table  G3 

FORTRAN  PROGRAM 
"ROSRED" 


PROGRAM  ROSRED 
C      PROGRAM  TO  OBTAIN  PRINCIPAL  STRESSES  AND  AXIS  ORIENTATION  FROM 
C      RECTANGULAR  ROSETTE  STRAIN  DATA 

OD  I  MENS  I  ON  EP1  ( »+00  )  ,EP2  ( i*00  )  ,  EP3  (  400  ) ,  SIGMAX(  400  J ,  SIGMINUOO  )  , 
1TAUMAXU00)  ,  PHIPUOO)  , GAGE (4 00) 
READ  1 ,N 
1  F0RMATU3) 

READ  2,  (GAGE( I ),EP1( I),  EP2( I),EP3(  I),    I=1fN) 
C      FOR  EACH  ROSETTE,  LOWEST  NUMBER  GAGE  IDENTIFIES  ROSETTE  AND 
C      GIVES  REFERENCE  POINT  FDR  ANGLE  PHIPR INCI PAL .WHERE  POSITIVE 
C      ANGLE  IS  MEASURED  TOWARDS  GAGE  2  OF  ROSETTE 
C      STRAIN  DATA  IS  MICRO-I HZ HES/ INCH 
-   2  FORMAT  (I3.3F8.0) 
DO  20  I  =  l.N 
A( I)=  EP1 ( I )  +  EP3( I) 
B(I )=  EP1 ( I)-  EP3( I  ) 

R(I)=  SQRTF(B(I)«*2  +  (EP2(I)*2.0  -A(I))««2) 
SIGMAX(I)  =  (10.3/2.0) »{ (A( I )/. 68)  +  (  1  .0/ 1 . 32) »R(  I )  ) 
SIGMIN(I)=  5.15  ♦( (A( I)/.68)-(1  .0/ 1 . 32) *R( I )  ) 
TAUMAX(I)=  ( 10.3/2.64)*R(I) 
20  PHIP(I)=  ( (ATANF( (2.0»EP2( I )  -A( I ) ) /B ( I ) ) ) /2 .0)  *  57.3 
PRINT  30 

30  FORMAT  (43H1  MAX        MIN        MAX         P) 
PRINT  31(GAGE( D.SIGMAXt I),SIGMIN(I ),TAUMAX( I), PHIP(I)  I  *1,N) 

31  FORMAT  (1X,I3,3F10.0, 1F10.2) 
STOP 

END 
END 


152 


Table  G*f 

EXPERIMENTAL  RESULTS 

(Output  of  "ROSRED"  Program) 


800Q  lb 

Test 

Rosette 
Number 

Cf  >Psi 
max 

(J    .  >Psi 

min 

max 

0^Deg, 

13 

446. 

-1082. 

764. 

-27.91 

16 

737. 

-464. 

600. 

-37.08 

19 

689. 

-356. 

522. 

-38.09 

22 

550. 

-474. 

512. 

-43.47 

28 

1088. 

-1043. 

1065. 

-40. 05- 

31 

1109. 

-624. 

867. 

-38.23 

3U 

960. 

-1005. 

983. 

-42.61 

37 

1069. 

-1038. 

1  053. 

.00 

40 

1  149. 

-1058. 

1  104. 

40.94 

43 

358. 

-631. 

494. 

-39.08 

46 

1079. 

-1049. 

1  064. 

-38.65 

49 

1333. 

-1243. 

1288. 

-40.47 

52 

1522. 

-1007. 

1264. 

-44.83 

55 

1566. 

-1566. 

1  566. 

42.72 

58 

1336. 

-1018. 

1  177. 

38.98 

66 

794. 

-673. 

734. 

44.70 

69 

1  118. 

-1087. 

1  102. 

41.  14 

72 

1032. 

-1093. 

1  062. 

-39.93 

75 

717. 

-762. 

740. 

-42.73 

83 

393. 

-424. 

408. 

41.71 

86 

459. 

-519. 

489. 

40.87 

89 

391  . 

-527. 

459. 

38.87 

92 

322. 

-534. 

428. 

37.60 

95 

652. 

-637. 

645. 

43.44 

101 

682. 

-561. 

621. 

41.40 

104 

848. 

-772. 

810. 

42.66 

107 

820. 

•       -835. 

827. 

42.16 

110 

823. 

-777. 

800. 

40.94 

113 

10U3. 

-997. 

1  020. 

43.36 

1  16 

1015. 

-1045. 

1  030. 

41.52 

119 

1037. 

-976. 

1007. 

41  .44 

122 

1  173. 

-1128. 

1  150. 

41.79 

125 

1  126. 

-1202. 

1  164. 

42.21 

128 

1  104. 

-1225. 

1  164. 

41.35 

131 

1  176. 

-1297. 

1236. 

-36.19 

134 

1  170. 

-1624. 

1  397. 

.16 

137 

1616. 

-1494. 

1  555. 

-5.20 

140 

1684. 

-1593. 

1  639. 

-.41 

143 

1664. 

-1725. 

1695. 

-1.19 

146 

1057. 

-1042. 

1050. 

-.53 

149 

1280. 

-1265. 

1272. 

-2.37 

152 

1435. 

-1208. 

1  322. 

-42.04 

155 

1502. 

-1396. 

1449. 

-1.16 

158 

1871. 

-1295. 

1583. 

8.01 

163 

130. 

6. 

62. 

-9.22 

166 

108. 

-184. 

146. 

-10.19 

169 

1428. 

-1050.. 

1239. 

35.16 

172 

948. 

-1978. 

1463. 

1.99 

175 

1869. 

-1884. 

1877. 

.42 

178 

1057. 

-1042. 

1  050. 

-.53 

181 

997. 

-1133. 

1  065. 

44.48 

184 

1212. 

-1303. 

1258. 

1.42 

187 

1389. 

-798. 

1094. 

2.76 

190 

2117. 

-739. 

1428. 

-7.04 

153 


(Table  Gh   Continued) 


8000  lb 

Test 

RosS'ett 
Number 

e  (j        ,t>si 

^max ' . 

Cj    •    ,^si 

max7  - 

0p,  Deg. 

193 

1028. 

-1012. 

1020. 

1  .64 

196 

848. 

-879. 

863. 

-40.06 

199 

1  168. 

-926. 

1047. 

5.15 

204 

1  184. 

-684. 

934. 

41.04 

207 

329. 

-101. 

215. 

7.90 

^210 

340. 

-446. 

393. 

-5.43 

216 

1092. 

-1002. 

1047. 

38.76 

219 

975. 

-1141. 

1  058. 

1.16 

222 

15U0. 

-1358. 

1449. 

2.47 

225 

532. 

-532. 

532. 

42.90 

228 

855. 

-946. 

901  . 

-4.49 

234 

1284. 

-1618. 

1451. 

8.93 

237 

596. 

-551. 

574. 

-44.81 • 

240 

664. 

-831  . 

7U7. 

2.25 

243 

798. 

-586. 

692. 

39.47 

2U6 

238. 

-102. 

170. 

4.62 

"  249 

295. 

-385. 

340. 

-4.62 

252 
255 

1277. 

-12C1  . 

12  39. 
1  266. 

41.30 
35.84 

2023. 

-508. 

258 

-344. 

-1216. 

436. 

-17.77 

264 

378. 

-408. 

393. 

-3.42 

267 

907. 

-7U0. 

823. 

3.95 

270 

812. 

-509. 

660. 

3.39 

273 

242. 

-212. 

227. 

-43.03 

276 

130. 

-615. 

373. 

9.79 

279 

.  790. 

-532. 

661. 

43.14 

282 

1356. 

-977. 

1  167. 

-44.72 

285 

50. 

26. 

12. 

35.79 

288 

-4. 

-71  . 

34. 

-27.23 

291 

1402. 

-1251. 

1327. 

-44.67 

294 

398. 

-367.. 

382. 

-.58 

297 

257. 

-333. 

295. 

-3.42 

300 

63. 

-139. 

101. 

-37.  18 

303 

581. 

-55  1 . 

566. 

-3.56 

306 

-210. 

•  -1017. 

403. 

42.51 

309 

269. 

-269. 

269. 

5.00 

312 

225. 

-134. 

180. 

-43.76 

318 

597. 

-506. 

552. 

-40.94 

333 

1362. 

-1362. 

1  362. 

-41.22 

3U2 

660. 

-872. 

766. 

43.54 

355 

924. 

-591  . 

758. 

-39.06 

367 

820. 

-2U4. 

532. 

42.90 

394 

1434. 

-1  146. 

1  290. 

-38.09 

409 

1151. 

-1409. 

1280. 

-42.64 

454 

778. 

-536. 

657. 

-37.08 

457 

784. 

-663. 

724. 

-39.41 

466 

379. 

-909. 

644. 

-34.91 

469 

349. 

-395. 

372. 

32.93 

472 

764. 

-40  1. 

582. 

-36.03 

475 

448. 

-584. 

516. 

-36.88 

478 

461. 

-643. 

552. 

40.94 

481 

494. 

-903. 

698. 

27.05 

484 

925. 

-789. 

857. 

31.89 
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(Table  G*f  Continued) 


7000  lb 

Test 

Rossette  ^    psi 

Number   vmax" 

(J    .  >Psi 
wmin 

T   ?pa 

max 

0  ,  Deg. 

13 

387. 

-933. 

660. 

-28.53 

16 

647. 

-390. 

518. 

-36.47 

569. 

-342. 

455. 

-37.81 

22 

453. 
251. 

-453. 
-175. 

453. 

-44.02 

28 

213. 

-17.28 

31 

946. 

-537. 

741. 

-38.46 

3U 

812. 

-918. 

865. 

-42.55 

37 

936. 

-906. 

921  . 

.00 

40 

995. 

-950. 

972. 

41  .20 

43 

298. 

-566. 

442. 

-39. 15 

46 

927. 

-972. 

949. 

-38.47 

1+9 

1  141  . 

-1095. 

1  118. 

-40.49 

52 

1321. 

-927. 

1  124. 

-44.80 

55 

1359. 

-1374. 

1  366. 

42.63 

_.  58 

1154. 

-927. 

1041  . 

39.05 

66 

633. 

-724. 

679. 

44.67 

69 

984. 

-1014. 

999. 

41.41 

72 

887. 

-1054. 

970. 

-40.03 

75 

629. 

-750. 

689. 

-42.41 

83 

381. 

-366. 

373. 

41.70 

86 

450. 

-405. 

427. 

42.12 

89 

386. 

-431. 

409. 

39.78 

92 

336. 

-442. 

389. 

38.33 

95 

605. 

-559. 

582. 

43.66 

98 

637. 

-501. 

569. 

42.05 

101 

642. 

-490. 

566. 

41  .44 

104 

757. 

-697. 

727. 

43.16 

107 

740. 

-724. 

732. 

42.71 

no 

754. 

-694. 

724. 

41.29 

113 

921  . 

-860. 

890. 

43.75 

116 

762. 

-1020. 

891. 

41.86 

119 

777. 

-974. 

875. 

41.80 

122 

955. 

-940. 

947. 

39.90 

125 

864. 

-1212. 

1  038. 

42.52 

128 

833. 

-1151. 

992. 

42.  18 

131 

932. 

-1250. 

1091. 

-36.37 

134 

1343. 

-1570. 

1  456. 

15.93 

137 

11495. 

-1268., 

1  382. 

-5.78 

140 

1562. 

-1365. 

1  464. 

-.99 

143 

1589. 

-1422. 

1 505.  ' 

-1.86 

146 

1040. 

-874. 

957. 

-1.29. 

149 

12U6. 

-1064. 

1  155. 

-3.30 

152 

1298. 

-1025. 

1162. 

-41.92 

155 

1350. 

-1229. 

1289. 

-1.56 

158 

1677. 

-1  147. 

1412. 

8.27 

163 

95. 

-34. 

64. 

-7.02 

166 

115. 

-176. 

146. 

-7.76 

169 

1230. 

-942. 

1086. 

35.46 

172 

892. 

-1649. 

1270. 

2.82 

175 

1651  . 

-1697. 

1  674. 

.20 

178 

946. 

-946. 

946. 

-1.66 

181 

884. 

-1036. 

960. 

44.30 

184 

1  103. 

-1164. 

1  134. 

1.78 

187 

1270. 

-634. 

952. 

.70 

190 

1515. 

-1606. 

1561. 

.43 

193 

957. 

-91  1. 

934. 

1.56 

196 

790. 

-790. 

790. 

-40.46 

199 

1052. 

-840. 

946. 

5.47 

204 

1062. 

-623. 

842. 

40.34 

207 

334. 

-6  1. 

197. 

9.22 

i  21-0 

332. 

-378. 

355. 

-6.03 

213 

910. 

-728. 

819. 

39.51 

216 

1005. 

-838. 

.  921  . 

39.01 

219 

853. 

-1005. 

92  9. 

.72 

222 

1320. 

-1 199. 

1260. 

2.13 

225 

535. 

-551. 

543. 

43.56 

228 

740. 

-892. 

816. 

-4.96 

234 

1097. 

-1430. 

1263. 

8.99 

237 

553. 

-493. 

523. 

.00 

155 


(Table  Gk   Continued) 


7000  lb 

Test 

Rosette 

Nnrnhpr 

(J        »?si 

CIA  x  * 

^min'?si 

Wsi 

0,  >  D*s 

243 

1466. 

-1329. 

1  397., 

42.68 

246 

183.  ' 

-138. 

160. 

2.09 

249 

247. 

-368. 

308. 

-4.37 

252 

1086. 

-1006. 

1  086. 

41.29 

255 

1836. 

-518. 

1  1-77. 

36.03 

25.8 

494. 

-72  1  . 

608. 

2.76 

261+ 
2W 

341. 

-356. 

348. 

-2.25 

812. 

-645. 

728. 

3.85 

270 

685. 

.  -458. 

572. 

4.  12 

273 

196. 

-196. 

196. 

-42.72 

276 

115. 

-493. 

304. 

7.82 

279 

734. 

-477. 

605. 

43.71 

282 

1223. 

-314. 

1018. 

-44.89 

285 

38. 

22. 

8. 

.00 

288 

5. 

-50. 

28. 

-22.50 

291 

1261. 

-1080. 

1  171. 

-44.62 

294 

347. 

-301. 

324. 

-1 .04. 

297 

1077. 

331  . 

373. 

-34.27 

300 

68. 

-159. 

114. 

-37.03 

303 

566. 

-505. 

535. 

-3.77 

306 

329. 

-314. 

321. 

5.25 

309 

631  . 

233. 

199. 

44.44 

312 

236. 

-130. 

183. 

44.39 

315 

539. 

-630. 

585. 

6.56 

321 

61  1  . 

-414. 

513. 

-39.52 

324 

1224. 

-860. 

1042. 

-40.70 

327 

81. 

-81. 

81. 

8.35 

333 

1333. 

-696. 

1  014. 

-44.78 

336 

371  . 

-325. 

348. 

1.61 

339 

543. 

-512. 

527. 

3.83 

342 

623. 

-790. 

706. 

44.53 

345 

938. 

-772. 

855. 

-.92 

349 

332. 

-378. 

355. 

-.31 

352 

•   589. 

-574. 

581. 

-.19 

355 

830. 

-557. 

693. 

-39.16 

358 

795. 

-1068. 

931. 

-2.16 

361 

•504. 

-473. 

488. 

38.54 

364 

554. 

-296. 

425. 

37.28 

367 

719. 

-204. 

461. 

43.06 

370 

642. 

-400. 

521. 

40.70 

373 

504. 

-50. 

277. 

40.14 

376 

298. 

-813. 

555. 

-35.15 

379 

621. 

-697. 

659. 

-35.14 

382 

839. 

-823. 

831. 

-39.45 

385 

1049. 

-822. 

936. 

-38.13 

388 

80. 

-1  10. 

95. 

-4.73 

391 

325. 

-371. 

348. 

-1.61 

394 

1251  . 

-1039. 

1  145. 

-38.51 

397 

1398. 

-1519. 

1  458. 

2.76 

400 

1280. 

-1310. 

1295. 

.17 

403 

1724. 

i   -1618. 

1  671. 

3.55 

406 

172  6. 

-1817. 

1771. 

.38 

409 

1003. 

-1200. 

1  101. 

-42.26 

412 

2006. 

-1400. 

1703. 

8.39 

415 

5186. 

1812. 

1  687. 

-43.81 

418 

"  37. 

-62. 

60. 

5.66 

421 

11. 

-102. 

57. 

7.97 

424 

64. 

-34. 

49. 

-9.22 

427 

-7. 

-23.  • 

8. 

.00 

430 

71  . 

-177. 

124. 

12.07 

433 

59. 

-150. 

105. 

13.28 

436 

57. 

-4  2. 

50. 

-22.50 

439 

165. 

1. 

82. 

-12.67 

442 

658. 

-597. 

627. 

10.57 

448 

1278. 

-1127. 

1202. 

44.07 

451 

1301. 

-1 149. 

1225. 

44.64 

